matlab习题参考答案(胡良剑.pdf

第一章 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=[1 2 3 4 5];[xx,s]=func1(x) m 2、求满足 的最小 m 值 ln(1 ) +100 n ∑ n 0 s=0; n=0; while(s=100) s=s+log(1+n); n=n+1; end n,s 3、用循环语句形成 Fibonacci 数列 F F F F F k 1, , + 3,4, 。并验证极限 1 2k k −k −1 2 F 1 5+ k → （提示：计算至两边误差小于精度1e-8 为止） F 2 k −1 解： 求 Fibonacci 数列的函数文件： function f=fun(n) if n=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; while abs(x-a)e k=k+1; F(k)=F(k-1)+F(k-2); x=F(k)/F(k-1); end a,x,k 6 10 3 4、分别用 for 和 while 循环结构编写程序，求出K ∑ i ，并考虑一种避免循环语句的程序设