按关键词阅读: 报告 实验 matlab 微分方程 数值 上机
10、的兴趣 。
偏微分方程数值解法上机报告(四)一、 实验题目:设G是一个十字形区域 , 有五个相等的单位正方形组成 , 用五点差分格式求下列边值问题的数值解:二、 实验目的:通过本次上机 , 掌握椭圆型方程的有限差分法 , 熟悉其计算过程与基本的思想 。
三、 实验代码:h=0.125;
A=zeros(6,14);
for i=1:14A(1,i)=0;
A(6,i)=0;
endfor j=1:6A(j,1)=0;
A(j,14)=0;
endA(2,11)=0;
A(2,12)=0;
A(2,13)=0;
A(3,12)=0;
A(3,13)=0;
A(4,13)=0;
n=0;
for i=2:5for j=2:9+i-1A(I,j)= 。
11、h2/4+(A(i-1,j)+A(I,j-1)+A(i+1,j)+A(I,j+1)/4;
n=n+1;
endendAN四、运行结果:A =Columns 1 through 80 0 0 0 0 0 0 00 0.0039 0.0049 0.0051 0.0052 0.0052 0.0052 0.00520 0.0049 0.0063 0.0068 0.0069 0.0069 0.0069 0.00690 0.0051 0.0068 0.0073 0.0075 0.0075 0.0075 0.00750 0.0052 0.0069 0.0075 0.0076 0.0077 0.0077 0.00770 0 0 0 0 0 0 0Columns 9 through 140 0 0 0 0 00.0052 0.0052 0 0 0 00.0069 0.0069 0.0056 0 0 00.0075 0.0075 0.0072 0.0057 0 00.0077 0.0077 0.0076 0.0072 0.0057 00 0 0 0 0 0N =42五、总结:合理地作网格剖分是处理好边值问题非常重要的一个步骤和前提 , 通过本次上机 , 我懂得了椭圆型方程的有限差分法的基本思想和方法 , 熟悉了其计算过程和步骤 , 并且进一步增强了学习的兴趣 , 提高了编程能力 。
。

稿源:(未知)
【傻大方】网址:/a/2021/0902/0024074160.html
标题:微分方程|偏微分方程数值解上机实验报告(matlab做的)( 二 )