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

打靶法+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程序).doc
(打靶法,有限差 实验目的 分法) 运用 Matlab/C/C++/Java/Maple/Mathematica 等其中 实验要求 一种语言完成 实验内容 微分方程组边值问题数值算法 (打靶法, ...
打靶法+MATLAB程序.doc
打靶法+MATLAB程序_理学_高等教育_教育专区。求解微分方程的打靶法程序及原
打靶法求边值问题.doc
3.1 常微分方程的解法 ___ 166 3.2 打靶法matlab 实现 ___ 23 致谢:...
...(打靶法,有限差分法)”实验报告(内含matlab程序).doc
数学实验“微分方程组边值问题数值算法(打靶法,有限差分法)”实验报告(内含matlab程序)_数学_自然科学_专业资料。本文档介绍了微分方程组边值问题数值算法(打靶法,...
打靶法求边值问题...doc
3.2 打靶法Matlab 实现 3.2.1 打靶法算法 我们考虑两个边界条
matlab数值计算程序(matlab源代码).pdf
MATLAB/Simulink/C++等编程援助,请访问我的博客 编程博客:http://matlabsite....线性打靶法功能: 求区间[a,b]上的边值问题 x''=p(t)*x'(t)+q(t)*...
二阶非线性常微分方程的打靶法.ppt
综上,实现打靶法对非线性方程的拟合。 ? ? ? 接下来是每一步的代码: ...分别用Bvp4c 函数和打靶法解微分方程: 与matlab自带函数的比较 ? 得到结论:...
Blasius方程的数值解matlab程序.doc
Blasius方程的数值解matlab程序_计算机软件及应用_IT/计算机_专业资料。程序说明该...在用打靶法解题 过程中需要先选定两个初始值(即 f '' 的值)用于叠代计算,...
常微分方程的MATLAB求解_图文.ppt
14.4 一阶微分方程初值问题的数值解 1.欧拉法及其MATLAB实现对于一阶微分方程...14.6 边值问题的数值解 1.打靶法打靶法也称为试射法,其基本思想是把边值...
常微分方程的MATLAB求解_图文.ppt
14.4 一阶微分方程初值问题的数值解 1.欧拉法及其MATLAB实现对于一阶微分方程...14.6 边值问题的数值解 1.打靶法打靶法也称为试射法,其基本思想是把边值...
打靶法的一个c程序.txt
打靶法的一个c程序_电脑基础知识_IT/计算机_专业资料。利用c语言编写的一个...打靶法(含Matlab程序) 10页 1下载券 打靶法+MATLAB程序 2页 1下载券 ...
第14章常微分方程的MATLAB求解_图文.ppt
14.4 一阶微分方程初值问题的数值解 1.欧拉法及其MATLAB实现对于一阶微分方程...14.6 边值问题的数值解 1.打靶法打靶法也称为试射法,其基本思想是把边值...
MATLAB7课件(数学建模培训)-14.7.12-3_图文.ppt
迭加法 ? 打靶法 ? 松弛法 Matlab求解常微分方程初值问题方法 将待求解转化为...z, z2 ? z1 ' 此程序有什么错误? 边值问题的求解 z ?? ? (? ? 2q...
MATLAB在化工中的应用-第5讲常微分方程数值解_图文.ppt
迭加法 ? 打靶法 ? 松弛法 Matlab求解常微分方程初值问题方法 将待求解转化为...z, z2 ? z1 ' 此程序有什么错误? 边值问题的求解 z 求解: ?? ? (?...
实现纳齐谢姆-斯威格特方法的Matlab程序_论文.pdf
实现纳齐谢姆-斯威格特方法的Matlab程序_物理_自然科学_专业资料。流体输运中的...关键词 : 边界层微分方程 ; 打靶 法; 纳齐谢姆 一 斯威 格法 ;tb言 Maa...
Matlab边值问题.pdf
Matlab边值问题_自然科学_专业资料。Matlab边值问题 第2卷第I3期26( 1年 1...on 3 打靶法编程首先任意选取一个: 求解初值方程()()得到足够好的迭代初始值...
MATLAB在化工中的应用-第5讲常微分方程数值解 2_图文.ppt
迭加法 ? 打靶法 ? 松弛法 Matlab求解常微分方程初值问题方法 将待求解转化为...z, z2 ? z1 ' 此程序有什么错误? 边值问题的求解 z 求解: ?? ? (?...
MATLAB7课件(常微分方程数值解)-14.7.12-3_图文.ppt
迭加法 ? 打靶法 ? 松弛法 Matlab求解常微分方程初值问题方法 将待求解转化为...z, z2 ? z1 ' 此程序有什么错误? 边值问题的求解 z ?? ? (? ? 2q...
二阶常微分方程边值问题的数值解法.doc
对线性边值问题,我们总结了 两类常用的数值方法,即打靶法和有限差分方法,对每种方法都列出了详细的计算步骤 和 Matlab 程序代码,通过具体的算例对这两类方法的...
基于MATLAB的固定床煤气化过程动力学模型的模拟研究_李....pdf
MATLAB 软件编写模拟程序求解常微分方程组 ,采 用打靶法对固定床气化炉