MATLAB习题参考答案(胡良剑 .pdfVIP

第一章MATLAB入门

4、求近似解

解：x=-2:0.05:2;y=x.^4-2.^x

两个近似解：y1=f(-0.85)=-0.0328;y2=f(1.250)=0.0630

第二章MATLAB编程与作图

1、设x是数组，求均值和方差

解：函数文件如下：

function[xx,s]=func1(x)

n=length(x);

xx=sum(x)/n;

s=sqrt((sum(x.^2)-n*xx^2)/(n-1));

命令窗口：

x=[12345];[xx,s]=func1(x)

m

∑

2、求满足ln(1+n)100的最小m值

n=0

s=0;

n=0;

while(s=100)

s=s+log(1+n);

n=n+1;

end

n,s

3、用循环语句形成Fibonacci数列F=F=1,F=F+F,k=3,4,

。并验证极限

12kk−1k−2

F1+5

k→（提示：计算至两边误差小于精度1e-8为止）

F2

k−1

解：求Fibonacci数列的函数文件：

functionf=fun(n)

ifn=2

f=1;

else

f=fun(n-1)+fun(n-2);

end

验证极限的函数文件：

function[k,a]=funTest(e)

a=abs(1-(1+sqrt(5))/2);

k=2;

while(ae)

k=k+1;

a=abs(fun(k)/fun(k-1)-(1+sqrt(5))/2);

end

命令行：

[k,a]=funTest(10^-8)

k=

21

a=

9.7719e-009

或者M文件如下：

clear;F(1)=1;F(2)=1;k=2;x=0;

e=1e-8;a=(1+sqrt(5))/2;

whileabs(x-a)e

k=k+1;F(k)=F(k-1)+F(k-2);x=F(k)/F(k-1);

end

a,x,k

1063

4、分别用for和while循环结构编写程序，求出K=∑，并考虑一种避免循环语句的程序设

2i

i=1

计，比较各种算法的运行时间。

解：循环结构：M文件loop.m

k=0;

fori=1:10^6

k=k+sqrt(3)*2^-i;

end

k

非循环结构：M文件nonLoop