matlab中求中心差分的方法？

matlab中求中心差分的方法？~

调用filter函数解差分方程。

1）yn=filter(B,A,xn)是计算输入向量xn的零状态响应输出信号yn，yn与xn长度相等，
其中B=[b0,b1，……bn], A=[a0,a1,……an]。其中a0=1。

2）yn=filter(B,A,xn,xi)是计算全响应的函数。xi是等效初始条件的输入序列，xi能由初始条件确定。此时需要调用filtic函数。xi=filtic(B,A,ys,xs)。其中ys,xs是初始条件向量。
其中ys=[y(-1),y(-2)……y(-N)]，xs=[x(-1),x(-2),……x(-M)]
另外若xn为因果序列xs=0可缺省。

举个例子若y(n)-0.8y(n-1)=x(n)，初始条件y(-1)=1。
a=0.8,ys=1;
xn=[1,zeros(1,30)];
B=1, A=[1,-a];
xi=filtic(B,A,ys);
yn=filter(B,A,xn,xi);
%以下是解出来的yn图像
n=0:length(yn)-1;
stem(n,yn,'.');

调用filter函数解差分方程。

1）yn=filter(B,A,xn)是计算输入向量xn的零状态响应输出信号yn，yn与xn长度相等，
其中B=[b0,b1，……bn], A=[a0,a1,……an]。其中a0=1。

2）yn=filter(B,A,xn,xi)是计算全响应的函数。xi是等效初始条件的输入序列，xi能由初始条件确定。此时需要调用filtic函数。xi=filtic(B,A,ys,xs)。其中ys,xs是初始条件向量。
其中ys=[y(-1),y(-2)……y(-N)]，xs=[x(-1),x(-2),……x(-M)]
另外若xn为因果序列xs=0可缺省。

举个例子若y(n)-0.8y(n-1)=x(n)，初始条件y(-1)=1。
a=0.8,ys=1;
xn=[1,zeros(1,30)];
B=1, A=[1,-a];
xi=filtic(B,A,ys);
yn=filter(B,A,xn,xi);
%以下是解出来的yn图像
n=0:length(yn)-1;
stem(n,yn,'.');

第8章 MATLAB数值积分与微分

8.1 数值积分

8.2 数值微分

<!--[if !supportEmptyParas]--> <!--[endif]-->

8.1 数值积分

8.1.1 数值积分基本原理

求解定积分的数值方法多种多样，如简单的梯形法、辛普生(Simpson) 法、牛顿－柯特斯(Newton-Cotes)法等都是经常采用的方法。它们的基本思想都是将整个积分区间[a,b]分成n个子区间[xi,xi+1]，i=1,2,…,n，其中x1=a，xn+1=b。这样求定积分问题就分解为求和问题。

<!--[if !supportEmptyParas]--> <!--[endif]-->

8.1.2 数值积分的实现方法

1．变步长辛普生法

基于变步长辛普生法，MATLAB给出了quad函数来求定积分。该函数的调用格式为：

[I,n]=quad('fname',a,b,tol,trace)

其中fname是被积函数名。a和b分别是定积分的下限和上限。tol用来控制积分精度，缺省时取tol=0.001。trace控制是否展现积分过程，若取非0则展现积分过程，取0则不展现，缺省时取trace=0。返回参数I即定积分值，n为被积函数的调用次数。

<!--[if !supportEmptyParas]--> <!--[endif]-->

例8-1 求定积分。

(1) 建立被积函数文件fesin.m。

function f=fesin(x)

f=exp(-0.5*x).*sin(x+pi/6);

(2) 调用数值积分函数quad求定积分。

[S,n]=quad('fesin',0,3*pi)

S =

0.9008

n =

77

差分计算方法如下：
一阶差分：gradient命令
二阶差分：del2命令
下面的希望对你有用
http://www.baisi.net/viewthread.php?tid=757131
http://zhseda.bloghome.cn/posts/176364.html

---

正蓝旗19259899841： matlab中求中心差分的方法? [？
汉具头孢： 一阶差分:gradient命令 二阶差分:del2命令L = DEL2(U), when U is a matrix, is a discrete approximation of0.25*del^2 u = (d^2u/dx^2 + d^2/dy^2)/4. The matrix L is the samesize as U, with each element equal to the difference between an...

正蓝旗19259899841： matlab怎样对一个矩阵中的元素用中央差分公式求解后得到另外一个矩阵 - ？
汉具头孢： 你这除以1是要干吗?此外,这里a是等差数列,bi都会是0.代码如下:b = a(1:end-2)+a(3:end)-2*a(2:end-1);

正蓝旗19259899841： 怎么利用MAtlab求离散点的导数 - ？
汉具头孢： 离散点求导都是用差分来进行的,差分一般分前向差分、后向差分和中心差分,形式分别为 ux(i) = u(i)-u(i-1); ux(i) = u(i+1)-u(i); ux(i) = (u(i+1)-u(i-1))/2;

正蓝旗19259899841： matlab计算差分序列 - ？
汉具头孢： 用数值差分指令diff 用法 diff(x')'

正蓝旗19259899841： 数值解法,全隐式有限差分求解,这方程我不会求啊,怎么迭代啊,matlab新手,这周才接触的. - ？
汉具头孢： 老实说,你的那个图片我根本看不清楚,模模糊糊看起来是热传导方程,时间关于向前差分,空间是中心差分,后面好像有些模糊的非齐次部分.然后边界条件看起来是第三种ROBIN边界条件.时间和空间步长都给定了. 这个问题不是你想的那样简单.MATLAB不是万能的.最起码你的方程组要交给MATLAB(不管是线性的还是非线性的).如果你的方程组出来了,那么对于MATLAB就是一个命令或者直接输入人家现成的程序里就有结果出来了. 仅供参考,欢迎批评指正.

正蓝旗19259899841： MATLAB中求导所用的数学函数是说什么?就叫diff函数吗?？
汉具头孢： 1、利用MATLAB库函数diff,这里要注意的是diff求的是差分,即diff求的是dy,要求导数还要除以dx2、利用符号计算工具箱,函数同样是diff你可以参考matlab的帮助文件,或者把正确的函数表达式重新发一下,我帮你算

正蓝旗19259899841： matlab一阶向前差分可以用diff,向后差分用什么?中点公式用什么? - ？
汉具头孢： 向后差分的话,你可以先把你的数组反过来,然后用diff,得到的结果再反一下,数组是横向的话用fliplr反转,数组纵向的话用flipud反转.median是求数组中值的函数.

正蓝旗19259899841： MATLAB 波动方程有限差分方法 - ？
汉具头孢： DX=[]; for x=0:0.001:1 f=1./cos(x).^2; DX=[DX,f]; end

正蓝旗19259899841： matlab图像差分算法,用采集到的主被动图像进行差分运算 - ？
汉具头孢： 图像数据矩阵相减就可以了啊,差分就是减法,把相同的背景区域的值变成0

正蓝旗19259899841： matlab怎么用有限差分法求解方程 - ？
汉具头孢： FungEcho 删除举报此信息 如图中,需要解出力的大小.