应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 使用fsolve求解非线性方程问题
101/1返回列表
查看: 2295  |  回复: 9
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

chnwangnannan

铁虫 (小有名气)

[求助] 使用fsolve求解非线性方程问题

急需求助dbb627版主

e=50;
rho=1.225;
c=2;
s=7.5;
xf=0.48*c;
m=100;
ec=0.23;
aw=2*pi;
M_theradot=-1.2;
I_k=m*s^3*c/3;
I_ktheta=((m*s^2)/2)*(c^2/2-c*xf);
I_theta=m*s*(c^3/3-c^2*xf+c*xf^2);
Kk=I_k*(5*2*pi)^2;
Ktheta=I_theta*(10*2*pi)^2;
b11=rho*c*s^3*aw/6;
b12=0;
b21=-rho*ec*c^2*s^2*aw/4;
b22=-rho*c^3*s*M_theradot/8;
c11=0;
c12=rho*c*s^2*aw/4;
c21=0;
c22=-rho*ec*c^2*s*aw/2;
e11=Kk;
e22=Ktheta+0.75*e*A^2;
a1=I_k*I_theta-I_ktheta^2;
a2=-V^2*(I_theta*c11+I_k*c22-I_ktheta*c12-I_ktheta*c21+b11*b22-b12*b21)-(I_k*e22+I_theta*e11);
a3=V^4*(c11*c22-c12*c21)+V^2*(c11*e22+e11*c22)+e11*e22;
b1=I_ktheta*b12+I_ktheta*b21-I_k*b22-I_theta*b11;
b2=b11*c22+c11*b22-b21*c12-b12*c21;
b3=b11*e22+e11*b22;
F=(b3+b2*V^2)^2*a1-b1*(b3+b2*V^2)*a2+b1^2*a3=0;
当给定参数A=0时，可以求出速度V=154.35。
小于此速度时A的值为复数，大于此值时是实数，我的目的就是当给定不同的速度V时，求解出幅值A的值，因为这个方程是A的四次方的方程，所以用fsolve求解直接找到其实数值，但是我一直无法求解，请指教，谢谢
这个问题的V-A曲线如图所示。

1.JPG
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2012-10-06 21:07:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chnwangnannan

铁虫 (小有名气)
怎么没人回复呢？各位大侠请帮帮忙
赞一下 回复此楼
2楼2012-10-07 16:05:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

感谢参与，应助指数 +1
不知道你要什么结果
F=(b3+b2*V^2)^2*a1-b1*(b3+b2*V^2)*a2+b1^2*a3;
A=solve(F,A)
ezplot(A(1),[140,200])
就可以得到你这个图，如果你要得到数据那么用
for V=140:200
    C(V-139,=subs(A1)';
end
就可以
赞一下 回复此楼
showmethemoney
3楼2012-10-08 11:51:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chnwangnannan

铁虫 (小有名气)
引用回帖:
3楼: Originally posted by csgt0 at 2012-10-08 11:51:35
不知道你要什么结果
F=(b3+b2*V^2)^2*a1-b1*(b3+b2*V^2)*a2+b1^2*a3;
A=solve(F,A)
ezplot(A(1),)
就可以得到你这个图，如果你要得到数据那么用
for V=140:200
    C(V-139,=subs(A1)';
end
就可以...

非常感谢您的回复。
我就想画出上面的那个幅值随速度的变化曲线图。
我用ezplot进行绘制时，出现一下错误：
Error in ==> ezplot>ezplot1 at 420
error(nargchk(4,6,nargin,'struct'));

??? Output argument "hp" (and maybe others) not assigned during call to "C:\Program
Files\MATLAB\R2010a\toolbox\matlab\specgraph\ezplot.m>ezplot1".

Error in ==> ezplot at 144
    [hp,cax] = ezplot1(cax,f{1},vars,labels,args{:});

Error in ==> sym.ezplot at 56
      h = ezplot(char(f),varargin{:});

Error in ==> cal1 at 36
ezplot(A(1),[140,200])

这说明用ezplot不对，请问您画出上面的图了吗？
赞一下 回复此楼
4楼2012-10-08 14:31:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chnwangnannan

铁虫 (小有名气)
引用回帖:
3楼: Originally posted by csgt0 at 2012-10-08 11:51:35
不知道你要什么结果
F=(b3+b2*V^2)^2*a1-b1*(b3+b2*V^2)*a2+b1^2*a3;
A=solve(F,A)
ezplot(A(1),)
就可以得到你这个图，如果你要得到数据那么用
for V=140:200
    C(V-139,=subs(A1)';
end
就可以...

这个ezplot无法绘制出图的原因是不是因为求解的结果A中既有复数又有实数，所填无法绘制？
赞一下 回复此楼
5楼2012-10-08 14:39:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

A=solve(F,A);
ezplot(A(1),[140,200])
4个A第1个和第3个都是部分实数部分复数的可以画
第2个和第4个全部复数好像不能画。当然你可以取实部或虚部或模画。
赞一下 回复此楼
showmethemoney
6楼2012-10-08 15:42:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★
dbb627: 金币+1, 感谢应助 2012-10-10 16:59:28
我能画啊，可能你的版本跟我不用，你可以试下A(2),A(3),A(4)，我的是A(1)和A(3)能画，A(2),A(4)不行
赞一下 回复此楼
showmethemoney
7楼2012-10-08 15:44:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

chnwangnannan

铁虫 (小有名气)
引用回帖:
7楼: Originally posted by csgt0 at 2012-10-08 15:44:06
我能画啊，可能你的版本跟我不用，你可以试下A(2),A(3),A(4)，我的是A(1)和A(3)能画，A(2),A(4)不行

谢谢您的回复，我的目的是，那些为复数的值使其置为0，实数仍然保留
赞一下 回复此楼
8楼2012-10-14 14:49:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

你画了或者输出来看看就知道就知道，4次方程4个解，两个是复数画不了图，另两个实数可以画啊。
算出来A(2)，A(4)全部是复数，A(1),A(3)部分是复数。
复数改为零更好办啊
y =[  1.0000             2.0000             3.0000             4.0000 + 1.0000i]
y(imag(y)~=0)=0

y=1     2     3     0
赞一下 回复此楼
showmethemoney
9楼2012-10-15 09:44:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

csgt0

荣誉版主 (著名写手)

彩色挂图

【答案】应助回帖

★
dbb627: 金币+1, 感谢应助 2012-10-15 20:03:54
如果你觉得画不出来那你看我画的
ezplot(A(1),[140,200])

画A(3)就正负相反。

untitled.jpg
赞一下 回复此楼
showmethemoney
10楼2012-10-15 09:48:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 chnwangnannan 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿985，本科211，0817化学工程与技术319求调剂 +3 Liwangman 2026-03-15 3/150 2026-03-15 18:16 by JourneyLucky
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 308求调剂 +4 是Lupa啊 2026-03-09 4/200 2026-03-14 02:06 by tranquil_ya
[考研] 085600求调剂 +3 a邵星池 2026-03-09 3/150 2026-03-14 01:32 by JourneyLucky
[考研] 考研材料与化工，求调剂 +8 戏精丹丹丹 2026-03-09 8/400 2026-03-14 01:14 by JourneyLucky
[考研] 求调剂，一志愿江南大学环境工程085701 +3 Djdjj12 2026-03-10 4/200 2026-03-14 00:31 by JourneyLucky
[考研] 26考研调剂 +3 ying123. 2026-03-10 3/150 2026-03-14 00:18 by JourneyLucky
[考研] ［0860］321分求调剂，ab区皆可 +4 宝贵热 2026-03-13 4/200 2026-03-13 22:01 by 星空星月
[考研] 材料工程调剂 +4 咪咪空空 2026-03-11 4/200 2026-03-13 19:57 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-12 5/250 2026-03-13 10:56 by houyaoxu
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 求调剂 资源与环境 285 +3 未名考生 2026-03-10 3/150 2026-03-13 10:31 by houyaoxu
[考研] 08食品或轻工求调剂，本科发表3篇sci一区top论文，一志愿南师大食品科学与工程 +3 我是一个兵， 2026-03-10 3/150 2026-03-13 10:21 by Yuyi.
[考研] 工科0856专硕化学工程269能调剂吗 +10 我想读研11 2026-03-10 10/500 2026-03-13 10:14 by Yuyi.
[考研] 321求调剂（食品/专硕） +3 xc321 2026-03-12 6/300 2026-03-13 08:45 by xc321
[考博] 26读博 +4 Rui135246 2026-03-12 10/500 2026-03-13 07:15 by gaobiao
[考博] 26申博求助 +3 跳跃饼干 2026-03-10 4/200 2026-03-10 21:15 by Tntcnn
[考研] 一志愿：武汉理工，材料工程，英二数二 总分314 +3 2202020125 2026-03-10 4/200 2026-03-10 13:54 by xiongyaxuan
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
[考研] 一志愿山东大学，总分327，英语二79，有论文，有竞赛，已过四六级 +3 木木目目1 2026-03-09 3/150 2026-03-09 19:52 by yuningshan
信息提示
关闭
请填处理意见
关闭