Question

非线性二元隐式方程的MATLAB绘图

我的方程式来源于力学里面的模型简化，具体就只有一下两个方程：
F=afa*(deta1+60)*E*sita1/(deta1-2*miu*deta3)+(deta1-deta3)*E*sita1/((3^0.5)*(deta1-2*miu*deta3))
deta1=E*sita1*exp(-(F/7.280388)^0.2314)-13.8768444
其中E=19446，afa=36.45，deta3=30，miu=-0.2312是已知参数
F，deta1和sita1是未知参数，怎么样通过matlab绘制出他的deta1与sita1的图呢？
求大神帮忙，我试过了eplot不得行，感觉是不是要编M文件啊

Answer 1

你把E=19446，afa=36.45，deta3=30，miu=-0.2312代入F

然后把F还有E代入deta1=E*sita1*exp(-(F/7.280388)^0.2314)-13.8768444

得

deta1=19446*sita1*exp(-((36.45*(deta1+60)*19446*sita1/(deta1-2*-0.2312*30)+(deta1-30)*19446*sita1/((3^0.5)*(deta1-2*-0.2312*30)))/7.280388)^0.2314)-13.8768444

再用ezplot

 

Answer 2

Deta=[];
Sita=[];
for deta1=0:0.01:1 %范围及步长根据需要定
for sita1=0:0.01:1
F=afa*(deta1+60)*E*sita1/(deta1-2*miu*deta3)+(deta1-deta3)*E*sita1/((3^0.5)*(deta1-2*miu*deta3));
deta1=E*sita1*exp(-(F/7.280388)^0.2314)-13.8768444;
Deta=[Deta; deta1];
Sita=[Sita; sita1];
end
end
plot(Deta, Sita);