您所在位置：主题：matlab 导数等零的求曲线图
nni120
[专家分：0] 发布于 2010-05-02 13:44:00
matlab 导数等零的求曲线图
syms s M2 p Tc a a1 q Th;
D=1-M2*(p*((Th/Tc)^4-1)+Th/Tc-1);
T=sqrt(((1+a)*Th+a*Tc-a1*D)^2-4*(1+a)*a*Th*Tc);
q=s*D*(a1*D+(a-1)*Th-a*Tc+T)/(2*a*Th+a*(a1*D+(a-1)*Th-a*Tc+T));
qdTh=diff(q,Th)
由上式子 求得 了导数，现在我想求导数等于零的 Th 与 a1 的曲线图
a=0.1;
s=0.8
p=1;
M2=0.0028;
Tc=300;
Th纵坐标 a1 横坐标 谢谢
syms s M2 p Tc a a1 q Th;
D=1-M2*(p*((Th/Tc)^4-1)+Th/Tc-1);
T=sqrt(((1+a)*Th+a*Tc-a1*D)^2-4*(1+a)*a*Th*Tc);
q=s*D*(a1*D+(a-1)*Th-a*Tc+T)/(2*a*Th+a*(a1*D+(a-1)*Th-a*Tc+T));
qdTh=diff(q,Th)
由上式子 求得 了导数，现在我想求导数等于零的 Th 与 a1 的曲线图
a=0.1;
s=0.8
p=1;
M2=0.0028;
Tc=300;
Th纵坐标 a1 横坐标 谢谢