7楼: Originally posted by 815292578 at 2014-03-06 19:41:20
感谢您的回复啊，将=改成==后还是不行。对matlab不熟悉，所以不太会修改。...

gamma=2.856;  
syms psaiI theta1 x psaiII
for psaiI=45.6*pi/1800:45.6*pi/1800:45.6*pi/180
    m1=1+(gamma+1)^2/(gamma*tan(psaiI));
    theta1=atan(tan(psaiI)/(1+gamma*(1+tan(psaiI)*tan(psaiI))));
    temp=solve(tan(theta1)==(m1*m1*sin(x)*cos(x)-cot(x))/(m1*m1*((gamma+1)/2-sin(x)*sin(x))+1), x);  
    psaiII=temp-theta1;
end

9楼: Originally posted by 815292578 at 2014-03-06 21:27:51
感谢这么详细的回复，我用2007和2010的版本怎么运行还是显示错误啊，你是用哪个版本的？
(金币不多，送您红花吧，请见谅)...

function psaiII = t
gamma=2.856;  
syms psaiI theta1 x psaiII
for psaiI=45.6*pi/1800:45.6*pi/1800:45.6*pi/180
    m1=1+(gamma+1)^2/(gamma*tan(psaiI));
    theta1=atan(tan(psaiI)/(1+gamma*(1+tan(psaiI)*tan(psaiI))));
    temp=solve(tan(theta1)==(m1*m1*sin(x)*cos(x)-cot(x))/(m1*m1*((gamma+1)/2-sin(x)*sin(x))+1), x);
    psaiII=temp-theta1;
end

4楼: Originally posted by 沙漠海鱼 at 2014-03-05 02:17:21
你可以在你现在用的matlab版本上help solver！

6楼: Originally posted by hppdyx at 2014-03-06 16:31:27
solve函数中的要使用==来表示等号，所以你把solve括号里面的=改为==试一下

8楼: Originally posted by hppdyx at 2014-03-06 21:01:48
1、程序里面的gamma没有用到
2、solve等号前面最好不要用x（因为你本身的未知量就是x，第一次循环过后就和第二个的x冲突了），我将其改为temp（见后面的代码）
3、solve函数里面等式要用==
4、solve函数里面的等 ...