| 查看: 731 | 回复: 2 | |||
huangshuang1木虫 (著名写手)
|
[求助]
Matlab ode45求解非线性方程 已有1人参与
|
|
function f=mmm(t,y) f=zeros(4,1); f(1)=y(2); f(2)=(50*pi-y(2))/((50*pi-y(2))^2+595.8566)*1.858*10^6-1.52*10^5*(0.0625*y(1)-0.149*y(3))-1.43*10^7*(0.0625*y(1)-0.149*y(3))^2-2.86*10^6*(0.0625*y(1)-0.149*y(3))^3-32.36*(0.0625*y(2)-0.149*y(4)); f(3)=y(4); f(4)=-86.994-2.437*10^4*(0.149*y(3)-0.0625*y(1))-2.29*10^6*(0.149*y(3)-0.0625*y(1))^2-4.58*10^5*(0.149*y(3)-0.0625*y(1))^3-5.185*(0.149*y(4)-0.0625*y(2)); y0=[0,50*pi,0,50*pi]; [t,y]=ode45(@mmm,[0:0.002:20],y0); plot(t,y(:,2),t,y(:,4),':') 用ODE45一直处于busy状态，然后我换成ODE23S，可是结果振荡得很厉害，不收敛。 请问各位，这个该怎么解？？ |
回复此楼» 猜你喜欢
为什么网上很多人说本周 12号出结果
已经有4人回复
-
基金中了
已经有15人回复
-
帮忙看看fileCode
已经有8人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有3人回复
-
FileCode能看出啥？
已经有15人回复
-
filecode
已经有8人回复
-
确定了，国自然21号放榜
已经有7人回复
-
综述论文作为代表作会不会影响评审专家的印象分？
已经有11人回复
-
关于代码变化问题，想知道的进来
已经有24人回复
-
静等基金结果
已经有16人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 求助如何用matlab求非线性方程组的解
已经有8人回复
- 1stopt 语句修改
已经有23人回复
- 求Matlab高手指导
已经有12人回复
- Matlab用fsolve求解非线性方程，卡了三四天了跪求帮助！
已经有35人回复
- MATLAB解非线性方程组
已经有9人回复
- matlab 运行ode45出错，不知道什么原因
已经有4人回复
- 求大神帮忙拟合一个非线性方程，求出模型参数
已经有15人回复
- 用matlab求解方程
已经有6人回复
- 用matlab求解非线性方程组说无解，一定是方程组本身无解，还是有可能程序有问题呢？
已经有11人回复
- matlab如何求解一个非线性微分方程组
已经有8人回复
- 6元方程求解----绘制等高线
已经有5人回复
- MATLAB求解非线性方程组
已经有5人回复
- Matlab用龙格库塔法求二阶非线性方程
已经有13人回复
- 求Matlab高手解决线性方程组的迭代求解问题
已经有17人回复
- matlab求解非线性方程组，错误提示怎么解决
已经有5人回复
- matlab求解方程中的参数
已经有21人回复
- 求助matlab---fsolve解非线性方程组
已经有6人回复
- 1stopt或matlab如何求解以下的非线性方程并拟合出相应曲线？
已经有13人回复
- matlab二阶微分方程求解程序
已经有5人回复
- matlab求解非线性方程组
已经有16人回复
- 求助matlab解积分式
已经有8人回复
- 【求助】向各位大侠求助matlab求解微分方程组遇到的一个问题
已经有21人回复
- 【求助】matlab怎么求解偏微分方程组啊，先谢谢了
已经有13人回复
- 【求助】 求MAtlab求解程序！！！
已经有20人回复
huangshuang1
木虫 (著名写手)
- 应助: 6 (幼儿园)
- 金币: 2762.2
- 帖子: 2263
- 在线: 464小时
- 虫号: 1886296
- 注册: 2012-07-10
- 性别: MM
- 专业: 机械动力学
2楼2015-01-29 14:11:15
CelestialCYJ
木虫 (小有名气)
- 应助: 33 (小学生)
- 金币: 2482.6
- 红花: 12
- 帖子: 144
- 在线: 43.3小时
- 虫号: 3277506
- 注册: 2014-06-16
- 专业: 结构工程
3楼2015-01-31 20:49:30
100