| 查看: 2531 | 回复: 4 | ||
| 【悬赏金币】回答本帖问题，作者白寒冰将赠送您 20 个金币 | ||
[求助]
matlab的fsolve函数求解隐函数绘图，遇到复数解怎么办 已有1人参与
|
||
|
用fsolve求解隐函数绘图，有的显示无解，这个图应该是蝶形的，但是不对，代码如下。可能是因为忽略了复数，因此又对结果进行判断，复数单独绘制，但是还是不对，求大神解答。 %用fsolve绘图 clc clear all; theta0=pi/1000:pi/1000:2*pi; for i=1:length(theta0) theta=theta0(i); C=3; D=(5/36)*pi; a=2; P3=20; E=8; K1=9*(1-E)^2; K2=-12*(1-E)^2+6*(1-E^2)*cos(2*theta); K3=10*(1-E)^2*(cos(2*theta))^2-4*(1-E)^2*(sin(D))^2*(cos(2*theta))^2-2*(1-E)^2*(sin(2*theta))^2-4*(1-E^2)*cos(2*theta)+(1+E)^2; K4=-4*(1-E)^2*cos(4*theta)+2*(1-E^2)*cos(2*theta)-4*(1-E^2)*(sin(D))^2*cos(2*theta)-(4*C*(1-E)*sin(2*D)*cos(2*theta))/P3; K5=(1-E)^2-(sin(D))^2*(1+E+(2*C*cos(D))/(P3*sin(theta)))^2; fun=@(r) K1*(a/r)^8+K2*(a/r)^6+K3*(a/r)^4+K4*(a/r)^2+K5; r_sln=fsolve(fun,pi/10); if isreal(r_sln)==1 r0(i)=r_sln; else r1(i)=abs(r_sln); bdelta = angle(r_sln); end end polar(theta0,r0,'. r'),hold on polar(bdelta,r1,'. y')  12.png  13.png  14.png  微信截图_20210127110318.png |
回复此楼» 猜你喜欢
297，工科调剂?
已经有4人回复
-
恳请有学校收留
已经有7人回复
-
291求调剂
已经有9人回复
-
300求调剂
已经有11人回复
-
22专硕求调剂
已经有12人回复
-
材料相关专业344求调剂双非工科学校或课题组
已经有25人回复
-
急需调剂
已经有7人回复
-
求调剂
已经有10人回复
-
一志愿华中农业071010，320求调剂
已经有16人回复
-
304求调剂
已经有5人回复
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17814.5
- 红花: 53
- 帖子: 2879
- 在线: 606.3小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
2楼2021-01-27 18:31:54
3楼2021-01-28 12:29:00
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17814.5
- 红花: 53
- 帖子: 2879
- 在线: 606.3小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
4楼2021-01-28 12:54:20
5楼2021-01-28 15:05:32
