1404的三分之二次方用计算机怎么算,微分第二次大实验.doc

微分第二次大实验

贾梦雪 1100500139

微分方程数值解第二章实验

题目:P129 第二题

网格剖分

进行正方形网格剖分,步长分别为0.2、0.1、0.04, 即J=10、20、50时对应的步长。

建立五点差分格式

uxx =(ui-1,j -2*uij +ui+1,j )/h^2;

uyy= (ui,j-1 -2*uij+ui,j+1 )/h^2;

则五点差分格式为:

-1/h^2*(ui-1,j+ui+1,j-4*uij+ui,j-1+ui,j+1)=

2*pi^2*sin(pi*x)*sin(pi*y)

-1

3.求解方程组

其中A为系数矩阵,f=2*pi^2*sin(pi*x)*sin(pi*y)

求解初值u

4.由PR迭代格式求得迭代b次后的u的值,与初值进行比较

1.J=10时 h=0.2

以下为相关程序

>>f=[];

>> for i=1:9

for j=1:9

x=-0.8:0.2:0.8;

y= -0.8:0.2:0.8;

f(end+1)=2*pi*sin(pi*x(i)).*sin(pi*y(j))

end

end

>> f

>> f

f =

Columns 1 through 7

2.1708 3.5124 3.5124 2.1708 0 -2.1708 -3.5124

Columns 8 through 14

-3.5124 -2.1708 3.5124 5.6832 5.6832 3.5124 0

Columns 15 through 21

-3.5124 -5.6832 -5.6832 -3.5124 3.5124 5.6832 5.6832

Columns 22 through 28

3.5124 0 -3.5124 -5.6832 -5.6832 -3.5124 2.1708

Columns 29 through 35

3.5124 3.5124 2.1708 0 -2.1708 -3.5124 -3.5124

Columns 36 through 42

-2.1708 0 0 0 0 0 0

Columns 43 through 49

0 0 0 -2.1708 -3.5124 -3.5124 -2.1708

Columns 50 through 56

0 2.1708 3.5124 3.5124 2.1708 -3.5124 -5.6832

Columns 57 through 63

-5.6832 -3.5124 0 3.5124 5.6832 5.6832 3.5124

Columns 64 through 70

-3.5124 -5.6832 -5.6832 -3.5124 0 3.5124 5.6832

Columns 71 through 77

5.6832 3.5124 -2.1708 -3.5124 -3.5124 -2.1708 0

Columns 78 through 81

2.1708 3.5124 3.5124 2.1708

(2)

>>for i=1:1:9^2

for j=1:1:9^2

if i==j

a(i,j)=4;

end

if (i==j)

if (mod(i,9)==0)&(mod(j,9)==0)

a(i,j+1)=0;

a(i+1,j)=0;

else

a(i,j+1)=-1;

a(i+1,j)=-1;

end

end

if abs(i-j)==9

a(i,j)=-1;

end

end

end

for i=1:9^2

for j=1:9^2

H(i,j)=a(i,j);

end

end

以上程序求得系数矩阵H (由于H的数据太多,就没有粘贴过来)

(3)

>> u=inv(1/0.04*H)*f' %左式求得的u为uh(0)

此程序则求得方程的81个解

>> u'

ans =

Columns 1 through 9

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值