各位朋友，如何用matlab编写这个方程，并求解？ - 仿真模拟 - MATLAB

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

11楼2012-07-31 09:28:16

yan_fei

银虫 (小有名气)

应助: 0 (幼儿园)
金币: 1406.5
帖子: 81
在线: 218.6小时
虫号: 896603
注册: 2009-11-07
专业: 交通工程

引用回帖:

10楼: Originally posted by dbb627 at 2012-07-31 09:27:24
这个确实要估计下t，可以做图估计
function yan
global s phi lam alpha
s=0.5;
phi=9;
lam=0.1389;
A=;
for i=1:11
alpha=A(i) ;
=fsolve(@sumt,16+i*1.9);
end
X
FVAL
t=0.1:0.1:50;
for i=1:le ...

为什么初值给的是16+i*1.9？这点不太理解~~~

12楼2012-07-31 09:44:06

dbb627

荣誉版主 (著名写手)

专家经验: +4
仿真EPI: 6
应助: 289 (大学生)
贵宾: 0.589
金币: 24640.4
散金: 551
红花: 61
沙发: 1
帖子: 1246
在线: 1794.8小时
虫号: 149791
注册: 2005-12-29
性别: GG
专业: 污染控制化学
管辖: 计算模拟

根据图估计的

The more you learn, the more you know, the more you know, and the more you forget. The more you forget, the less you know. So why bother to learn.

13楼2012-08-02 21:18:53

chnwangnannan

铁虫 (小有名气)

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

引用回帖:

13楼: Originally posted by dbb627 at 2012-08-02 21:18:53
根据图估计的

版主您好，我正在求解一个方程（b3+b2*V^2）^2*a1-b1*(b3+b2*V^2)*a2+b1^2*a3=0。其中a2 、a3和b3是关于参数A的表达式，给定不同的V值可以求解得到不同的A值，在解该方程的时候，当A=0时，可求解出V=15.35m/s，小于此速度时A是复数，大于此速度时是实数，用fslove求解可直接接触实数值，但是给定不同速度时，初始值A的设定怎么求，您上面的例子我搞懂了，可是我这个问题弄了好几天了，一直没有求解出结果，希望得到您的帮助

14楼2012-10-06 17:03:01

chnwangnannan

铁虫 (小有名气)

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

引用回帖:

11楼: Originally posted by dbb627 at 2012-07-31 09:28:16
计算结果如图所示。
80/fe/149791_1343698097_782.jpg

版主，你好，我刚才的数据不方便您的查看，现将各个参数列出：
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求解直接找到其实数值，但是我一直无法求解，请指教，谢谢

1.JPG

15楼2012-10-06 17:22:33

gxm1979429

金虫 (正式写手)

应助: 1 (幼儿园)
金币: 162.2
散金: 1013
帖子: 452
在线: 90.1小时
虫号: 815384
注册: 2009-07-25
性别: GG
专业: 自动化

可以到卡尔曼滤波群里提问，63578911

16楼2015-05-30 08:55:31