ZKX's LAB

求用MATLAB解多元差分方程的通解或者具体数值 常见差分格式解抛物型方程matlab

2020-07-20知识8

matlab怎么求解偏微分方程 Matlab偏微分方程工具箱应用简介1.概述本文只给出该工具箱的函数列表,读者应先具备偏微分方程的基本知识,然后根据本文列出的函数查阅Matlab的。matlab怎么解偏微分方程,detool是matla的一个重要的工具箱,它可以用数值解法来求解各种繁琐的偏微分方程问题,并且操作非常便捷。它能够画出解的三维图像,更形象具体的。用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*ones(n,1),1)+diag(-beta*ones(n,1),-1);Q=dt/gou/c*ones(n+1,1);边界条件A(1,1)=1;A(1,2)=0;A(end,end)=1;A(end,end-1)=0;T0=25*log(2*pi*x(:));Tseriers=T0;leg_info{1}='t=0';T=T0;i=1;for t=0:dt:tendi=i+1;right=T+Q;边界条件right(1)=0;right(end)=0;T=A\\right;Tseriers=[Tseriers,T];leg_info{i}=['t=',num2str(t)];endplot(x,Tseriers)legend(leg_info)plot(x,T,x,2*exp(-pi*tend/4)*sin(2*pi*x),'r*')legend({['T=',num2str(tend)],'精确解'})

#dx#偏微分方程#差分方程#matlab

随机阅读

qrcode
访问手机版