吉林大学 计算方法实验课 实验0.2
代码用octave语言编写。
题目:
实验代码:
a = [1,2,3,4,5,6,7,8];%e矩阵
r = [1,2,3,4,5,6,7,8];%给的矩阵
m = [1,2,3,4,5,6,7,8];%误差矩阵
rr = [1,2,3,4,5,6,7,8];%自己算的误差矩阵
a(1) = 1- exp(-1);
for n = 2:8
a(n) = 1- n*a(n-1);
%a(n) = round(a(n)*1000)/1000;%不会写函数
end
disp(a);
m(1)=0.00005;
for n = 2:8
m(n) = n*m(n-1);
end
disp(m);
r(1) = 0.632;
for n = 2:8
r(n) = 1-n*r(n-1);
r(n) = round(r(n)*1000)/1000;
end
disp(r);
for n = 1:8
rr(n) = r(n)- a(n);
end
disp(rr);
plot(rr);
%hold on;
%plot(m);
实验现象: