系统仿真导论-答案.docVIP

Chapter 0习题 MATLAB语言编程基础 PAGE 5/ NUMPAGES 6 PAGE PAGE 5 Chapter 0习题 MATLAB语言编程基础 目的： （1）熟悉MATLAB的使用。 （2）掌握MATLAB的编程的基本循环、判断语句、画图语句等。 （3）学会如何充分利用MATLAB的强大的帮助功能。 （4）掌握MATLAB函数的编写和调用。 用MATLAB语句输入矩阵和 ， 前面给出的是的矩阵，如果给出的命令将得出什么结果？ 【解答】 A=[1,2,3,4; 4,3,2,1; 2,3,4,1; 3,2,4,1] A = 1 2 3 4 4 3 2 1 2 3 4 1 3 2 4 1 A(5,6)=5 A = 1 2 3 4 0 0 4 3 2 1 0 0 2 3 4 1 0 0 3 2 4 1 0 0 0 0 0 0 0 5 B=[1+4j, 2+3j, 3+2j, 4+1j; 4+1j, 3+2j, 2+3j, 1+4j; 2+3j, 3+2j, 4+1j, 1+4j; 3+2j, 2+3j, 4+1j, 1+4j] B = 1.0000 + 4.0000i 2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i 4.0000 + 1.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 1.0000 + 4.0000i 2.0000 + 3.0000i 3.0000 + 2.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i 3.0000 + 2.0000i 2.0000 + 3.0000i 4.0000 + 1.0000i 1.0000 + 4.0000i 试用符号元素工具箱支持的方式表达多项式，并令，将替换成的函数，并写出最后的结果。（注：subs()函数可进行变量替换，用simple()函数进行结果化简）。 【解答】 syms x f=x^5+3*x^4+4*x^3+2*x^2+3*x+6 f = x^5 + 3*x^4 + 4*x^3 + 2*x^2 + 3*x + 6 syms s f=subs(f,x,(s-1)/(s+1)) f = (3*(s - 1))/(s + 1) + (2*(s - 1)^2)/(s + 1)^2 + (4*(s - 1)^3)/(s + 1)^3 + (3*(s - 1)^4)/(s + 1)^4 + (s - 1)^5/(s + 1)^5 + 6 simple(f) ans = 19 - (72*s^4 + 120*s^3 + 136*s^2 + 72*s + 16)/(s + 1)^5 采用循环的形式求出。 【解答】 循环形式： format long s=0; for i=0:63 s=s+2^i; end s s = 1.844674407370955e+019 符号运算形式： syms k symsum(2^k,0,63) ans = 18446744073709551615 用MATLAB语言实现下面的分段函数 。 【解答】 MATLAB M文件：getf.m function y = getf( x ) h=1;D=2; if xD y=h; elseif x-D y=-h; else y=h./(D.*x); end 测试getf.m的MATLAB命令及结果： getf(3) ans = 1 getf(1) ans = 0.500000000000000 getf(-1) ans = -0.500000000000000 getf(-3) ans = -1 编写一个矩阵相加函数mat_add( )，使得该函数能够接受任意多个矩阵，且其调用格式为。 【解答】 MATLAB M文件：mat_add.m function A=mat_add(varargin) A=0; for i=1:length(var