清楚的
无图=0;
对于h=0:0.005:0.05
如果(h==0)|(h==0.005)|(h==0.025)|(h==0.05)
标志=h;
虽然标志~=0
m=2.0/h;
对于i=1:m+1
U(i,1)=经验(h*(i-1));
U(1,i)=exp(h*(i-1));
Ua(i,1)=经验(h*((i-1)+(1-1));
Ua(1,i)=exp(h*((1-1)+(i-1));
结束
平均误差=0;
对于i=1:m
对于j=1:m
Ua(i+1,j+1)=经验(h*(i+j));
U(i+1,j+1)=(-U(i,j)+U(i+1,j)+U(i,j+1)+...
((h^2.0)/4.0)*(U(i,j)+U(i+1,j)+U(i,j+1))/(1.0-(h^2.0)/4.0);
误差(i+1,j+1)=abs(U(i+1,j+1)-Ua(i+1,j+1))/abs(Ua(i+1,j+1));
平均误差=误差(i+1,j+1)+平均误差;
结束
结束
标志=0;
结束
标志=h;
虽然标志~=0
打印('\nH=%1.4f',小时)
打印('\nUa(0.25,0.25)=%2.7f,U(0.25,0.25)=%2.7f,错误(0.25,0.25)=%2.7e',...
Ua(圆形((0.25+h)/h)、圆形((0.25+h)/h)、U(圆形((0.25+h)/h)、圆形((0.25+h)/h)),...
误差(圆((0.25+h)/h)、圆((0.25+h)/h)
打印('\nUa(0.5,0.5)=%2.7f,U(0.5,0.5)=%2.7f,错误(0.5,0.5)=%2.7e',...
Ua(圆形((0.5+h)/h)、圆形((0.5+h)/h)、U(圆形((0.5+h)/h)、圆形((0.5+h)/h)),...
误差(圆((0.5+h)/h)、圆((0.5+h)/h)
打印('\nUa(0.75,0.75)=%2.7f,U(0.75,0.75)=%2.7f,错误(0.75,0.75)=%2.7e',...
Ua(圆形((0.75+h)/h)、圆形((0.75+h)/h)、U(圆形((0.75+h)/h)、圆形((0.75+h)/h)),...
误差(圆((0.75+h)/h)、圆((0.75+h)/h)
打印('\nUa(1.0,1.0)=%2.7f,U(1.0,1.0)=%2.7f,错误(1.0,1.0)=%2.7e',...
Ua(圆形((1.0+h)/h)、圆形((1.0+h)/h)、U(圆形((1.0+h)/h)、圆形((1.0+h)/h)),...
误差(圆((1.0+h)/h)、圆((1.0+h)/h)
打印('\nUa(1.25,1.25)=%2.7f,U(1.25,1.25)=%2.7f,错误(1.25,1.25)=%2.7e',...
Ua(圆形((1.25+h)/h)、圆形((1.25+h)/h)、U(圆形((1.25+h)/h)、圆形((1.25+h)/h)),...
误差(圆((1.25+h)/h)、圆((1.25+h)/h)
打印('\nUa(1.5,1.5)=%2.7f,U(1.5,1.5)=%2.7f,错误(1.5,1.5)=%2.7e',...
Ua(圆形((1.5+h)/h)、圆形((1.5+h)/h)、U(圆形((1.5+h)/h)、圆形((1.5+h)/h)),...
误差(圆((1.5+h)/h)、圆((1.5+h)/h)
打印('\nUa(1.75,1.75)=%2.7f,U(1.75,1.75)=%2.7f,错误(1.75,1.75)=%2.7e',...
Ua(圆形((1.75+h)/h)、圆形((1.75+h)/h)、U(圆形((1.75+h)/h)、圆形((1.75+h)/h)),...
误差(圆((1.75+h)/h)、圆((1.75+h)/h)
打印('\nUa(2.0,2.0)=%2.7f,U(2.0,2.0)=%2.7f,错误(2.0,2.0)=%2.7e',...
Ua(圆形((2.0+h)/h)、圆形((2.0+h)/h)、U(圆形((2.0+h)/h)、圆形((2.0+h)/h)),...
误差(圆((2.0+h)/h)、圆((2.0+h)/h)
打印('\n***********************************');
打印('***********************************');
打印('\n');
打印('错误');
打印('对于AM方法(线性)');
打印('\nh:%2.4f',小时);
打印('\n***********************************');
打印('***********************************');
打印('\n | |');
打印(“X”);
打印('\n | Y | 0.5 1.0');
打印('1.5 2.0');
打印('\n | |********************************');
打印('***********************************');
打印('\n | 0.5 |%2.7e%2.7e%2.7e%2.7e',...
误差(圆((0.5+h)/h)、圆((0.5+h)/h)、误差(圆((0.5+h)/h)、圆((1+h)/h)),...
误差(圆((0.5+h)/h)、圆((1.5+h)/h)、误差(圆((0.5+h)/h)、圆((2+h)/h));
打印('\n | 1.0 |%2.7e%2.7e%2.7e%2.7e',...
误差(圆((1+h)/h)、圆((0.5+h)/h)、误差(圆((1+h)/h)、圆((1+h)/h)),...
误差(圆((1+h)/h)、圆((1.5+h)/h)、误差(圆((1+h)/h)、圆((2+h)/h));
打印('\n | 1.5 |%2.7e%2.7e%2.7e%2.7e',...
误差(圆((1.5+h)/h)、圆((0.5+h)/h)、误差(圆((1.5+h)/h)、圆((1+h)/h)),...
误差(圆((1.5+h)/h)、圆((1.5+h)/h)、误差(圆((1.5+h)/h)、圆((2+h)/h));
打印('\n | 2.0 |%2.7e%2.7e%2.7e%2.7e\n',...
误差(圆((2+h)/h)、圆((0.5+h)/h)、误差(圆((2+h)/h)、圆((1+h)/h)),...
误差(圆((2+h)/h)、圆((1.5+h)/h)、误差(圆((2+h)/h)、圆((2+h)/h));
打印('\naverage错误=%2.7e\n\n\n',(平均误差)/(m^2.0));
标志=0;
结束
无图=无图+1;
如果h==0
图(无图)
[X,Y]=网格(0:2/10:2,0:2/10:2);
W=exp(X+Y);
网格(X,Y,W)
埃尔塞夫(h==0.25)|(h==0.5)
图(无图)
[X,Y]=网格(0:2/4:2,0:2/4:2);
W=U(1:m/4:m+1,1:m/4:m+1);
网格(X,Y,W)
其他的
图(无图)
[X,Y]=网格(0:2/10:2,0:2/10:2);
W=U(1:m/10:m+1,1:m/10:m+1);
网格(X,Y,W)
结束
xlabel公司('x轴')
伊拉贝尔('y轴')
兹拉贝尔('z轴')
标题('使用AM方法解决线性问题')
结束
结束
1条评论
直接链接到此评论
//www.tianjin-qmedu.com/matlabcentral/answers/283743如何暂停网格命令#comment_
直接链接到此评论
//www.tianjin-qmedu.com/matlabcentral/answers/283743如何暂停网格命令#comment_
登录发表评论。