MATLAB上机实验2.doc

MATLAB上机实验（2） 学号： 姓名： MATLAB程序设计 计算当 n = 100 时,其误差为多少? function F=f(n) F=0; for x=1:n F=F+1/((2*x-1)^2*(2*x+1)^2); End 调用 >> y=f(100) y = 0.1169 >> z=sqrt(16*y+8) z = 3.1416 >> wucha=abs(pi-z) wucha = 5.2263e-008 编写一个分段函数 ，调用并求出x=[-2,0,-2]的值。 function y=fenduan(x) if x<-1 y=1/(1+x); elseif x>=-1 && x<=1 y=1/(1+x.^3); else y=1/(1+x.^2); end 调用 >> a=fenduan(-2) a = -1 >> b=fenduan(0) b = 1 >> c=fenduan(2) c = 0.2000 自己编写一个 MATLAB 函数myhilb(m,n)，使它能自动生成一个m*n的 Hilbert 矩阵Hi,j =1/(i+j-1)，若只给出一个输入参数，则会自动生成一个方阵 。 function A=myhilb(m,n) if nargin ==1 A=zeros(m,m); for i=1:m for j=1:m A(i,j)=1/(i+j+1); end end else A=zeros(m,n); for i=1:m for j=1:n A(i,j)=1/(i+j+1); end end end 使用牛顿迭代法（见绪论，或者查资料）求方程 在区间[0,1]的一个实根，要求精确到4位有效数字。 >> clear >> x=0.0; >> for i=1:20 xl=x; x=x-(2*exp(-x)-sin(x))/(- 2/exp(x) - cos(x)); if abs(x-xl)<0.0001 break x=x end endy = 0.9210