| 查看: 1659 | 回复: 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 |
回复此楼» 猜你喜欢
with editor 两个月了，什么原因？
已经有17人回复
-
江苏理工学院要慎选
已经有3人回复
-
博士招生-211农业院校
已经有7人回复
-
材料方向24博士申请/一作SCI三篇
已经有5人回复
-
失之东隅 收之桑榆
已经有6人回复
-
真的好想读博！
已经有23人回复
-
基金开始函评了吗？
已经有17人回复
-
10043以后咋出去啊
已经有4人回复
-
研二光催化6月底四篇二区什么水平
已经有16人回复
-
Nature一直在编辑手里，考虑好几天了，是悬了吗
已经有13人回复
hzlhm
至尊木虫 (著名写手)
- 应助: 373 (硕士)
- 金币: 16185.8
- 红花: 52
- 帖子: 2571
- 在线: 572.3小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
2楼2021-01-27 18:31:54
3楼2021-01-28 12:29:00
hzlhm
至尊木虫 (著名写手)
- 应助: 373 (硕士)
- 金币: 16185.8
- 红花: 52
- 帖子: 2571
- 在线: 572.3小时
- 虫号: 1517335
- 注册: 2011-11-30
- 性别: GG
- 专业: 常微分方程与动力系统
4楼2021-01-28 12:54:20
5楼2021-01-28 15:05:32
5