胡冬冬

2017年3月27日

Matlab:双曲方程

摘要：效果图：阅读全文

posted @ 2017-03-27 18:40 胡冬冬阅读(1554) 评论(0) 推荐(0) 编辑

Matlab:正则Euler分裂

摘要：函数文件1：函数文件2：函数文件3：脚本文件：阅读全文

posted @ 2017-03-27 18:32 胡冬冬阅读(434) 评论(0) 推荐(0) 编辑

2017年3月19日

Mtlab:抛物型方程的交替方向隐格式(ADI)

摘要：取t=1，N=5,10,20,40,80；真解与数值解的结果图： N=5: N=10; N=20; N=40; N=80; 误差阶数：阅读全文

posted @ 2017-03-19 21:55 胡冬冬阅读(1874) 评论(1) 推荐(2) 编辑

2017年3月17日

Matlab:椭圆方程的导数边值问题

摘要： 1 tic; 2 clear 3 clc 4 N=8; 5 M=2*N; 6 h1=2/M; 7 h2=1/N; 8 x=0:h1:2; 9 y=0:h2:1; 10 fun=inline('exp(x)*sin(pi*y)','x','y'); 11 f=inline('(pi^2-1)*exp(x)*sin(pi*y)','x','y'... 阅读全文

posted @ 2017-03-17 22:44 胡冬冬阅读(1583) 评论(0) 推荐(0) 编辑

2017年3月16日

Matlab:Crank Nicolson方法求解线性抛物方程

摘要： 1 tic; 2 clear 3 clc 4 M=[10,20,40,80,160,320,640];%x的步数 5 K=M; %时间t的步数 6 for p=1:length(M) 7 hx=1/M(p); 8 ht=1/K(p); 9 r=ht/hx^2; %网格比 10 x=0:hx:1; 11 t=0:ht:1; 12 numerical=zeros(M(p)+1,K... 阅读全文

posted @ 2017-03-16 17:37 胡冬冬阅读(6495) 评论(0) 推荐(1) 编辑

2017年3月14日

Matlab:非线性高阶常微分方程的几种解法

摘要：一、隐式Euler：函数文件1：函数文件2：函数文件3：脚本文件：效果图：二、变步长的隐式Euler方法：函数文件1：函数文件2：函数文件3：脚本文件：效果图：中心差分法：函数文件1：函数文件2：函数文件3：脚本文件：效果图：阅读全文

posted @ 2017-03-14 20:40 胡冬冬阅读(14256) 评论(0) 推荐(0) 编辑

2017年3月8日

Matlab:高阶常微分三种边界条件的特殊解法（中心差分法，高精度导数边界处理）

摘要：函数文件1：函数文件2：函数文件3：脚本文件: 效果图：阅读全文

posted @ 2017-03-08 21:48 胡冬冬阅读(4662) 评论(7) 推荐(0) 编辑

Matlab:高阶常微分三种边界条件的特殊解法（隐式Euler）

摘要：函数文件1：函数文件2：函数文件3：脚本文件：效果图：阅读全文

posted @ 2017-03-08 21:36 胡冬冬阅读(1299) 评论(0) 推荐(0) 编辑