24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

chnwangnannan

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 75.6
散金: 66
帖子: 103
在线: 10.7小时
虫号: 513955
注册: 2008-02-28
性别: MM

[求助] 使用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人回复
求取一些关于纳米材料和纳米技术相关的英文PPT。已经有0人回复
物理学I论文润色/翻译怎么收费? 已经有94人回复
【复旦大学】二维材料方向招收2026年博士研究生1名已经有0人回复
北京纳米能源与系统研究所王中林院士/曹南颖研究员课题组2026级硕/博/博后招生已经有10人回复
荷兰Utrecht University超快太赫兹光谱王海教授课题招收2026 CSC博士生已经有22人回复
反铁磁体中的磁性切换：两种不同的机制已成功可视化已经有0人回复
26申博推荐：南京航空航天大学国际前沿科学研究院光学方向招收博士生！已经有0人回复
求标准粉末衍射卡号 ICDD 01-076-1802 已经有0人回复

» 本主题相关价值贴推荐，对您同样有帮助:

matlab非线性方程组求解，急！！！已经有5人回复
1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？已经有13人回复
matlab求解非线性方程组，求助！已经有6人回复
求助！matlab用fsolve函数求解非线性方程组的问题！已经有19人回复
如何使用matlab求解非线性方程组的所有整数解？已经有9人回复
matlab求解非线性方程组已经有16人回复
求高人指点用matlab求解非线性方程组，解决了追加100金币；已经有11人回复
matlab的fsove 命令求解非线性方程组已经有6人回复
MATLAB非线性方程组求解~~~反应釜模拟已经有9人回复
【求助】用mathematica 5.0求解一个非线性方程组失败，特发帖求助！已经有5人回复
【求助】matlab求解非线性方程组，并画图处理。要求y,z是实数解！已经有18人回复
【求助】用fsolve求解非线性一元方程已经有5人回复

1楼2012-10-06 21:07:49

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

chnwangnannan

铁虫 (小有名气)

应助: 0 (幼儿园)
金币: 75.6
散金: 66
帖子: 103
在线: 10.7小时
虫号: 513955
注册: 2008-02-28
性别: MM

引用回帖:

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不对，请问您画出上面的图了吗？