ZKX's LAB

一维抛物型反问题差分解法 有限差分法的差分方法的发展和应用

2021-03-08知识6

急求!!! 大学数学,用matlab解决问题,题目是一维抛物型偏微分方程差分解法 显式前向欧拉法源程序:function[u,x,t]=EF_Euler(A,xf,T,it0,bx0,bxf,M,N)解方程 A u_xx=u_t,0,0初值: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+1u(i,1)=it0(x(i));endfor j=1:N+1u([1 M+1],j)=[bx0(t(j));bxf(t(j))];endr=A*dt/dx/dx,r1=1-2*r;if(r>;0.5)disp('r>;0.5,unstability');endfor j=1:Nfor i=2:Mu(i,j+1)=r*(u(i+1,j)+u(i-1,j))+r1*u(i,j);(9.2.3)endendu=u';在MATLAB中编写脚本文件:A=0.5;方程系数it0=inline('sin(pi*x)','x');初始条件bx0=inline('0');bxf=inline('0');边界条件xf=2;M=80;T=0.1;N=100;[u1,x,t]=EF_Euler(A,xf,T,it0,bx0,bxf,M,N);figure(1),clf,mesh(u1)xlabel('x')ylabel('t')zlabel('U')title('r>;0.5')M=50;[u1,x,t]=EF_Euler(A,xf,T,it0,bx0,bxf,M,N);figure(2),clf,mesh(u1)xlabel('x')ylabel('t')zlabel('U')title('r)隐式后向欧拉法源程序:function[u,x,t]=IB_Euler(A,xf,T,it0,bx0,bxf,M,N)解方程 A1 u_xx=u_t,0,0初值:u(x,0)=it0(x)边界条件:u(0,t)=bx0(t),u(xf,t)=bxf(t)M:x 轴的。

请求各位大虾,帮忙指点一下,如何用MATLAB解一维抛物型方程?谢谢! 可以用命令pdetool进入pde工具箱来用数值方法求解

一维抛物线法求f(x)=x^2-lnx-5的极值,区间为(0.5,4),用matlab,求程序

请求各位大虾,如何用MATLAB解一维抛物型方程?

有限差分法的差分方法的发展和应用

一维抛物线法求f(x)=x^2-lnx-5的极值,区间为(0.5,4),用matlab,求程序 function hdhsyms x;f=x^2-log(x)-5;a1=0.5;a2=0.8;a3=4;k=0;while(abs(a3-a1)>;1.0e-6)f1=subs(f,x,a1);f2=subs(f,x,a2);f3=subs(f,x,a3);C1=(f3-f1)/(a3-a1);C2=((f2-f1)/(a2-a1)-C1)/(a2-a3);ap=0.5*(a1+a3-C1/C2);fp=subs(f,x,ap);if ap>;a2;if f2>;=fpa1=a2;f1=f2;a2=ap;f2=fp;elsea3=ap;f3=fp;endelseif f2>;=fp;a3=a2;f3=f2;a2=ap;f2=fp;elsea1=ap;f1=fp;endendk=k+1a=apff=subs(f,x,ap)end

一维抛物型反问题差分解法 有限差分法的差分方法的发展和应用

#一维抛物型反问题差分解法

随机阅读

qrcode
访问手机版