请教Matlab求解一元方程式

显示全部楼层 · 2009-1-5 21:01:19

方程式为:
arccos[(asin40°-1)/(a+1)]-arccos[(asin40°-1)/(a-1)]-12.5581°=0
我Matlab编程为：
solve('acos((t*sin(40*pi/180)-1)/(t+1))-acos((t*sin(40*pi/180)-1)/(t-1))-12.5581*pi/180')
我的计算结果为2.1685
而实际结算结果应为4.53635
我哪错了呢？请高手指点，谢谢大家！

千问 · 2009-1-5 21:01:19

你看一下这个图,你就知道为什么了. t=0:0.005:5;y=acos((t*sin(40*pi/180)-1)./(t+1))-acos((t*sin(40*pi/180)-1)./(t-1))-12.5581*pi/180;plot(t,y,t,zeros(size(t)))[0,5]之间有三个解(或两个,那个1不好说,我也没有判断就是了)我直接用你的公式计算,有时出现2.1有时出现4.5可能跟它迭代的初值有关....

千问 · 2009-1-5 21:01:19

t=0:0.005:5; y=acos((t*sin(40*pi/180)-1)./(t+1))-acos((t*sin(40*pi/180)-1)./(t-1))-12.5581*pi/180; plot(t,y,t,zeros(size(t)))...