matlab求二阶微分方程数值解，万分感谢是什么意思。顺便能给下matlab解微分方程数值解

来源：蜘蛛抓取(WebSpider) 时间：2016-10-06 13:50 标签： 万分感谢 英文
常微分方程数值解及其Matlab实现
本文对常微分方程初值问题的数值解法作了进一步探讨,并运用Matlab对其中较难以求解的改进尤拉法和四阶龙格库塔法进行编程实现,程序简洁、直观,求解速度快,方法实用性较强.
襄樊学院数学系
华东交通大学
年，卷(期)
机标分类号
国家自然科学基金
本文读者也读过
互动百科相关词
加载中，请稍候
．客服电话
&&8:00-11:30,13:00-17:00(工作日)如何用matlab求解二阶微分方程，以及程序实例
微分方程的求解主要用到一个dsolve函数，如下面的“Uc1=dsolve('D2y+1000*Dy+10^6*y=0','y(0)=10','Dy(0)=0','t');”，可以看出，函数的第一部分是所要求解的微分方程，其次是初始条件，最后是对自变量的说明。
下面给出的程序实例是用于分析一个最简单零输入的二阶电路。其中C=1uf,L=1H。R是不确定的，他的值的选取将会直接影响到方程解的形式以及最后画出的曲线形状，在此我取R的值分别为00欧姆。
while (R&=3000)
Uc1=dsolve('D2y+1000*Dy+10^6*y=0','y(0)=10','Dy(0)=0','t');
It1=-1*diff(Uc1)*(1e-6);
Ul1=diff(It1);
elseif R==2000
Uc2=dsolve('D2y+2000*Dy+10^6*y=0','y(0)=10','Dy(0)=0','t');
It2=-1*diff(Uc2)*(1e-6);
Ul2=diff(It2);
Uc3=dsolve('D2y+3000*Dy+10^6*y=0','y(0)=10','Dy(0)=0','t');
It3=-1*diff(Uc3)*(1e-6);
Ul3=diff(It3);
%while i&=3
figure (1)
xlabel('t')
ylabel('Uc(t)')
p1=ezplot(Uc1,[0,0.02]);
p2=ezplot(Uc2,[0,0.02]);
p3=ezplot(Uc3,[0,0.02]);
title('电容电压Uc(t)变化曲线')
set(p1,'LineStyle',':');
set(p2,'LineStyle','--');
legend('uc1','uc2','uc3');
set(p1,'Color','r');
set(p2,'Color','k');
legend('uc1','uc2','uc3');
figure (2)
xlabel('t')
ylabel('I(t)')
p4=ezplot(It1,[0,0.02]);
p5=ezplot(It2,[0,0.02]);
p6=ezplot(It3,[0,0.02]);
title('电流I（t）变化曲线')
set(p4,'LineStyle',':');
set(p5,'LineStyle','--');
legend('uc1','uc2','uc3');
set(p4,'Color','r');
set(p5,'Color','k');
legend('It1','It2','It3');
figure (3)
xlabel('t')
ylabel('Ul(t)')
p7=ezplot(Ul1,[0,0.02]);
p8=ezplot(Ul2,[0,0.02]);
p9=ezplot(Ul3,[0,0.02]);
title('电感电压Ul(t)变化曲线')
set(p7,'LineStyle',':');
set(p8,'LineStyle','--');
legend('uc1','uc2','uc3');
set(p7,'Color','r');
set(p8,'Color','k');
legend('Ul1','Ul2','Ul3');
1.有很大帮助，看了之后把问题解决了
2.有一点帮助
3.没有帮助
投票已截止
登录之后投票，你也可以出现在此
已投稿到：
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。【图文】微分方程数值解-matlab_百度文库
两大类热门资源免费畅读
续费一年阅读会员，立省24元！
评价文档：
微分方程数值解-matlab
上传于||文档简介
&&里​面​具​体​介​绍​微​分​方​程​的​M​a​t​l​a​b​求​解
大小：1.60MB
登录百度文库，专享文档复制特权，财富值每天免费拿！
你可能喜欢君，已阅读到文档的结尾了呢~~
二阶常微分方程matlab的数值解和解析解分析总报告 精心收集的各类精品文档，欢迎下载！
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
二阶常微分方程matlab的数值解和解析解分析总报告
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer-4.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口查看: 2524|回复: 10|关注: 0
matlab如何求解结构动力学二阶微分方程组
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
关注者： 2
结构动力学微分方程形式如下图：
2.png (6.87 KB, 下载次数: 0)
10:36 上传
现在的关键问题是方程中的[M],[C],[K],[F]是时变的，而且矩阵阶数比较大，每一时刻的对应[M],[C],[K],[F]可以通过建立一个函数.m文件调用生成，对于这样的问题，这个微分方程组该如何求解？用newmark法的话相当于产生下一时刻点的一个等效的[M],[C],[K],[F]来求解，但是这与本身时变的[M],[C],[K],[F]不符，所以求教，谢谢您的回答！
<h1 style="color:# 麦片财富积分
关注者： 2
本帖最后由 shuidilife 于
18:53 编辑
每一步利用newmarkbeta法求解即
<h1 style="color:# 麦片财富积分
我也遇到这样的问题，能否看看你的程序？谢谢
<h1 style="color:# 麦片财富积分
关注者： 2
我也遇到这样的问题，能否看看你的程序？谢谢
A=zeros(n+1,1);
&&Q1(:,1)=A;
&&Q2(:,1)=A;
&&Q3(:,1)=A;
&&v=196/3.6;
for t=0:0.001:L/v
&&[M,C,K,F]=case21(t,n);
&&dt=0.01;
& & q11=Q1(:,k);
& & q21=Q2(:,k);
& & q31=Q3(:,k);
& & [q1,q2,q3]=Newmark1(F,M,C,K,q11,q21,q31,dt);
& & Q1(:,k+1)=q1;
& & Q2(:,k+1)=q2;
& & Q3(:,k+1)=q3;
& & k=k+1;
其中 [M,C,K,F]=case21(t,n); [q1,q2,q3]=Newmark1(F,M,C,K,q11,q21,q31,dt);是另外两个函数文件，根据每个人的情况不同编写
<h1 style="color:# 麦片财富积分
A=zeros(n+1,1);
非常感谢！
<h1 style="color:# 麦片财富积分
好的，非常感谢
<h1 style="color:# 麦片财富积分
楼主，请问你的问题解决没，能否指导下，我也卡住了，万分感谢，毕业设计要用很着急:)
<h1 style="color:# 麦片财富积分
原来很多人卡在动力学这个问题上了，楼主，我也是在动力学求逆这里一直出问题，看来这是一个大坑那。
<h1 style="color:# 麦片财富积分
楼主，请问你的问题解决没，能否指导下，我也卡住了，万分感谢，毕业设计要用很着急 ...
你解决了没？:'(
<h1 style="color:# 麦片财富积分
请教一下，这个编程思路是怎么样的额，若看到，请加我扣扣，谢谢
站长推荐 /2
Powered by

我要回帖

更多关于 万分感谢 英文 的文章

 

随机推荐