亚洲在线久爱草,狠狠天天香蕉网,天天搞日日干久草,伊人亚洲日本欧美

為了賬號安全,請及時綁定郵箱和手機立即綁定
已解決430363個問題,去搜搜看,總會有你想問的

為什么用solve怎么也弄不出來,要怎么才能弄出來。詳情如下?

為什么用solve怎么也弄不出來,要怎么才能弄出來。詳情如下?

b = tand(70);t = (2*3068)/(2*0.01256*0.00351);d = sqrt(13)*t-(2/3)*t*b;syms p q;p = -20000000:20000000:120000000;a = ((3*q-4*d*b+sqrt(9*(q^2)-24*d*b*q-24*(b^2)*p*q-16*(b^2)*(q^2)))/(4*(b^2)));r = (sqrt((2/(3*q))*(p-a)));solve('sprt(((p-a)^2)+(r^2)*(q^2))+r*(d+a*b) = 0','q'
查看完整描述

1 回答

?
叮當貓咪

TA貢獻1776條經驗 獲得超12個贊

由于給出的方程比較復雜,用solve函數求解,得到是無解,所以應改用vpasolve函數求解。具體實現方法如下:
1、由于vpasolve函數不適應求解數組形式的數據,所以只能逐個求解。
2、應使用for循環語句來完成求解。即
for i=1:8
S=vpasolve(eval(eq(i))) %eq為具體的表達式
end
3、實現代碼:
syms q
%assume(q>0);
b = tand(70);
t = (2*3068)/(2*0.01256*0.00351);
d = sqrt(13)*t-(2/3)*t*b;
p = (-20000000:20000000:120000000);
a = ((3*q-4*d*b+sqrt(9*(q^2)-24*d*b*q-24*(b^2).*p*q-16*(b^2)*(q^2)))/(4*(b^2)));
r = (sqrt((2/(3*q)).*(p-a)));
eq=sqrt(((p-a).^2)+(r.^2)*(q.^2))+r.*(d+a.*b);
for i=1:length(r)
S=vpasolve(eval(eq(i)))
end
4、運行上述代碼,可以得到如下結果。








查看完整回答
反對 回復 2023-04-23
  • 1 回答
  • 0 關注
  • 457 瀏覽
慕課專欄
更多

添加回答

舉報

0/150
提交
取消
微信客服

購課補貼
聯系客服咨詢優惠詳情

幫助反饋 APP下載

慕課網APP
您的移動學習伙伴

公眾號

掃描二維碼
關注慕課網微信公眾號