达芬方程的分叉图

王福昌 2018-4-19 23:11

以达芬振子为例 dafen = @(t,x,flag,force) ; ff=0.23:0.002:0.40;%参数 axis( ); options=odeset('RelTol',1e-7);tt=2*pi/1.2;hold on for j=1:length(ff) =ode45(dafen, , ,options,[],ff(j)); i=2000:100:3000; plot(ff(j),x(i,2),'.') en ...

5301 次阅读|没有评论

Laplace逆变换数值方法的MATLAB程序

热度 2 王福昌 2018-4-19 22:25

从mathworks 网站上下载的Laplace逆变换数值方法的MATLAB程序，先保留这里，以后用着的时候下载。 WEEKS_METHOD_Exchange_Dec15.zip invlaplace.zip Numerical_Inverse_Laplace_Transform.zip

5631 次阅读|3 个评论 热度 2

校园秋日随拍

热度 1 王福昌 2017-10-30 23:02

2139 次阅读|1 个评论 热度 1

涞源白石山手机随拍

王福昌 2017-10-29 23:16

2065 次阅读|没有评论

中国工业与应用数学学会第15届年会青岛20171012-15

王福昌 2017-10-22 10:43

中国工业与应用数学学会第15届年会于2017年10月12-15日在青岛黄海饭店举行，去凑了下热闹，碰到了几个大学同学，与大学同学小聚一下。

1819 次阅读|没有评论

国庆节香山随拍

王福昌 2017-10-6 23:29

1978 次阅读|没有评论

怀柔喇叭沟之秋随拍

王福昌 2017-10-6 23:21

图1 门口 图2 人工种植的菊科植物 &nbs ...

2120 次阅读|没有评论

用MATLAB编写预估-校正法程序求分数阶常微分方程组数值解

王福昌 2017-10-5 10:49

主函数比较 % 根据文献自己编写 q = 0.5; %分数阶阶数 fdefun = @(t,y) ; %一元微分方程 y0 =0;%初值y0为列向量 h = 2^(-6);%步长 tspan = ; = fdewfc(q,fdefun,tspan,y0,h); figure(1) plot(t,y(1,1:end)) ; xlabel('t'); ylabel('y(t)'); hold on plot(t,t.^2-t+0.5,'r-.') ...

10700 次阅读|没有评论

本页有 2 篇博文因作者的隐私设置或未通过审核而隐藏