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

打靶法+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 程序代码,通过具体的算例对这两类方法的...
matlab程序设计例题及答案
matlab程序设计例题及答案_其它_高等教育_教育专区。14道典型的例题,为您解决matlab编程困难1.编写程序:计算 1/3+2/5+3/7+……+10/21 法一: s=0; for ...
matlab部分程序
z 2 ? 0(二 重零点) ,试用 MATLAB 绘出该系统在 0 ~ 1kHz 频率范围内的幅频 特性曲线(要求用频率响应的几何求解方法实现) ,并说明该系统的 作用。 f=0...
Matlab画图程序
Matlab画图程序_理学_高等教育_教育专区。画图时,非常有用奥!Matlab 二维绘图 2010-06-08 20:41 本节介绍 MATLAB 的两种基本绘图功能:二维平面图形和三维立体图形...
MATLAB程序大全_图文
MATLAB程序大全_计算机软件及应用_IT/计算机_专业资料。MATLAB程序大全 ...(i); %使用前趋表反向搜寻,借鉴 Dijsktra 中的松弛方 法 k=k+1; end Mi...
MATLAB应用程序与命令大全
MATLAB应用程序与命令大全_IT/计算机_专业资料。MATLAB应用程序,柱状图、极坐标系...(z,5); %利用递归法 surf(z_nearest); contour(z_nearest); %等值图 ...
常用算法的matlab的源程序
(X,Y,Z)立体散点图 MATLAB 中灰度直方图的显示 MATLAB 图象处理工具箱提供了 imhist 函数来计算和显示图象的直方图,imhist 函数的语 法格式为: imhist(I,n) ...
MatLab练习题-程序设计
MatLab练习题-程序设计_理学_高等教育_教育专区。MatLab 练习题( 4 ) - 程 ...2 -4 -4 -2 2 0 4 *12 ,编写函数文件 SQRT.m 用迭代法求 1 a xn ...
matlab程序大全
用筛选法求素数的基本思想是:要 找出 2~m 之间的全部素数,首先在 2~m 中...2 典型的直流电动机双闭环调速系统的结构图如图所示, 编写 MATLAB 程序 (或 ...
更多相关标签: