上机实验题

附录a上机实验题

实验3

实验6

实验2

可用命令polyval(p,x)进行检验

第二章上机实验题实验1

实验2

实验3

实验6 实验7

第三章上机实验题

果

在命令行中运行>> clear

>> t=linspace(0,2*pi,100);

>> x1=sqrt(5)*cos(t)+2;y1=sqrt(5)*sin(t)+3-2*x1;

>> x2=sqrt(2)*cos(t)+3; y2=6*sin(t);

>> plot(x1,y1,x2,y2)

>> grid on

图形

在命令行中运行显示结果

>>

f=inline('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]','x');

>> [x,y,h]=fsolve(f,[1.7 -2.8])

>> [x,y,h]=fsolve(f,[3.5 -5.8])

>> [x,y,h]=fsolve(f,[4 -4])

>> [x,y,h]=fsolve(f,[1.6 2])

1.7362 -

2.6929

3.4829 -5.6394

4.0287 -4.1171

1.6581

1.8936

实验6

程序ex3_6.m function x=ex3_6(fname,x0,x1,eps,n)

if nargin<5 n=50; end

if nargin<4 eps=1e-6; end

k=0;

while abs(x0-x1)>eps & k<=n

x=x1-feval(fname,x1)/(feval(fname,x1)-feval(fname,x0))*(x1-x0);

x0=x1;x1=x;

k=k+1;

end

实验7

实验8

第四章上机实验题实验1

>> x=-1:3;y=polyval(p,x);

>> nafit(x,y,3)

ans = 7.0000 -1.5714 -4.8571 0.3714

(2)

>> naorthfit(x,y,3)

ans = 7.0000 -1.5714 -4.8571 0.3714

实验2

函数作图>> fplot('1/(1+x^2)',[-5 5]);hold on

五次插值多项式作图>> x=linspace(-5,5,6);y=1./(1+x.^2);xx=linspace(-5,5,100);

>> yy=nalagr(x,y,xx);

>> plot(x,y,'o',xx,yy,'r')

十次插值多项式作图>> x=linspace(-5,5,11);y=1./(1+x.^2);xx=linspace(-5,5,100);

>> yy=nalagr(x,y,xx);

>> plot(x,y,'*',xx,yy,'g')

实验3

程序ex4_3.m function yy=ex4_3(x,y,xx)

% Newton Interpolation

nx=length(x);ny=length(y);

if nx~=ny

warning('The length of vector x and vector y must be

程序ex4_6.m clear;

t=(0:18)*pi/9;

x=[100 134 164 180 198 195 186 160 136 100 66 35 15 0 5 17

32 63 100];

y=[503 525 514.3 451 326.5 188.6 92.2 59.6 62.2 102.7 147.1

191.6 236 280.5 324.9 369.4 413.8 458.3 503];

sx=csape(t,x,'periodic');

tt=(0:0.1:18)*pi/9;

xx=ppval(sx,tt);

sy=csape(t,y,'periodic');

yy=ppval(sy,tt);

plot(x,y,'o',xx,yy)

在命令行中运行>> fun = inline('10-(10-x(1))*exp(-t/x(2))','x','t');

>> t=[0.5 1:5 7 9];

>> v=[6.36 6.48 7.26 8.22 8.66 8.99 9.43 9.63];

>> x=lsqcurvefit(fun,[1 1],t,v)

x = 5.5577 3.5002

注：初始值可选用满足其中任两个数据点的方程组解出

实验8

程序ex4_8.m clear;

t=0:24;

C=[15 14 14 14 14 15 16 18 20 22 23 25 28 31 32 31 29 27 25 24 22 20

18 17 16];

p2=polyfit(t,C,2),e2=norm(C-polyval(p2,t))^2

p3=polyfit(t,C,3),e3=norm(C-polyval(p3,t))^2

p4=polyfit(t,C,4),e4=norm(C-polyval(p4,t))^2

a=lsqcurvefit('ex4_8f',[20 1 15],t,C)

e=norm(C-ex4_8f(a,t))^2

tt=0:0.1:24;

plot(t,C,'o',tt,polyval(p2,tt),tt,polyval(p3,tt),tt,polyval(p4,tt),tt,ex4_8f(a,tt))

legend('数据','二次','三次','四次','函数')

程序ex4_8f.m function C=ex4_8f(a,t)

C=a(1)*exp(-a(2)*(t-a(3)).^2);

图形：

程序ex4_9.m clear;

x=linspace(0,1,31);

y=polyval(ones(1,6),x);

p5=nafit(x,y,5),p15=nafit(x,y,15)

p5=naorthfit(x,y,5),p15=naorthfit(x,y,15)

p5=polyfit(x,y,5),p15=polyfit(x,y,15)

实验10

第五章上机实验题

实验1