求助贴。非常感谢!
各位高手:工作中需要编个小程序帮着计算一下设备参数。其中有三个方程,知道用matlab解方程组比较方便。可是水平有限。自已编了,可是解不出来。请高手帮忙看一下,出错在哪里。不胜感谢!clc
format long
syms x y ct
A=input('输入A')
LMAX=input('输入LMAX')
LMIN=input('输入LMIN')
R=input('输入R')
D=input('输入D')
AE=(1-LMIN/LMAX)/log10(R)
CE=-1
for i=1:10
h=i/10
[x y ct]=solve('(y+x*tan(ct*pi/180))/LMAX=h','2*log(10)*x/AE/LMAX-sin(2*ct*pi/180)=0','x*x-15*x+130/pi*10^((h+CE)/AE)*cos(ct*pi/180)=0')
x=vpa(x,4)
y=vpa(y,4)
ct=vpa(ct,4)
end