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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 使用fsolve求解非线性方程问题
101/1返回列表
查看: 2481  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 今天务委会开完了，明天出结果吗 +8 angus9576 2026-08-25 8/400 2026-08-25 16:13 by mike1816
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +8 yuleib84 2026-08-25 9/450 2026-08-25 15:23 by Weed114114
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 18/900 2026-08-25 14:50 by Weed114114
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 今天基金会出结果吗？20260819 +17 kkkl_v 2026-08-19 18/900 2026-08-25 09:41 by windflowerwy
[基金申请] 只有每年这种时候来逛逛小木虫 +25 yaoyewhu2008 2026-08-20 27/1350 2026-08-25 08:01 by Equinoxhua
[基金申请] 明天应该可查了！？ +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 今天放榜吗？ +15 布布和一二 2026-08-19 16/800 2026-08-23 09:55 by 张春生
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
信息提示
关闭
请填处理意见
关闭