首頁猿問為什么用solve怎么也弄不出來，...

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

數(shù)據(jù)分析&挖掘

慕碼人8056858 2023-04-20 19:15:18

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條經(jīng)驗獲得超12個贊

由于給出的方程比較復(fù)雜，用solve函數(shù)求解，得到是無解，所以應(yīng)改用vpasolve函數(shù)求解。具體實現(xiàn)方法如下：
1、由于vpasolve函數(shù)不適應(yīng)求解數(shù)組形式的數(shù)據(jù)，所以只能逐個求解。
2、應(yīng)使用for循環(huán)語句來完成求解。即
for i=1:8
S=vpasolve(eval(eq(i))) %eq為具體的表達式
end
3、實現(xiàn)代碼：
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、運行上述代碼，可以得到如下結(jié)果。