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,M,N,即不在边界上,有(M-2)×(N-2))
-(2uij-ui-1,j-ui+1,j)-(2uij-ui,j-1-ui,j+1)=(pi^2-1)e^(2i/M)sin(pi(j/N))
边界处用边界条件(有M×2+2(N-2)个,保证一起有MN个方程求解MN个未知量)
uij=(pi^2-1)e^(2i/M)sin(pi(j/N))
那么整理就是一个矩阵方程Ax=b的形式,这里还有一个矩阵拉直的问题,因为上面写出来的是一个矩阵方程,而不是线性方程组。一般是按列拉直。
3求解这个方程就OK了,这部才是matlab干的活
不负责任的告诉你可以这么求,x=inv(A)b就OK啦。
因为离散点多了,matlab是解不了大规模的。加速算法就要看结构了
主要问题是列方程,系数矩阵A和右端b都有特殊结果吧,自己写吧
楼上乱回答啊,差分电文格式是指gps输出的gps数据的格式,格式定义了每段数据的开头,结尾,长度,还有哪部分是什么数据。
常用的标准有RTCM23,30,现在用30的比较多。也有其他Gps厂商自己的电文格式。
1、差分又名差分函数或差分运算,差分的结果反映了离散量之间的一种变化,是研究离散数学的一种工具。它将原函数f(x) 映射到f(x+a)-f(x+b) 。差分运算,相应于微分运算,是微积分中重要的一个概念。差分又分为前向差分、向后差分及中心差分三种。
2、差分方程(是一种递推地定义一个序列的方程式:序列的每一项目是定义为前一项的函数。某些简单定义的递推关系式可能会表现出非常复杂的(混沌的)性质,他们属于数学中的非线性分析领域。
扩展资料:
差分方程举例:
dy+ydx=0,y(0)=1 是一个微分方程, x取值[0,1](注:解为y(x)=e^(-x));
要实现微分方程的离散化,可以把x的区间分割为许多小区间 [0,1/n],[1/n,2/n],[(n-1)/n,1]
这样上述微分方程可以离散化为:y((k+1)/n)-y(k/n)+y(k/n)(1/n)=0, k=0,1,2,,n-1 (n 个离散方程组)
利用y(0)=1的条件,以及上面的差分方程,可以计算出 y(k/n) 的近似值了。
差分方程的性质
1、Δk(xn+yn)=Δkxn+Δkyn。
2、Δk(cxn)=cΔkxn。
3、Δkxn=∑(-1)jCjkXn+k-j。
4、数列的通项为n的无限次可导函数,对任意k>=1,存在η,有 Δkxn=f(k)(η)。
参考资料来源:百度百科-差分方程
参考资料来源:百度百科-差分
我对Poision方程得五点差分格式有一些了解,看过余德浩写的《微分方程数值解法》那本书,不过现在可能一时半会儿遍不出来程序,呵呵,发现你的问题太晚了!
等到了周末看有空了看能否帮你编出来!
以上就是关于求解微分方程数值解法问题,五点差分格式求解椭圆型偏微分方程。用MATLAB解决。全部的内容,包括:求解微分方程数值解法问题,五点差分格式求解椭圆型偏微分方程。用MATLAB解决。、RTK差分电文格式是什么意思、什么叫差分,差分方程是啥等相关内容解答,如果想了解更多相关内容,可以关注我们,你们的支持是我们更新的动力!
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)