(完整版)第4章MATLAB程序流程控制_习题答案
第4章 MATLAB程序流程控制习题4一、选择题1．下列关于脚本文件和函数文件的描述中不正确的是（）。
AA．函数文件可以在命令行窗口直接运行B．去掉函数文件第一行的定义行可转变成脚本文件C．脚本文件可以调用函数文件D．函数文件中的第一行必须以function开始2．下列程序的输出结果是（）。
Dy=10;if y==10y=20;elseif y>0y=30enddisp(y)A．1 B．30 C．10 D．203．有以下语句：a=eye(5);for n=a(2:end,:)for循环的循环次数是（）。
CA．3 B．4 C．5 D．104．设有程序段k=10;while kk=k-1end则下面描述中正确的是（）。
AA．while循环执行10次B．循环是无限循环C．循环体语句一次也不执行D．循环体语句执行一次5．有以下程序段：x=reshape(1:12,3,4);m=0;n=0;for k=1:4if x(:,k)<=6m=m+1;elsen=n+1;endend则m和n的值分别是（）。
CA．6 6 B．2 1 C．2 2 D．1 26．调用函数时，如果函数文件名与函数名不一致，则使用（）。
A A．函数文件名B．函数名C．函数文件名或函数名均可D．@函数名7．如果有函数声明行为“function [x,y,z]=f1(a,b,c)”，则下述函数调用格式中错误的是（）。
BA．x=f1(a,b,c) B．[x,y,z,w]=f1(a,b,c)C．[x,b,z]=f1(a,y,c) D．[a,b]=f1(x,y,z)8．执行语句“fn=@(x) 10*x;”，则fn是（）。
AA．匿名函数B．函数句柄C．字符串D．普通函数9．执行下列语句后，变量A的值是（）。
D>> f=@(x,y) log(exp(x+y));>> A=f(22,3);A．22,3B．22 C．3 D．2510．程序调试时用于设置断点的函数是（）。
AA．dbstop B．dbclear C．dbcont D．dbstack二、填空题1．将有关MA TLAB命令编成程序存储在一个扩展名为.m的文件中，该文件称为。
M文件2．有语句“for k=[12;34]”引导的循环结构，其循环体执行的次数为。
1 3．MA TLAB中用于控制不确定重复次数的循环语句为，若在循环执行过程中需要终止该循环时采用的语句为。
while…end，break4．函数文件由语句引导。
在函数定义时，函数的输入输出参数称为参数，简称。
在调用函数时，输入输出参数称为参数，简称。
function，形式，形参，实际，实参5．在MA TLAB中，函数文件中的变量是变量。
定义变量是函数间传递信息的一种手段，可以用命令定义。
局部，全局，global6．应用程序的错误有两类，一类是错误，另一类是运行时的错误，即233错误。
MATLAB 程序调试方法有两种，一是利用 进行程序调试，二是利用进行程序调试。
语法，逻辑，调试函数，调试工具 三、应用题1．写出下列程序的输出结果。
s=0;a=[12,13,14;15,16,17;18,19,20;21,22,23]; for k=a for j=1:4if rem(k(j),2)~=0 s=s+k(j); end end end s2．分别用if 语句和switch 语句实现以下计算，其中a 、b 、c 的值从键盘输入。
⎪⎪⎩⎪⎪⎨⎧<≤+<≤+<≤++=5.55.3,ln 5.35.1,sin 5.15.0,2x x c b x x b a x c bx ax y c if 语句:a=input('a='); b=input('b='); c=input('c='); x=input('x='); if 0.5<=x&x<1.5 y=a*x^2+b*x+c; elseif 1.5<=x&x<3.5 y=a*(sin(b))^c+x; elseif 3.5<=x&x<5.5 y=log(abs(b+c/x)); end yswitch 语句： a=input('a='); b=input('b=');c=input('c=');x=input('x=');switch fix(x*10)case num2cell(5:14)y=a*x^2+b*x+c;case num2cell(15:34)y=a*(sin(b))^c+x;case num2cell(35:54)y=log(abs(b+c/x));otherwisedisp('nsrdxcw')endy3．产生20个两位随机整数，输出其中小于平均值的偶数。
x=fix(10+rand(20)*89);a=sum(x)/20;for i=1:20if x(i)<aif rem(x(i),2)==0disp(x(i))endendend4．输入20个数，求其中最大数和最小数。
要求分别用循环结构和调用MA TLAB的max函数、min函数来实现。
循环结构:for a=1:20A(a)=input('请输入20个数：');endAmi=A(1);ma=A(1);for i=1:20if mi>A(i)mi=A(i);endif ma<A(i)ma=A(i);455end end disp(ma) disp(mi) max 、min 函数: for a=1:20A(a)=input('请输入20个数：'); end Ami=min(A),ma=max(A) 5．已知s=1+2+22+23+…+263分别用循环结构和调用MATLAB 的sum 函数求s 的值。
循环结构: s=0; for x=0:63 s=2^x+s; end disp(s) sum 函数: i=0:63; f=2.^i; sum(f)6．当n 分别取100、1000、10000时，求下列各式的值：(1) )2ln (1)1(41312111=+-++-+-+ n n (2) )4(7151311π=+-+-(3) )31(4164116141=+++++n (4) )2()12)(12()2)(2(756653443122πn n n n =⎪⎪⎭⎫⎝⎛+-⎪⎭⎫ ⎝⎛⨯⨯⎪⎭⎫ ⎝⎛⨯⨯⎪⎭⎫⎝⎛⨯⨯ 要求分别用循环结构和向量运算（使用sum 或prod 函数）来实现。
(1): 循环结构:n=input('请输入n 的值：');s=0;for a=1:ns=s+(-1)^(a+1)/a;ends向量运算:n=input('请输入n的值：');a=1:n;f=(-1).^(a+1)./a;s=sum(f)(2):循环结构:n=input('请输入n的值：');s=0;for a=1:n;s=s+(-1)^(a+1)/(2*a-1);ends向量运算:n=input('请输入n的值：');a=1:n;f=(-1).^(a+1)./(2*a-1);s=sum(f)(3):循环结构:n=input('请输入n的值：');s=0;for a=1:ns=s+(1/4)^a;ends向量运算:n=input('请输入n的值：');a=1:n;f=(1/4).^a;s=sum(f)(4):677循环结构:n=input('请输入n 的值：'); s=1; for a=1:nf=(2*a)^2/((2*a-1)*(2*a+1)); s=s*f; end s 向量运算:n=input('请输入n 的值：'); a=1:n;f=(2*a).^2./((2*a-1).*(2*a+1));s=prod(f)7．编写一个函数文件，求小于任意自然数n 的斐波那契（Fibnacci ）数列各项。
斐波那契数列定义如下：121211(2)nn n f f f f fn --=⎧⎪=⎨⎪=+>⎩function f=Fibnacci(n) for i=1:n if i==1 f=1;elseif i==2 f=1; elsef=Fibnacci(i-1)+Fibnacci(i-2); end endn=input('n='); f=Fibnacci(n)8．编写一个函数文件，用于求两个矩阵的乘积和点乘，然后在脚本文件中调用该函数。
function [C,D]=CJ(A,B) C=A*B; D=A.*B;8A=[1 2 3;4 5 6;7 8 9]; B=[1 1 1;1 1 1;1 1 1]; [C,D]=CJ(A,B)9．先用函数的递归调用定义一个函数文件求1nmi i=∑，然后调用该函数文件求∑∑∑===++10501210011k k k kkk 。
function f=xt4_10(m,n) if n<1 f=0; elsef=n^m+xt4_10(m,n-1); endxt4_10(1,100)+xt4_10(2,50)+xt4_10(-1,10) 10．写出下列程序的输出结果。
脚本文件exe.m ： global x x=1:2:5; y=2:2:6; sub(y); x y函数文件sub.m ：function fun=sub(z) global x z=3*x; x=x+z;。
MATLAB课后习题集附标准答案
第2章MATLAB概论1、与其他计算机语言相比较，MA TLAB语言突出的特点是什么？答：起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富，可扩展性强.2、MA TLAB系统由那些部分组成？答：开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装.第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可.4、MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？答：在MATLAB操作桌面上有五个窗口，在每个窗口的右下角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口称为独立的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面称为独立窗口，在独立窗口的view菜单中选择Dock，菜单项就可以将独立的窗口重新防止的桌面上.5、如何启动M文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动.在命令窗口中键入edit命令时也可以启动M文件编辑/调试器.6、存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可.7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中.8、如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file菜单中的Set Path菜单项来完成.在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上.9、在MA TLAB中有几种获得帮助的途径？答：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器.（2）help命令：在命令窗口键入“help”命令可以列出帮助主题，键入“help 函数名”可以得到指定函数的在线帮助信息.（3）lookfor命令：在命令窗口键入“lookfor 关键词”可以搜索出一系列与给定关键词相关的命令和函数.（4）模糊查询：输入命令的前几个字母，然后按Tab键，就可以列出所有以这几个字母开始的命令和函数.注意：lookfor和模糊查询查到的不是详细信息，通常还需要在确定了具体函数名称后用help命令显示详细信息.第3章 MATLAB 数值运算3.1在MA TLAB 中如何建立矩阵⎥⎦⎤⎢⎣⎡194375，并将其赋予变量a ？ 答：在Command Window 窗口输入操作:>> a=[5 7 3;4 9 1]3.2有几种建立矩阵的方法？各有什么优点？ 答：（1）直接输入法，如a=[1 2 3 4]，优点是输入方法方便简捷；（2）通过M 文件建立矩阵，该方法适用于建立尺寸较大的矩阵，并且易于修改； （3）由函数建立，如y=sin(x)，可以由MATLAB 的内部函数建立一些特殊矩阵； （4）通过数据文件建立，该方法可以调用由其他软件产生数据.3.3在进行算术运算时，数组运算和矩阵运算各有什么要求？答：进行数组运算的两个数组必须有相同的尺寸.进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数.渗釤呛俨匀谔鱉调硯錦。
matlab第四章课后作业解答
第四章习题解答1、求下列多项式的所有根，并进行验算。
（3）267235865x x x x-+-（4）4)32(3-+x 解：>> p=zeros(1,24);>> p(1)=5;p(17)=-6;p(18)=8;p(22)=-5;>> root=roots(p)root =0.97680.9388 + 0.2682i0.9388 - 0.2682i0.8554 + 0.5363i0.8554 - 0.5363i0.6615 + 0.8064i0.6615 - 0.8064i0.3516 + 0.9878i0.3516 - 0.9878i-0.0345 + 1.0150i-0.0345 - 1.0150i-0.4609 + 0.9458i-0.4609 - 0.9458i-0.1150 + 0.8340i-0.1150 - 0.8340i-0.7821 + 0.7376i-0.7821 - 0.7376i-0.9859 + 0.4106i-0.9859 - 0.4106i-1.0416-0.7927>> polyval(p,root)ans =1.0e-012 *-0.07120.0459 - 0.0081i0.0459 + 0.0081i-0.0419 + 0.0444i-0.0419 - 0.0444i0.0509 + 0.0929i0.0509 - 0.0929i-0.2059 + 0.0009i-0.2059 - 0.0009i-0.0340 + 0.0145i-0.0340 - 0.0145i0.1342 + 0.0910i0.1342 - 0.0910i0.0025 + 0.0027i0.0025 - 0.0027i-0.0077 + 0.4643i-0.0077 - 0.4643i-0.3548 - 0.1466i-0.3548 + 0.1466i-0.0251-0.0073(4) >> p1=[2 3];>> p=conv(conv(p1,p1),p1)-[0 0 0 4];>> root=roots(p)root =-1.8969 + 0.6874i-1.8969 - 0.6874i-0.7063>> polyval(p,root)ans =1.0e-014 *-0.7105 - 0.6217i-0.7105 + 0.6217i6、求解下列方程组在区域1,0<<βα内的解⎩⎨⎧-=+=.sin 2.0cos 7.0,cos 2.0sin 7.0βαββαα 解：以初值)5.0,5.0(),(00=βα进行求解>> fun=inline('[0.7*sin(x(1))+0.2*cos(x(2))-x(1),0.7*cos(x(1))-0.2*sin(x(2))-x(2)]');>> [x,f,h]=fsolve(fun,[0.5 0.5])Optimization terminated: first-order optimality is less than options.TolFun.x =0.5265 0.5079f =1.0e-007 *-0.1680 -0.2712h =1因而，该方程组的近似根为5079.0,5265.0==βα。
第4章 MATLAB程序流程控制
第4章 MATLAB程序流程控制
2．switch语句 switch语句根据表达式的取值不同，分别执行不 同的语句，其语句格式为：
第4章 MATLAB程序流程控制
(3) 多分支if语句： if 条件1 语句组1 elseif 条件2 语句组2 …… elseif 条件m 语句组m else 语句组n end
第4章 MATLAB程序流程控制
第4章 MATLAB程序流程控制
例4-4 输入一个字符，若为大写字母，则输出其对应 的小写字母；若为小写字母，则输出其对应的大写 字母；若为数字字符则输出其对应的数值，若为其 他字符则原样输出。
age is 17>> >> fprintf('age is %d\n',17); age is 17
第4章 MATLAB程序流程控制
例4-2 求一元二次方程ax2 +bx+c=0的根。 程序如下： a=input('a=?'); b=input('b=?'); c=input('c=?'); delta=b*b-4*a*c; x=[(-b+sqrt(delta))/(2*a),(-b-sqrt(delta))/(2*a)]; disp(['x1=',num2str(x(1)),',x2=',num2str(x(2))]);
M文件的扩展名为.mFra bibliotek第4章 MATLAB程序流程控制
1．建立新的M文件 为建立新的M文件，启动MATLAB编辑器有3种方法。 ① 在MATLAB主窗口选择“主页”选项卡，在“文件” 命令组中单击“新建脚本”命令按钮。 ② 在MATLAB命令行窗口输入命令： edit 文件名 启动MATLAB编辑器后，输入M文件的内容并存盘。 ③ 在命令历史窗口选中一些命令（按住Ctrl可同时选 择多条命令），然后从右键快捷菜单中选择“创建 脚本”命令，将会启动MATLAB编辑器，并在编辑 区中加入所选中的命令。
matlab课后习题答案第四章(内容参考)
符号法在本例中，只求出一个极值点。其余很多极值点无法秋初，更不可能得到最小值。
16
〖目的〗
学习如何把高阶微分方程写成一阶微分方程组。
ode45解算器的导数函数如何采用匿名函数形式构成。
如何从ode45一组数值解点，求指定自变量对应的函数值。
18
〖目的〗
展示特征值分解可能存在的数值问题。
condeig是比较严谨的特征值分解指令。
Jordan分解的作用。
〖解答〗
（1）特征值分解
A=gallery(5)
[V,D]=eig(A);
diag(D)'%为紧凑地显示特征值而写
A =
-9 11 -21 63 -252
70 -69 141 -421 1684
练习mat数据文件中数据的获取。
实验数据求导的后果
把两条曲线绘制在同一图上的一种方法。
〖解答〗
（1）从数据文件获得数据的指令
假如prob_data401.mat文件在当前目录或搜索路径上
clear
load prob_data401.mat
（2）用diff求导的指令
dt=t(2)-t(1);
yc=diff(y)/dt;%注意yc的长度将比y短1
求导会使数据中原有的噪声放大。
12
〖提示〗
指定区间内的积分函数可用cumtrapz指令给出。
在计算要求不太高的地方可用find指令算得。
〖目的〗
指定区间内的积分函数的数值计算法和cumtrapz指令。
find指令的应用。
〖解答〗
dt=1e-4;
t=0:dt:10;
第4章 MATLAB程序流程控制_习题答案
第4章 MATLAB程序流程控制习题4一、选择题1．下列关于脚本文件和函数文件的描述中不正确的是（）。
AA．函数文件可以在命令行窗口直接运行B．去掉函数文件第一行的定义行可转变成脚本文件C．脚本文件可以调用函数文件D．函数文件中的第一行必须以function开始2．下列程序的输出结果是（）。
Dy=10;if y==10y=20;elseif y>0y=30enddisp(y)A．1 B．30 C．10 D．203．有以下语句：a=eye(5);for n=a(2:end,:)for循环的循环次数是（）。
CA．3 B．4 C．5 D．104．设有程序段k=10;while kk=k-1end则下面描述中正确的是（）。
AA．while循环执行10次B．循环是无限循环C．循环体语句一次也不执行D．循环体语句执行一次5．有以下程序段：x=reshape(1:12,3,4);m=0;n=0;for k=1:4if x(:,k)<=6m=m+1;elsen=n+1;endend则m和n的值分别是（）。
CA．6 6 B．2 1 C．2 2 D．1 26．调用函数时，如果函数文件名与函数名不一致，则使用（）。
A A．函数文件名B．函数名C．函数文件名或函数名均可D．@函数名7．如果有函数声明行为“function [x,y,z]=f1(a,b,c)”，则下述函数调用格式中错误的是（）。
BA．x=f1(a,b,c) B．[x,y,z,w]=f1(a,b,c)C．[x,b,z]=f1(a,y,c) D．[a,b]=f1(x,y,z)8．执行语句“fn=@(x) 10*x;”，则fn是（）。
AA．匿名函数B．函数句柄C．字符串D．普通函数9．执行下列语句后，变量A的值是（）。
D>> f=@(x,y) log(exp(x+y));>> A=f(22,3);A．22,3B．22 C．3 D．2510．程序调试时用于设置断点的函数是（）。
matlab第四章课后答案
%Exerc ise 1(1)r oots([1 11])%Exer cise1(2)roots([3 0 -4 0 2 -1])%Exerc ise 1(3)p=zero s(1,24);p([1 17 1822])=[5 -6 8 -5];ro ots(p)%E xerci se 1(4)p1=[2 3];p2=conv(p1,p1);p3=co nv(p1, p2);p3(end)=p3(en d)-4; %原p3最后一个分量-4r oots(p3)%Exer cise2fun=inli ne('x*log(sqrt(x^2-1)+x)-sqrt(x^2-1)-0.5*x');fzer o(fun,2)%Exer cise3fun=inli ne('x^4-2^x');fplot(fun,[-2 2]);gr id on;fze ro(fu n,-1),fzer o(fun,1),f minbn d(fun,0.5,1.5)%Exe rcise 4fu n=inl ine('x*sin(1/x)','x');fp lot(f un, [-0.10.1]);x=z eros(1,10);fori=1:10, x(i)=fz ero(f un,(i-0.5)*0.01);end;x=[x,-x]%Ex ercis e 5f un=in line('[9*x(1)^2+36*x(2)^2+4*x(3)^2-36;x(1)^2-2*x(2)^2-20*x(3);16*x(1)-x(1)^3-2*x(2)^2-16*x(3)^2]','x');[a,b,c]=fso lve(f un,[0 0 0])%E xerci se 6fun=@(x)[x(1)-0.7*si n(x(1))-0.2*cos(x(2)),x(2)-0.7*cos(x(1))+0.2*sin(x(2))]; [a,b,c]=fsolv e(fun,[0.5 0.5])%E xerci se 7clear; clo se; t=0:p i/100:2*pi;x1=2+sqr t(5)*cos(t); y1=3-2*x1+sq rt(5)*sin(t);x2=3+s qrt(2)*cos(t);y2=6*sin(t);pl ot(x1,y1,x2,y2); gri d on; %作图发现4个解的大致位置，然后分别求解y1=fsolv e('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[1.5,2])y2=fsolv e('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[1.8,-2])y3=fsol ve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[3.5,-5])y4=fso lve('[(x(1)-2)^2+(x(2)-3+2*x(1))^2-5,2*(x(1)-3)^2+(x(2)/3)^2-4]',[4,-4])%Exerc ise 8(1)c lear;fun=inlin e('x.^2.*(x.^2-x-2)');fp lot(f un,[-2 2]);grid on;%作图观察x(1)=-2;x(2)=fminb nd(fu n,-1,-0.5);x(4)=fmi nbnd(fun,1,2);fun2=inlin e('-x.^2.*(x.^2-x-2)');x(3)=f minbn d(fun2,-0.5,0.5);x(5)=2feval(fun,x)%答案: 以上x(2)(4)是局部极小，x(1)(3)(5)是局部极大，从最后一句知道x(1)全局最大， x(4)最小。
matlab习题4答案
matlab习题4答案Matlab习题4答案Matlab作为一种强大的数学计算工具，广泛应用于科学研究、工程设计和数据分析等领域。
在学习和使用Matlab过程中，习题是提高技能和理解能力的重要途径。
本文将为大家提供Matlab习题4的答案，帮助读者更好地掌握Matlab 的应用。
习题1：编写一个Matlab程序，计算并输出1到100之间所有正整数的和。
解答：```matlabsum = 0;for i = 1:100sum = sum + i;enddisp(sum);```该程序使用了for循环结构，从1到100遍历所有正整数，并将它们累加到变量sum中。
最后使用disp函数输出sum的值。
习题2：编写一个Matlab程序，计算并输出1到100之间所有奇数的和。
解答：```matlabsum = 0;for i = 1:2:100sum = sum + i;enddisp(sum);```该程序同样使用了for循环结构，但是在每次迭代中，i的增量为2，即只遍历奇数。
其余部分与习题1的解答相同。
习题3：编写一个Matlab程序，生成一个3×3的随机矩阵，并输出该矩阵的转置。
解答：```matlabA = rand(3);disp(A);disp(A');```该程序使用了rand函数生成一个3×3的随机矩阵A，并使用disp函数分别输出A和其转置A'。
习题4：编写一个Matlab程序，计算并输出斐波那契数列的前20个数。
解答：```matlabfib = zeros(1, 20);fib(1) = 1;fib(2) = 1;for i = 3:20fib(i) = fib(i-1) + fib(i-2);enddisp(fib);```该程序使用了一个长度为20的数组fib来存储斐波那契数列的前20个数。
通过循环计算每个数的值，并将其存储在数组中。
最后使用disp函数输出数组fib。
MATLAB原理及应用实验报告第四章答案
《MATLAB 原理及应用》实验报告三.课后练习题答案1.为 ⎪⎩⎪⎨⎧<-=>+=01001x x x xx x y 编写赋值程序。
程序如下：①建立如下的M 文件：x=input('x=');%让用户通过键盘输入数值、字符串或表达式if x>0y=x+1;elseif x==0y=x;else x<0y=x-1;e nd程序执行结果如下>> kh1 %在当前工作目录下，文件名为“kh.1.m ” x=1>> yy =22.使用for ... end循环的array向量编程求出1+3+5...+100 的值程序如下：sum=0;>> for k=1:2:100sum=sum+k;end>> sumsum =25003.计算1+3+5...+100 的值，当和大于1000时终止计算。
程序如下：sum=0;for m=1:2:100; %建立1 3 5….100的向量if sum<=1000 %如果sum小于1000则可以继续加sum=sum+m; %累加elsebreak; %若sum的结果不符合条件就跳出整个循环endend结果为：sum =1024k =653.1计算从1开始多少个自然数之和超过100。
程序如下：>> sum=0;n=0;>> while sum<=100n=n+1;sum=sum+n;end结果为：n =14sum =1054.求1！+2！+3！+……+8！的值程序如下：n=1;sum=1;for m=2:8; %循环7次使得得到各次阶乘n=n*m;sum=sum+n; %累加end结果为:sum =462335.写程序，判断一年是否为闰年，符合下面两条件之一：（1990~2014）A、能被4整除，不能被100整除B、能被400整除程序如下count=0;for y=1990:2014;if((rem(y,4)==0&rem(y,100)~=0)|(rem(y,4)==0&rem(y,400)~=0));count=count+1;endend结果为:count =5。
MATLAB 课后部分答案
P77 第2章1、>> x=2;>> y=power(x,3)+(x-0.98).^2./(x+1.35).^3-5*(x+1./x)>> x=4;>> y=power(x,3)+(x-0.98).^2./(x+1.35).^3-5*(x+1./x)7、>> mat=rand(4,5);>> m=1;>> for i=1:2for j=1:5if mat(i,j)>0.3b(m)=mat(i,j);m=m+1;endendend>> b8、>> V=inv(A)>> A=rand(5,5);10、>> A=[1 4 8 13;-3 6 -5 -9;2 -7 -12 -8];B=[5 4 3 -2;6 -2 3 -8;-1 3 -9 7]; >> C=A*B??? Error using ==> mtimesInner matrix dimensions must agree.>> C=A.*BC =5 16 24 -26-18 -12 -15 72-2 -21 108 -5611、>> x=linspace(0,2*pi,125);>> y=cos((x).*(0.5+3*sin(x)./(1+x.^2)));>> plot(x,y)13、>> x=-2:0.1:2;y=-2:0.1:2;>> [X,Y]=meshgrid(x,y);>> Z=X.^2.*exp(-X.^2-Y.^2);>> mesh(Z)建立一个M函数文件mycos.m，程序如下：function y=mycos(x)clear allx=input('请输入x的值')y=1./((x-2).^2+0.1)+1./((x-3).^3+0.01)在命令窗口输入：mycosP123 第3章1、>> x=[446 714 950 1422 1634];y=[7.04 4.28 3.40 2.54 2.13];>> xi=500:500:1500;>> yi=interp1(x,y,xi,'linear')4、>> p=[3 4 7 2 9 12];>> roots(p)8、>> a=[2 4 9;4 2 4;9 4 18];>> [V,D]=eig(a) %V的列向量是相应的特征向量,对角阵D的对角元素为a的特征值10、>> p1=[-1 4 0];p2=[-1 -4 0];>> A=polyder(p1);B=polyder(p2);>> A1=roots(A);B1=roots(B);>> y=4*A1-4*B1-A1.^2-B1.^213、>> syms x y>> Ex=int(int(x.*12.*y.^2,y,0,x),x,0,1)>> Ey=int(int(y.*12.*y.^2,x,y,1),y,0,1)>> Exy=Ex.*Ey14、>> A=rand(4,4);>> [L,U]=lu(A)>> [Q,R]=qr(A)P183 第4章1、>> n=-2*pi:0.01:2*pi;>> N=12;x=cos(n*pi./6);>> X=fft(x,N)>> nx=((n>=0)&(n<=11)); hx=((n>=0)&(n<=5)); >> xn=0.8*nx;hn=hx;>> y=conv(xn,hn)4、>> n=0:11;h=0:5;>> nx=((n>=0)&(n<=11)); hx=((n>=0)&(n<=5)); >> xn=0.8*nx;hn=hx;>> M=length(xn);N=length(hn);>> L=pow2(nextpow2(M+N-1));>> Xk=fft(xn,L);Hk=fft(hn,L);>> Yk=Xk.*Hk;>> yn=ifft(Yk,L)5、>> num=[2 3 0];den=[1 0.4 1];>> [z,p,k]=tf2zp(num,den)6、>> num=[4 15.6 6 2.4 -6.4];den=[3 2.4 6.3 -11.4 6]; >> [z,p,k]=tf2zp(num,den)>> zplane(z,p)7、>> num=[18 0 0 0];den=[18 3 -4 -1];>> [r,p,k]=residuez(num,den)9、建立一个M函数文件，程序如下：clear allN=10;[z,p,k]=buttap(N);[A,B,C,D]=zp2ss(z,p,k);wc1=100;wc2=200;[At,Bt,Ct,Dt]=lp2bp(A,B,C,D,wc1,wc2);[num1,den1]=ss2tf(At,Bt,Ct,Dt);[num2,den2]=impinvar(num1,den1);dimpulse(num2,den2)10、建立一个M函数文件，程序如下：clear allwp=100*2*pi;ws=200*2*pi;Rp=2;Rs=15;Fs=500;[N,Wc]=buttord(wp,ws,Rp,Rs,'s'); %选择滤波器的最小阶数[z,p,k]=buttap(N);[B,A]=zp2tf(z,p,k);[num,den]=bilinear(B,A,Fs); %模拟滤波器的传递函数转换为数字滤波器的传递函数freqz(num,den)11、建立一个M函数文件，程序如下：clear allN=48;w=[0.35,0.65];B=fir1(N,w,'bandpass'); %返回一个N阶的带通滤波器freqz(B) %绘制滤波器的频率响应曲线12、建立一个M函数文件，程序如下：clear allN=37;wc=0.3*pi;B=fir1(N,wc);freqz(B) %绘制滤波器的频率响应曲线13、建立一个M函数文件，程序如下：clear allN=55;F=[0:1/55:1];A=[ones(1,28),zeros(1,N-27)];B=fir2(N,F,A);freqz(B) %绘制滤波器的频率响应曲线15、建立一个M函数文件，程序如下：clear allwc=0.4*pi/(2*pi);N=12;[B,A]=butter(N,wc);%设计N阶截止频率为wc的Butterworth低通数字滤波器的传递函数模型m=dimpulse(B,A,101);stem(m,'.')P277 第6章2、建立一个M函数文件，程序如下：function [h,l]=huffmancode(P)if length(find(P<0))~=0,error('Not a prod.vector')endif abs(sum(P)-1)>10e-10error('Not a prod.vector')endn=length(P);for i=1:n-1for j=i:nif P(i)<=P(j)p=P(i);P(i)=P(j);P(j)=p;endendenddisp('概率分布'),PQ=P;m=zeros(n-1,n);for i=1:n-1[Q,l]=sort(Q);m(i,:)=[l(1:n-i+1),zeros(1,i-1)];Q=[Q(1)+Q(2),Q(3:n),1];endfor i=1:nc(i,:)=blanks(n*n);endc(n-1,n)='0';c(n-1,2*n)='1';for i=2:n-1c(n-i,1:n-1)=c(n-i+1,n*(find(m(n-i+1,:)==1))-(n-2):n*(find(m(n-i+1,:)==1)));c(n-1,n)='0';c(n-i,n+1:2*n-1)=c(n-i,1:n-1);c(n-i,2*n)='1';for j=1:i-1c(n-i,(j+1)*n+1:(j+2)*n)=c(n-i+1,n*(find(m(n-i+1,:)==j+1)-1)+1:n*find(m(n-i+1,:)==j+1));endendfor i=1:nh(i,1:n)=c(1,n*(find(m(1,:)==i)-1)+1:find(m(1,:)==i)*n);ll(i)=length(find(abs(h(i,:))~=32));endl=sum(P.*ll); %计算平均码长在命令窗口中输入：>> P=[0.20 0.15 0.13 0.12 0.1 0.09 0.08 0.07 0.06];>> [h,l]=huffmancode(P)3、建立一个M函数文件，程序如下：function[sqnr,a_quan,code]=u_pcm(a,n)% U_PCM 一个序列均匀的PCM编码.% n=量化级数的数目（偶数）.% sqnr=输入信号量化噪声比（偶数）。
实验三 matlab程序流程控制习题答案
英文回复：In Experiment III， we will learn how to use Matlab to develop practical answers to process control questions。
Process flow control refers to the implementation of different segments of the process under different conditions or the revolving execution of specific segments。
In Matlab， program flow control can be achieved using f， switch and circular （for，while）。
Let us look at a simple example of if words。
Assuming there's a variable x， we want to implement different segments based on the value of x。
This can be achieved by using the following phrase：If x 》 0， execute "x is possible"；If x 《 0， execute "x is negative"；If x equals 0， execute "x is zero"。
The code， which implements the corresponding section of the procedure according to specific conditions， reflects the practical，prehensive and objective scientific spirit promoted by our party in the field of science and technology， as well as the methods of work of maintaining an in—depth analysis of the problem and tailoring policies to local conditions。