matlab dae,matlab解DAE遇到的问题
各位同事:
我在解DAE时,用示例1是可以调试出结果的,但是实例2却总是
Error using daeic12 (line 77)
This DAE appears to be of index greater than 1.
Error in ode15s (line 311)
[y,yp,f0,dfdy,nFE,nPD,Jfac] = daeic12(odeFcn,odeArgs,t,ICtype,Mt,y,yp0,f0,...
Error in sunli_Nsolve (line 12)
[t,x]=ode15s(@odefun,[0 10],x0,options);%
各位知道为什么吗?
示例1程序
function dx = rigid(t,x,k,l,m)
dx= zeros(4,1);
dx(1)=-0.2*x(1)+k*x(2)*x(3)+l*0.3*x(1)*x(2);
dx(2)=2*x(1)*x(2)-5*x(2)*x(3)-2*m*x(2)*x(2);
dx(3)=x(1)+x(2)+x(3)-1;
dx(4)=-x(4)-1+x(3);%微分方程组
end
M=[1 0 0 0;0 1 0 0;0 0 0 0;0 0 0 0;];%质量矩阵
options=odeset('mass',M);%对以DAE问题,mass属性必须设置
x0=[0.8;0.1;0.1;0.1];%初值
k=1;l=1;m=1;
[t,x]=ode15s(@rigid,[0 20],x0,options,k,l,m);%这里好像不能使用ode45
figure('numbertitle','off','name','DAE demo—by Matlabsky')
plot(t,x)
legend('x1(t)','x2(t)','x3(t)')
实例2程序:
function dx = odefun(t,x)
dx=zeros(6,1);
l=1.33;k0=19212;k2=0.000560;
c0=180;c1=10600;c2=59830;d1=500;d2=3000;tao=17;
rb=96;Dw=820;ut=76;Dsw=20;hec=2000;
dx(1)=1/c1*((hec-d1)*Dw+(l*x(2)-hec)*Dsw+(d1-l*x(2))*x(6)+k0*rb);%pm
dx(2)=1/c2*((hec-d2)*Dw+(l*x(2)-hec)*Dsw+(d2-l*x(2))*x(6)+k0*rb);%hm
dx(3)=k2*(x(6)*x(5)-x(5)*hec);
dx(4)=x(1)-(x(6)-Dsw)^2/(74554*(x(1)-x(4))^1.3);
dx(5)=(l*x(2)*(x(6)-Dsw)+Dsw*hec)/x(6);
dx(6)=ut*4877.6*x(4)*(-854+13.47*x(4)+0.4*x(5)-0.003*x(4)*x(5));%微分方程组
end
M=[1 0 0 0 0 0;0 1 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0;0 0 0 0 0 0];%质量矩阵
options=odeset('mass',M);%对以DAE问题,mass属性必须设置
x0=[25;2800;950;20;2500;850]
[t,x]=ode15s(@odefun,[0 10],x0,options);%这里好像不能使用ode45
figure('numbertitle','off','name','DAE demo—by Matlabsky')
plot(t,x)
legend('x1(t)','x2(t)','x3(t)')
非常疑惑,基本一样,就加了几个参数
matlab dae,matlab解DAE遇到的问题相关推荐
- 2021-01-13 Matlab求解微分代数方程 (DAE)
Matlab求解微分代数方程 (DAE) 什么是微分代数方程? 微分代数方程是一类微分方程,其中一个或多个因变量导数未出现在方程中.方程中出现的未包含其导数的变量称为代数变量,代数变量的存在意味着不能 ...
- matlab ode45 二阶微分,matlab关于ode45解二阶微分方程的困惑
matlab关于ode45解二阶微分方程的困惑 matlab关于ode45解二阶微分方程的困惑 一个二阶微分方程: y''+y'+y=sin(t) 初始条件为y(0)=5,y'(0)=6. 过程: 先 ...
- 小波基函数构造matlab,五种常见小波基函数及其matlab实现全解.docx
五种常见小波基函数及其matlab实现全解 与标准的傅里叶变换相比,小波分析中使用到的小波函数具有不唯一性,即小波函数 具有多样性.小波分析在工程应用中,一个十分重要的问题就是最优小波基的选择问题,因 ...
- 数字信号处理实验matlab版答案,数字信号处理习题答案及matlab实验详解.pdf
数字信号处理习题答案及matlab实验详解.pdf 第一章 参考答案: 1 (1)2 2 14 ,有理数,所以周期为 14 0 3 3 7 2 2 (2 ) 12 ,无理数,非周期 ...
- 非局部相似性 matlab,非局部均值滤波(NLM)和MATLAB程序详解视频教程保持图像细节...
[内容简介]<非局部均值滤波与应用和MATLAB程序详解视频>共6章28节视频,总学时698分钟,合11.6小时.主要内容包括:非局部均值滤波类算法入门,基于滤波参数自适应的非局部均值滤波 ...
- 完全手册-MATLAB使用详解:基础、开发及工程应用
[书名]完全手册-MATLAB使用详解:基础.开发及工程应用 [作者]董霖 编著 [ISBN]978-7-121-07397-7 [出版社]电子工业出版社 [出版日期]2009年1月 [内容简介] M ...
- 汽车理论matlab编程,汽车理论课后作业matlab编程详解(带注释)[试题学习]
<汽车理论课后作业matlab编程详解(带注释)[试题学习]>由会员分享,可在线阅读,更多相关<汽车理论课后作业matlab编程详解(带注释)[试题学习](11页珍藏版)>请在 ...
- 汽车理论matlab编程,汽车理论课后作业matlab编程详解带注释[10页]
<汽车理论课后作业matlab编程详解带注释[10页]>由会员分享,可在线阅读,更多相关<汽车理论课后作业matlab编程详解带注释[10页](11页珍藏版)>请在读根文库上搜 ...
- 汽车理论课后习题matlab程序,汽车理论课后作业matlab编程详解(带注释).doc
汽车理论课后作业matlab编程详解(带注释) 1.3matlab程序: (1)%驱动力-行驶阻力平衡图%货车相关参数. m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;i ...
最新文章
- (转)Blend操作入门: 别站在门外偷看,快进来吧!
- Delphi的类与继承
- 4、Mysql 主从复制报错[ERROR] [MY-013117] 踩坑
- poj-2115 C Looooops(扩展欧几里得)
- 机器学习之单标签多分类及多标签多分类
- java怎么安装_Windows、Linux、Mac下安装JDK
- 谷歌浏览器如何设置flash访问权限
- JavaScript之jQuery够用即可(each循环、位置偏移量和大小、克隆外层标签、嵌套返回顶部、集体绑定事件)
- 计算机网络第四章思维导图_初级会计实务的第四章所有者权益的思维导图丨初级讲堂...
- linux gdb基本调试及其多线程调试
- ORACLE的程序包1-程序包的基
- Android下基于UVC的UsbCam的开发
- pandas根据现有列新添加一列
- java毕设查重,已经毕业的给毕设查重的人一点建议
- 收藏级干货——Auto CAD历史版本功能大盘点(上)
- 学习要有但行好事,莫问前程的心态
- 10个有效的Linux VPS/服务器性能检测与稳定性测评工具/命令
- 家庭教育:怎样安慰不想上学的人
- prometheus+alertmanager 企业微信告警
- Vue Get请求数据