bvp4c_example.m

来自「matlab7.0从入门到精通(光盘内容)」· M 代码 · 共 13 行

%bvp4c_example.m
%Mathieu方程在边界条件下的解
lambda = 15;                                                 %未知参数猜测解
solinit = bvpinit(linspace(0,pi,10),@Mathieu_initfun,lambda);%求初始值
sol = bvp4c(@bvp_Mathieufun,@Mathieu_bcfun,solinit);         %求解方程
xint = linspace(0,pi);                                       %画图
Sxint = deval(sol,xint);
plot(xint,Sxint(1,:))
axis([0 pi -1 1.1])
title('Mathieu方程在边界条件下的解'); 
xlabel('x')
ylabel('y')
text(0.1,1,['4阶Mathieu方程的特征值 = ' num2str(sol.parameters)]);