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

打靶法+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程序_数学_自然科学_专业资料。Euler,隐式Euler法,改进Euler法Matlab程序1.Euler法 function [x,y]=naeuler(dyfun,xspan,y0,h) x=xspan(1):h...
MATLAB程序大全_图文
MATLAB程序大全_计算机软件及应用_IT/计算机_专业资料。MATLAB程序大全 ...(i); %使用前趋表反向搜寻,借鉴 Dijsktra 中的松弛方 法 k=k+1; end Mi...
数值计算方法matlab程序
列选主元的高斯消元法: 数值计算方法 matlab 程序 function [A,x]=gauss5(A,b) %本算法用列选主元的高斯消元法求解线性方程组 n=length(b); A=[A,b...
常用算法的matlab的源程序
(X,Y,Z)立体散点图 MATLAB 中灰度直方图的显示 MATLAB 图象处理工具箱提供了 imhist 函数来计算和显示图象的直方图,imhist 函数的语 法格式为: imhist(I,n) ...
数值分析matlab程序
数值分析matlab程序_数学_自然科学_专业资料。牛顿迭代法求解方程的根利用牛顿法求非线性方程 9*x^2-sin(x)-1=0 的根。 牛顿迭代公式如下: xk ?1 ? xk ...
matlab编程实例100例
matlab编程实例100例_计算机软件及应用_IT/计算机_专业资料。经典的100例编程,很给力的哦!希望对你们有用!1-32 是:图形应用篇 33-66 是:界面设计篇 67-84 ...
matlab部分程序
z 2 ? 0(二 重零点) ,试用 MATLAB 绘出该系统在 0 ~ 1kHz 频率范围内的幅频 特性曲线(要求用频率响应的几何求解方法实现) ,并说明该系统的 作用。 f=0...
蒙特卡洛方法matlab程序
蒙特卡洛方法 matlab 程序 % Example Monte Carlo Simulation in Matlab % ...matlab 蒙特卡洛法估计... 6页 免费 实验十五: MATLAB的蒙特... 3页 免费...
单纯形法matlab程序
单纯形法matlab程序_自然科学_专业资料。用matlab程序写单纯形法只用最基本的循环,加减乘除运算,没有打包程序的应用便于比较运算量和存储空间 ...
qdpsk程序matlab
Matlab 程序: %---QDPSK---clear all; close all; clc; %---%参数设置 Fc = 5e6; data_len=10000; %---...
更多相关标签: