matlab 用向后差分格式求二维抛物型方程

如题所述

第1个回答  2014-10-19
从田野里跃起的云雀
信心
不丑,也不好看
丁零零宝石和金属的碰撞声
鸟儿几乎都停息了歌声,
的命的美丽哈哈本回答被提问者和网友采纳

用matlab求解抛物型方程,急啊!!用最简隐格式(向后差分格式)求解抛物...
你的精确定绝对有问题。你自己将精确解代入那个泛定方程,或者初值都不符的。一维热传导方程的差分格式 k=1\/16;xleft=0;xright=1;tend=0.2;%时间终值 dx=0.1;dt=0.05;n=(xright-xleft)\/dx;x=xleft:dx:xright;beta=k*dt\/dx\/dx;A=diag((1+2*beta*ones(n+1,1)))+diag(-beta*o...

解二维抛物型方程的高精度显式差分格式
一般也只有二阶精度吧."显式"就是说对时间用"欧拉向前差分"

一维薛定谔方程的差分格式
如果采用向前差分格式,可以得到以下形式:如果采用前后向同时差分,可以得到以下形式:首先,我们使用公式4进行计算。观察公式4,我们可以发现,这是一个显式格式。这意味着,如果我们已知前一时刻各点的值以及空间边界处的值,我们就可以得到下一时刻某点的值。通过这种方式,我们可以依次得到各个时刻的波...

...matlab解决问题,题目是一维抛物型偏微分方程差分解法
解方程 A u_xx = u_t , 0 <= x <= xf, 0 <= t <= T 初值: u(x,0) = it0(x)边界条件: u(0,t) = bx0(t), u(xf,t) = bxf(t)M : x 轴的等分段数 N : t 轴的等分段数 dx = xf\/M; x = [0:M]*dx;dt = T\/N; t = [0:N]'*dt;for i= 1:M ...

MATLAB常用指令(按字母表索引)
删去串尾部的空格符 dec2base 十进制转换为X进制 dec2bin 十进制转换为二进制 dec2hex 十进制转换为十六进制 deconv 多项式除、解卷 delaunay Delaunay 三角剖分 del2 离散Laplacian差分 demo Matlab演示 det 行列式 diag 矩阵对角元素提取、创建对角阵 diary Matlab指令窗文本内容记录 diff 数值差分、符号微分 digits...

求解微分方程数值解法问题,五点差分格式求解椭圆型偏微分方程。用MATLAB...
1.划分网格 把[0,2]×[0,1]的区域划分成M×N个点,即离散的xi=2i\/M,yj=j\/N,i=0,...,M-1,j=0,...,N-1 在网格上的u的值就是uij,就是问题的解 2.离散方程 uxx=2uij-ui-1,j-ui+1,j,uyy=2uij-ui,j-1-ui,j+1(这就是五个点)在第i,j个点的方程为(i,j≠0,...

MATLAB函数库
fft2 二维离散Fourier变换fftn 高维离散Fourier变换fftshift 直流分量对中的谱fieldnames 构架域名figure 创建图形窗fill3 三维多边形填色图find 寻找非零元素下标findobj 寻找具有指定属性的对象图柄findstr 寻找短串的起始字符下标findsym 机器确定内存中的符号变量finverse 符号计算中求反函数fix 向零取整flag 红白蓝黑...

用MATLAB怎么实现离散信号的卷积啊 编程?
实现差分方程,先从简单的说起:filter([1,2],1,[1,2,3,4,5]),实现y[k]=x[k]+2*x[k-1]y[1]=x[1]+2*0=1 (x[1]之前状态都用0)y[2]=x[2]+2*x[1]=2+2*1=4 (2)y=conv(x,h)是用来实现卷级的,对x序列和h序列进行卷积,输出的结果个数等于x的长度与h的长度...

关于matlab filter这个函数到底怎么用?
2、算法实现 filter函数通过直接 II 型转置结构实现,可得下面的时域表达式:y(n)=b(1)*x(n)+b(2)*x(n-1)+...+b(nb+1)*x(n-nb)-a(2)*y(n-1) - ... -a(na+1)*y(n-na)其中 n-1是滤波器阶数,na 是后向反馈的滤波器阶数,nb 是前向反馈的滤波器阶数。

马明书科研成果及奖励
1. **高精度差分格式**:马明书研究了不同维度(一维、二维、三维)的抛物型方程,并开发了相应的高精度差分格式,包括显格式和隐格式。这些格式应用于方程求解,显著提高了计算精度和稳定性。2. **分支稳定性**:在多项研究中,马明书强调了分支稳定性的重要性。他设计的高精度显式差分格式不仅提高...

相似回答