当前位置:首页 >> 理学 >>

打靶法+MATLAB程序


function dy=shootingfun(t,y); % 定义打靶法的微分方程 % y''+ty'-4y=12t^2-3t 0<t<1 % y(0)=0,y(1)=2 dy=[y(2);4*y(1)-t.*y(2)+12*t.^2-3*t];
复制内容到剪贴板 代码:

% 关于线性打靶法的 GUI 文件,函数方程用 shootingfun.m 文件 close all figure; axes('position',[0.06,0.16,0.7,0.6]); p='y''(0)=';P='\Delta(y(1))';Q=[]; Sl=uicontrol(gcf,'style','slider',... 'unit','normalized','position',[0.93,0.02,0.03,0.9],... 'BackgroundColor',[0.9 0.9 0.9],'ForegroundColor','r',... 'fontsize',14,'SliderStep',[0.01,0.01]); set(Sl,'callback',['a=str2num(get(Ed1,''string''));',... 'b=str2num(get(Ed0,''string''));',... 'y0=b+(a-b)*get(Sl,''value'');q=num2str(y0);',... 'set(Te,''string'',[p,q]);eval(SS);']); Te=uicontrol(gcf,'style','text',... 'unit','normalized','position',[0.78,0.92,0.24,0.05],... 'BackgroundColor','w','ForegroundColor','r',... 'fontsize',10,'string',['y''(0)=',... num2str(get(Sl,'value'))]); Ed1=uicontrol(gcf,'style','edit',... 'unit','normalized','position',[0.81,0.86,0.1,0.05],... 'BackgroundColor','w','ForegroundColor','r',... 'string','2'); Ed0=uicontrol(gcf,'style','edit',... 'unit','normalized','position',[0.81,0.02,0.1,0.05],... 'BackgroundColor','w','ForegroundColor','r',... 'string','0'); Te0=uicontrol(gcf,'style','text',... 'unit','normalized','position',[0.11,0.8,0.4,0.2],... 'BackgroundColor',[0.8 0.8 0.9],'ForegroundColor','b',... 'string',{'differential equations:',... 'y''''+ty''-4y=12t^2-3t 0<t<1','y(0)=0,y(1)=2'},... 'fontsize',14,'HorizontalAlignment','left'); plot([0,1],[0,2],'r*');hold on;plot([0,1],[0,2],'rs'); SS=['[t,y]=ode45(@shootingfun,[0,1],[0,y0]);',... 'set(h,''xdata'',t);set(h,''ydata'',y(:,1));',... 'set(DD,''string'',{P,num2str(y(end,1)-2)});'];

h=plot([0,1],[0,0]); DD=uicontrol(gcf,'style','text',... 'unit','normalized','position',[0.55,0.8,0.18,0.14],... 'BackgroundColor','w','ForegroundColor','r',... 'string',{P,Q});


赞助商链接
相关文章:
...(打靶法,有限差分法)”实验报告(内含matlab程序)
数学实验“微分方程组边值问题数值算法(打靶法,有限差分法)”实验报告(内含matlab程序)_数学_自然科学_专业资料。本文档介绍了微分方程组边值问题数值算法(打靶法,...
二阶常微分方程边值问题的数值解法
对线性边值问题,我们总结了 两类常用的数值方法,即打靶法和有限差分方法,对每种方法都列出了详细的计算步骤 和 Matlab 程序代码,通过具体的算例对这两类方法的...
微分方程数值解 打靶法求边值问题
3、实验总结: (实验体会,知识掌握情况,对实验环境的要求等) 本次实验的函数间调用关系较为复杂, 需要仔细分析打靶法的原理与运算的过程才能弄清程序的 编写方法...
数值分析程序目录
数值分析程序目录 - 研究生数值分析matlab源程序全集... 第三篇 数值分析 MATLAB 程序 目录 目录 第一章 ...打靶法及其 MATLAB 程序 10.8.2 有限差分方法及其 ...
二阶常微分方程边值问题的数值解法
对线性边值问题,我们总结了 两类常用的数值方法,即打靶法和有限差分方法,对每种方法都列出了详细的计算步骤 和 Matlab 程序代码,通过具体的算例对这两类方法的...
常微分方程边值问题的数值解法
Matlab 程序 function itaconv=Blasius(n,m,beta,dita,maxita,tol,options) ...定解问题 , 通常主要有两大类求解方法 , 即: 打靶法算法和有限差分方 法....
数学模型期末作业
重庆大学数学模型 数学模型 课程名称 学院 用打靶法求周期解 物理学院 年级 ...四,程序实现:在 MATLAB 中编辑公式如下: function ys=dbf(f,a,b,alfa,beta...
两点边值问题的两种数值解法
介绍常微分方程组两点边值问题的解法--打靶法和配置法常微分方程组两点边值问题...4 y = 0 y ( 0) = 1 y (1) = 3 方法一:配置法 Matlab 程序: ...
2013级数值计算方法课程设计题目
编写 Newton 插值多项式求函 数的近似值的 Matlab 程序,然后解决一个具体实际...20、边值问题的数值解法 设计要求:分析求解边值问题数值解的相关算法(打靶法、...
打靶修正
打靶问题 暂无评价 1页 1财富值 打靶法 5页 1财富值 打靶试题 6页 5财富值...matlab 打靶修正程序matlab 打靶修正程序隐藏>> 对于边值问题 线性边值问题时,p...
更多相关标签: