实验0.3

最新推荐文章于 2018-01-03 22:47:02 发布

JLUGujw

最新推荐文章于 2018-01-03 22:47:02 发布

阅读量784

点赞数 2

分类专栏：计算方法实验课文章标签：数值计算方法

本文链接：https://blog.csdn.net/GJW20160612/article/details/78778868

版权

计算方法实验课专栏收录该内容

8 篇文章 13 订阅

订阅专栏

吉林大学计算方法实验课实验0.3

代码用octave语言编写。

实验题目：

实验代码及现象：

（1）选择函数，f(x) = sin(x);

（2）同样的步长，比较两种方法。

b = [1,2,3,4,5,6,7,8];%导数1
c = [1,2,3,4,5,6,7,8];%导数2
%f(x)=sin(x)
jd = 10^-5;
for n = 1:8
  b(n) = (sin(n+jd)-sin(n))/jd;
end;
for n = 1:8
  c(n) = (sin(n+jd)-sin(n-jd))/(2*jd);
end;
disp(b);
disp(c);

（3）同样的方法，比较不同步长。

jd=[1,2,3,4,5,6,7,8,9,10,11,12,13,14];
a=[1,2,3,4,5,6,7,8,9,10,11,12,13,14];
b=[1,2,3,4,5,6,7,8,9,10,11,12,13,14];
%f(x)=sin(x)
%取单点值=5
for n = 1:14
  jd(n) = 10^(-n);
  end
%disp(jd);已获得精度值
%选择方法1
for n = 1:14
  a(n) = (sin(5+jd(n))-sin(5))/(jd(n))-cos(5);
  b(n) = (sin(5+jd(n))-sin(5-jd(n)))/(2*jd(n))-cos(5);
end;
disp(a);
plot(a);
hold on;
plot(b);