| 查看: 2586 | 回复: 9 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
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 |
回复此楼» 猜你喜欢
澳门大学IAPME李海峰教授常年招收博士研究生
已经有0人回复
-
深圳大学应用技术学院招聘凝聚态物理博士后
已经有123人回复
-
物理学I论文润色/翻译怎么收费?
已经有207人回复
-
求助中医院临床中医中标的青年标书学习，谢谢
已经有5人回复
-
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有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人回复
chnwangnannan
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 75.6
- 散金: 66
- 帖子: 103
- 在线: 10.7小时
- 虫号: 513955
- 注册: 2008-02-28
- 性别: MM
8楼2012-10-14 14:49:31
chnwangnannan
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 75.6
- 散金: 66
- 帖子: 103
- 在线: 10.7小时
- 虫号: 513955
- 注册: 2008-02-28
- 性别: MM
2楼2012-10-07 16:05:45
=subs(A1)';
3楼2012-10-08 11:51:35
chnwangnannan
铁虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 75.6
- 散金: 66
- 帖子: 103
- 在线: 10.7小时
- 虫号: 513955
- 注册: 2008-02-28
- 性别: MM
|
非常感谢您的回复。 我就想画出上面的那个幅值随速度的变化曲线图。 我用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
