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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 使用fsolve求解非线性方程问题
北京石油化工学院2026年研究生招生接收调剂公告
101/1返回列表
查看: 2312  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料专硕(0856) 339分求调剂 +6 哈哈哈鹅哈哈哈 2026-04-04 6/300 2026-04-05 09:53 by 杨光于青云
[考研] 353求调剂 +10 MayUxw1 2026-04-03 10/500 2026-04-05 09:23 by 无际的草原
[考研] 材料与化工371求调剂 +8 陪琳看海 2026-04-04 8/400 2026-04-05 09:10 by guoweigw
[考研] 考研调剂 +5 四川王涛 2026-04-04 5/250 2026-04-04 22:18 by 啵啵啵0119
[考研] 282电子信息0854专硕调剂 +4 202451007219 2026-04-02 6/300 2026-04-04 21:55 by laoshidan
[考研] 材料调剂 +12 一样YWY 2026-04-02 13/650 2026-04-04 20:49 by 蓝云思雨
[考研] 321求调剂 +13 认真求上学 2026-04-02 13/650 2026-04-04 18:23 by macy2011
[考研] 277工科求调剂 +7 1915668 2026-04-04 7/350 2026-04-04 17:21 by 啊俊！
[考研] 319求调剂 +4 星星不眨眼喽 2026-04-03 4/200 2026-04-04 16:25 by 中飞院空管学院�
[考研] 085600材料与化工调剂 +26 kikiki7 2026-03-30 27/1350 2026-04-04 09:18 by qlm5820
[考研] 学硕机械工程303求调剂 +6 无名所以叫吴明 2026-03-30 7/350 2026-04-03 16:48 by asdfzly
[考研] 334求调剂 +9 Trying] 2026-03-31 9/450 2026-04-03 15:18 by 琢珥丶
[考研] 08工科275分求调剂 +14 AaAa7420 2026-03-31 14/700 2026-04-03 11:13 by cocolv
[考研] 专硕 351 086100 也是考的材科基 本科也是材料 +8 202451007219 2026-04-02 8/400 2026-04-03 09:50 by 蓝云思雨
[考研] 一志愿北交大材料工程，总分358 +4 cs0106 2026-04-01 4/200 2026-04-02 07:42 by 尚水阁主
[考博] 26年申博 +3 staryer 2026-03-30 4/200 2026-04-01 23:21 by ai4pharm
[考研] 生物与医药考研调剂 +5 铁憨憨123425 2026-03-31 5/250 2026-04-01 18:01 by syh9288
[考研] 调剂 +3 好好读书。 2026-04-01 3/150 2026-04-01 17:06 by zhouyuwinner
[考研] 求0861交通运输专硕or材料专硕调剂 +4 勒布朗@ 2026-03-31 4/200 2026-04-01 09:54 by 一只好果子?
[考研] 求调剂：一志愿：南京大学 专业：0705 总分320 ，本科985，四六级已过 +3 lfy760306 2026-03-31 3/150 2026-04-01 01:57 by Creta
信息提示
关闭
请填处理意见
关闭