ririwoking-CSDN博客

原创牛顿插值法 second version

for i=1:10 x(i)=i+1993endy=[67.052,68.008,69.803,72.024,73.400,72.063,74.669,74.487,74.065,74.777]m=2010p(1)=1for i=2:length(y) p(i)=p(i-1)*(m-x(i-1))enddisp(p(2))Q=zeros(length(y))for i=1:length(y) Q(i,1)=y(i)endfor i=2:length(y) .

2021-06-07 09:43:11 89

原创牛顿插值法的应用

function yi=Newton(x,y,xi)n=length(x);m=length(y);if n~=m error('The lengths of X ang Y must be equal!'); return;endY=zeros(n);Y(:,1)=y';for k=1:n-1 for i=1:n-k if abs(x(i+k)-x(i))<eps error('the DATA is error!');.

2021-06-07 09:40:40 1031

原创 MATLAB基础语句

for 循环基本格式for 循环变量=起始值：步长：终止值endwhile 表达式while 表达式循环体End如果使用；可以不显示具体的值如果不加；的话所有的值均会被显示出来Matlab画图>> x=0:pi/100:2*pi;>> y=sin(x);>> plot(x,y)矩阵的提取对角矩阵...

2021-06-07 09:38:42 174

原创 Jacobi迭代解方程

function [x]=jacobi(A,b,x0)%x0表示初始解n=length(b);%知道方程的个数max=500;%最大迭代次数x=zeros(n,1);k=1;eps=1e-3;while k<max for i=1:n x(i)=(b(i)-(A(i,[1:i-1,i+1:n])*x0([1:i-1,i+1:n],1)))/A(i,i); end if norm(x-x0,inf)<eps break; .

2021-06-07 09:36:06 200

原创牛顿法解方程

clear;clc;eps=1e-3;x0=-1;x1=f13(x0);while (abs(x1-x0)>=eps) x0=x1; x1=f13(x1);enddisp(x1);

2021-06-07 09:34:58 100

原创不动点迭代法解方程

clear;clc;eps=1e-3;x0=-1;max=100;x1=f12(x0);for i=1:1:100 if abs(x1-x0)<=eps break; else x0=x1; x1=f12(x1); endendif i==100 disp('达到最大迭代次数');else disp(x1);end

2021-06-07 09:33:33 707

原创二分法解方程

clear;clc;a=-1;b=0;e=1e-3;N = fix(log2((b - a) / e)) + 1;%迭代次数n = 1;x = 0;fa = f1(a);fb = f1(b);if fa * fb > 0 disp('根不在给的区间内');else while n <= N x = (a + b) / 2; if f1(a)*f1(x) > 0 a = x; el...

2021-06-07 09:32:28 62

原创 seidel迭代

function [x]=seidel(A,b,x0)%x0表示初始解n=length(b);%知道方程的个数max=500;%最大迭代次数x=zeros(n,1);k=1;eps=1e-3;while k<max for i=1:n if i==1 x(1)=(b(1)-(A(1,2:n)*x0(2:n,1)))/A(1,1); elseif i==n x(n)=(b(n)-(A(n,1:n-1.

2021-06-07 00:55:52 180

原创 lagrange 插值

function yh=lagrange (x,y,xh)n = length(x);m = length(xh);yh = zeros(1,m); c1 = ones(n-1,1);c2 = ones(1,m);for i=1:n xp = x([1:i-1 i+1:n]); yh = yh + y(i)*prod((c1*xh-xp'*c2)./(x(i)-xp'*c2));endyh=vpa(yh,3);end

2021-06-07 00:53:45 72

ririwoking的博客