| 查看: 2799 | 回复: 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 |
回复此楼» 猜你喜欢
2026国自然函评费到账
已经有9人回复
-
放榜前的不淡定
已经有14人回复
-
建议基金发布提前给出明确的时间点
已经有13人回复
-
今日不放榜？网传国自然预计 8 月 27 日可查结果
已经有20人回复
-
科研孤儿太难了
已经有18人回复
-
看来今天不会放榜了？
已经有11人回复
-
时间戳又变了
已经有19人回复
-
filecode，4个jtjc了
已经有14人回复
-
让我中一个面上吧！
已经有13人回复
-
人气不行了
已经有7人回复
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17854.5
- 红花: 53
- 帖子: 2880
- 在线: 606.9小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
2楼2021-01-27 18:31:54
3楼2021-01-28 12:29:00
hzlhm
至尊木虫 (著名写手)
- 应助: 387 (硕士)
- 金币: 17854.5
- 红花: 53
- 帖子: 2880
- 在线: 606.9小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
4楼2021-01-28 12:54:20
5楼2021-01-28 15:05:32
40
