前向欧拉公式matlab实现,欧拉法matlab程序.doc
欧拉法matlab程序
1.Euler法function [x,y]=naeuler(dyfun,xspan,y0,h)x=xspan(1):h:xspan(2);y(1)=y0;for n=1:length(x)-1 y(n+1)=y(n)+h*feval(dyfun,x(n),y(n));endx=x';y=y';x1=0:0.2:1;y1=(1+2*x1).^0.5;plot(x,y,x1,y1)
>> dyfun=inline('y-2*x/y');[x,y]=naeuler(dyfun,[0,1],1,0.2);[x,y]
ans =
0 1.0000 0.2000 1.2000 0.4000 1.3733 0.6000 1.5315 0.8000 1.68111.0000 1.8269
2.隐式Euler法function [x,y]=naeulerb(dyfun,xspan,y0,h)x=xspan(1):h:xspan(2);y(1)=y0;for n=1:length(x)-1 y(n+1)=iter(dyfun,x(n+1),y(n),h);endx=x';y=y';x1=0:0.2:1;y1=(1+2*x1).^0.5;plot(x,y,x1,y1)function y=iter(dyfun,x,y,h)y0=y;e=1e-4;K=1e+4;y=y+h*feval(dyfun,x,y);y1=y+2*e;k=1; while abs(y-y1)>e y1=y; y=y0+h*feval(dyfun,x,y); k=k+1; if k>K error('迭代发散'); endend>> dyfun=inline('y-2*x/y');[x,y]=naeulerb(dyfun,[0,1],1,0.2);[x,y]
ans =
0 1.0000 0.2000 1.1641 0.4000 1.3014 0.6000 1.4146 0.8000 1.5019 1.0000 1.5561
3.改进Euler法 function [x,y]=naeuler2(dyfun,xspan,y0,h) x=xspan(1):h:xspan(2); y(1)=y0; for n=1:length(x)-1 k1=feval(dyfun,x(n),y(n)); y(n+1)=y(n)+h*k1; k2=feval(dyfun,x(n+1),y(n+1)); y(n+1)=y(n)+h*(k1+k2)/2; end x=x';y=y'; x1=0:0.2:1;y1=(1+2*x1).^0.5; plot(x,y,x1,y1)
>> dyfun=inline('y-2*x/y');[x,y]=naeuler2(dyfun,[0,1],1,0.2);[x,y]
ans =
0 1.0000 0.2000 1.1867 0.4000 1.3483 0.6000 1.4937 0.8000 1.6279 1.0000 1.7542
前向欧拉公式matlab实现,欧拉法matlab程序.doc相关推荐
- pid算法matlab仿真程序和c程序,pid算法matlab仿真程序和c程序.doc
pid算法matlab仿真程序和c程序.doc 1 增量式 PID 控制算 法 Matlab 仿真程序 设一被控对象 G(s) =50/(0.125s^2+7s),用增量式 PID 控制 算法编写仿真 ...
- dmc matlab程序,matlab编的DMC程序.doc
matlab编的DMC程序 clear all; % close all; %系统模型建立 num=[0.8]; den=[225 1]; [a,b,c,d]=tf2ss(num,den); % st ...
- 微机保护算法matlab编程,微机继电保护matlab算法仿真(有源程序).doc
微机继电保护matlab算法仿真(有源程序) 微 机 继 电 保 护 算 法 仿 真 电控学院 两点乘算法 matlab中编写的程序 N=12; t=(0:0.02/N:0.02)'; m=size( ...
- matlab欧拉法程序,欧拉法matlab程序.doc
欧拉法matlab程序.doc 1.Euler法function[x,y]=naeuler(dyfun,xspan,y0,h)x=xspan(1):h:xspan(2);y(1)=y0;forn=1: ...
- 基于matlab的车牌识别系统程序,基于matlab的车牌识别系统的设计(附程序).doc
基于matlab的车牌识别系统的设计(附程序).doc 1车牌识别系统的设计1.摘要:汽车牌照自动识别系统是制约道路交通智能化的重要因素,包括车牌定位.字符分割和字符识别三个主要部分.本文首先确定车辆 ...
- matlab控制算法C语言,PID算法Matlab仿真程序和C程序
<PID算法Matlab仿真程序和C程序>由会员分享,可在线阅读,更多相关<PID算法Matlab仿真程序和C程序(6页珍藏版)>请在人人文库网上搜索. 1.增量式PID控制算 ...
- 用MATLAB实现转子的故障诊断,单盘转子故障诊断报告(附matlab程序).doc
单盘转子故障诊断报告(附matlab程序).doc 单盘转子故障诊断报告1数据TEST1TXT,数据长度1556486,2组振动数据第一.二列为测点1涡流传感器X.Y振动数据:第三.四列为测点2涡流传 ...
- matlab能否独立做程序,如何将MATLAB程序编译成独立可执行的程序
如何将MATLAB程序编译成独立可执行的程序?如何将编译好的独立可执行程序发布在没有安装MATLAB的电脑上?下面将一步步实现: 一.生成独立可执行的程序(exe文件)步骤 1.安装编译器.可有多种选 ...
- MATLAB工具箱,应用程序,软件和资源的精选清单
精选的MATLAB工具箱,应用程序,软件和资源的精选清单. # Awesome MATLAB [![Awesome](https://cdn.rawgit.com/sindresorhus/aweso ...
最新文章
- RocketMQ之Pull消费者客户端启动
- 科大星云诗社动态20210419
- 电脑小常识----文件名长度过长解决办法
- JavaScript中的Date对象在Safari与IOS中的“大坑”
- U盘专杀工具,U盘防御软件,U盘病毒防火墙--UDiskSyS
- 单片机读tf卡c语言程序,单片机读写U盘闪盘超精简C源程序
- php用什么电脑,我要学php了买一台什么配置的电脑最好?
- 虚拟资源拳王公社:小白从0到1搭建个人私域流量池的实操方法,6招玩转流量裂变法
- python中get啥意思_理解Python中的.get()方法
- php过滤微信特殊字符方案--》2017新版
- RLE压缩及优化--图片压缩
- php pdf添加水印图片,php pdf添加水印(中文水印,图片水印)
- jQuery特效:实现简易轮播图
- python surprise库_Python推荐系统库--Surprise实战
- linux系统官方版下载 百度云,百度网盘linux版
- 斗鱼实名认证 mysql_我告诉你斗鱼实名认证怎么取消
- php 公众号推送图片尺寸,微信公众号发图文消息图片的尺寸是多少为好?
- JavaScript制作日历
- 金蝶系统登录显示云服务器繁忙,金蝶系统登录显示云服务器繁忙
- Delphi TreeView树型框节点拖拽停靠、按条件开始拖拽