您所在位置：主题：Matlab解方程组
wolfnet
[专家分：0] 发布于 2008-03-31 04:41:00
求助,最近小弟在用Matlab解方程组,遇到一个难题,恳请高手指教.
四个未知变量 VRFe1 RFe2 Fluft Ftot
四个方程
VRFe1=sign(Ftot)*l*H/A*abs(Ftot);
RFe2=l*H/A*abs(Fluft)/(abs(Fluft)+10^(-9));
Fluft=Ftot-i/Ro+VRFe1/Ro
[i-VRFe1]*[1/Ro+1/(Rluft+RFe2)]=Ftot;
其中l,H,A,Ro,Rluft,都看作为1.
i如果是常数的话,这个方程能解,但是i是从-400到400递增的,我设i=(-400:1:400),总是出错.
我用的方法是
syms VRFe1 RFe2 Fluft Ftot
[Fluft,Ftot,RFe2,VRFe1]=solve(l*H/A*Ftot-VRFe1,l*H/A*abs(Fluft)/(abs(Fluft)+10^(-9))-RFe2,Ftot-i/Ro+VRFe1/Ro-Fluft,[i-VRFe1]*[1/Ro+1/(Rluft+RFe2)]-Ftot)
请高手指教,最后要的是i和Ftot的图
简单来说就是方程组x-y=i,x+y=0,i是从1,2,3...到10,求i和x的图.
四个未知变量 VRFe1 RFe2 Fluft Ftot
四个方程
VRFe1=sign(Ftot)*l*H/A*abs(Ftot);
RFe2=l*H/A*abs(Fluft)/(abs(Fluft)+10^(-9));
Fluft=Ftot-i/Ro+VRFe1/Ro
[i-VRFe1]*[1/Ro+1/(Rluft+RFe2)]=Ftot;
其中l,H,A,Ro,Rluft,都看作为1.
i如果是常数的话,这个方程能解,但是i是从-400到400递增的,我设i=(-400:1:400),总是出错.
我用的方法是
syms VRFe1 RFe2 Fluft Ftot
[Fluft,Ftot,RFe2,VRFe1]=solve(l*H/A*Ftot-VRFe1,l*H/A*abs(Fluft)/(abs(Fluft)+10^(-9))-RFe2,Ftot-i/Ro+VRFe1/Ro-Fluft,[i-VRFe1]*[1/Ro+1/(Rluft+RFe2)]-Ftot)
请高手指教,最后要的是i和Ftot的图
简单来说就是方程组x-y=i,x+y=0,i是从1,2,3...到10,求i和x的图.