Matlab:五点差分方法求解椭圆方程非导数边值问题

【Matlab:五点差分方法求解椭圆方程非导数边值问题】的更多相关文章

Matlab:五点差分方法求解椭圆方程非导数边值问题

差分格式脚本文件: tic; clear clc M=32;%x的步数 N=16;%y的步数 h1=1/M;%x的步长 h2=1/N;%y的步长 x=0:h1:1; y=0:h2:1; u=zeros(M+1,N+1);%给数值解分配内存单元 U=u;%给精确解分配内存单元 u(1,:)=y.^3;%y边值 u(M+1,:)=1+y.^3;%y边值 u(:,1)=x.^3;%x边值.uo u(:,N+1)=1+x.^3;%x边值.un for i=1:M+1 for j=1:N+1 Accura…

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

tic; clear clc M=[,,,,,,];%x的步数 K=M; %时间t的步数 :length(M) hx=/M(p); ht=/K(p); r=ht/hx^; %网格比 x=:hx:; t=:ht:; numerical=zeros(M(p)+,K(p)+); numerical(:,)=exp(x); %初始值 numerical(,:)=exp(t); %边值 numerical(M(p)+,:)=exp(t+); %边值 a=-r/*ones(M(p)-,);b=(+r)*on…

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

tic; clear clc N=; M=*N; h1=/M; h2=/N; x=:h1:; y=:h2:; fun=inline('exp(x)*sin(pi*y)','x','y'); f=inline('(pi^2-1)*exp(x)*sin(pi*y)','x','y'); lamda1=inline(','y'); lamda2=inline('2*y','y'); lamda3=inline('2*x','x'); lamda4=inline('x^2','x'); kesai1=i…

【matlab】MATLAB程序调试方法和过程

3.8 MATLAB程序的调试和优化在MATLAB的程序调试过程中,不仅要求程序能够满足设计者的设计需求,而且还要求程序调试能够优化程序的性能,这样使得程序调试有时比程序设计更为复杂.MATLAB提供了强大的程序调试功能,合理的运用MATLAB提供的程序调试工具尤其重要.本节从MATLAB程序调试的方法和过程开始介绍,先让用户懂得合理运用MATLAB的程序调试功能,再总结MATLAB程序优化的方法,从而达到实现提高程序性能的目的. 3.8.1 MATLAB程序调试方法和过程(1) MATL…

Complete the sequence! POJ - 1398 差分方法找数列规律

参考链接:http://rchardx.is-programmer.com/posts/16142.html vj题目链接:https://vjudge.net/contest/273000#status/kongbursi/L/0/ 题目给出了一个数列的前若干项,要求推测后面的项.我们很容易想到拉格朗日插值法,但是精度就变成了一个大问题. 这个问题虽然保证了所有的值都是整数,但是并没有保证其多项式的系数也是整数,因此在计算方面存在很大的困难. 除了插值法,求解这种数列问题我们有更好的差分方法,…

FESTUNG — 3. 采用 HDG 方法求解对流问题

FESTUNG - 3. 采用 HDG 方法求解对流问题[1] 1. 控制方程线性对流问题控制方程为 \[\begin{array}{ll} \partial_t c + \nabla \cdot f = h, & \mathrm{in} \; J \times \Omega \\ c(x, 0) = c_0(x), & \mathrm{on} \; {0} \times \Omega \\ c(x, t) = c_D, & \mathrm{on} \; J \times \pa…

彻底理解了call()方法，apply()方法和bind()方法

javascript中的每一个作用域中都有一个this对象,它代表的是调用函数的对象.在全局作用域中,this代表的是全局对象(在web浏览器中指的是window).如果包含this的函数是一个对象的方法,this指向的就是这个对象.因此在上面例子中就不用直接写对象的名字,而是使用this代替它,例如: var human = { name: '霍林林', sayName: function(){ console.log(this.name); } } human.sayName(); 下面这个…