CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 512  |  回复: 2

linan910202

木虫 (正式写手)

齐天大圣

[求助] 纯菜鸟,求助matlab大神帮忙,急用,感谢大神们

已知角sitac的范围,sitac=0:0.017453:0.69812
方程是:beta.*I=1.41421*cos(sitaa).*cos(phic)+sin(sitaa).*O
想求解出sitaa? 如何求解啊?

syms sitaa;
sitac=0:0.017453:0.69812
C=asin(1.41421*sin(sitac/2));
I=2*ellipticE(C,0.5)-ellipticF(C,0.5);
phic=acos(sqrt(sin(sitaa-sitac)));
O=0.8471-2*ellipticE(phic,0.5)+ellipticF(phic,0.5);
beta=sqrt(cos(sitac)./sin(sitaa-sitac));
sitaa=solve('beta.*I=1.41421*cos(sitaa).*cos(phic)+sin(sitaa).*O','sita');

[ Last edited by 月只蓝 on 2017-11-15 at 09:26 ]
回复此楼
蜗牛慢慢爬
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

somomo91

专家顾问 (职业作家)


jjdg: 金币+1, 感谢参与 2017-11-20 17:19:35
这个好像没有解析解,
另外,你程序最后以行,应该改为
sitaa=solve(beta.*I==1.41421*cos(sitaa).*cos(phic)+sin(sitaa).*O, sitaa)
2楼2017-11-19 02:41:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linan910202

木虫 (正式写手)

齐天大圣

引用回帖:
2楼: Originally posted by somomo91 at 2017-11-19 02:41:14
这个好像没有解析解,
另外,你程序最后以行,应该改为
sitaa=solve(beta.*I==1.41421*cos(sitaa).*cos(phic)+sin(sitaa).*O, sitaa)

这位大神,最后一行改完还是不行啊,那请问没有解析解的话,数值解应该怎么解那??
蜗牛慢慢爬
3楼2017-11-20 12:03:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 linan910202 的主题更新
信息提示
请填处理意见