下载文档原格式
15、今有多项式P1(x)=x4-2x+1,P2(x)=x2+4x-0.5，要求先求得P(x)=P1(x)+P2(x)，然后计算xi=0.2*i各点上的P(xi)（i=0,1,2,…,5）值。
p1=[1.0 0.0 0.0 -2.0 1.0];>> p2=[0.0 0.0 1.0 4.0 -0.5];>> p1x=poly2sym(p1);p2x=poly2sym(p2);>> p=p1x+p2xp =x^4+2*x+1/2+x^2>> x=0:5;>> x.^4+2*x+1/2+x.^2ans =0.5000 4.5000 24.5000 96.5000 280.5000 660.50001、试个MATLAB的工作空间中建立以下2个矩阵：A＝[1 2]1234B⎡⎤=⎢⎥⎣⎦，求出矩阵A和B的乘积，并将结果赋给变量C。
>> A=[1 2]A =1 2>> B=[1 23 4]B =1 23 4>> C=A*BC =7 102、利用MATLAB提供的帮助信息，了解inv命令的调用格式，并作简要说明。
help invINV Matrix inverse.INV(X) is the inverse of the square matrix X.A warning message is printed if X is badly scaled ornearly singular.See also SLASH, PINV, COND, CONDEST, LSQNONNEG, LSCOV. Overloaded methodshelp gf/inv.mhelp zpk/inv.mhelp tf/inv.mhelp ss/inv.mhelp lti/inv.mhelp frd/inv.mhelp sym/inv.mhelp idmodel/inv.m3、使用help命令查询函数plot的功能以及调用方法，然后利用plot命令绘制函数y=sin(x)的图形，其中0xπ≤≤。
P 第一次实验答案1． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clc clearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)y1=diff(y) subplot(2,1,1) plot(x,y)subplot(2,1,2) plot(x(1:150),y1)2绘制极坐标系下曲线（a,b,n 自定数据）clc clear a=10; b=pi/2; n=5;theta=0:pi/100:2*pi; rho=a*cos(b+n*theta); polar(theta,rho)3. 列出求下列空间曲面交线的程序clc clearx=[-5:0.5:5];[X,Y]=meshgrid(x); z1=X.^2-2*Y.^2;z2=X.*2-Y.*3; xlabel('x') ylabel('y') zlabel('z') surf(X,Y,z1) hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5); x1=X(k) y1=Y(k)z3=x1.^2-2*y1.^2 hold onplot3(x1,y1,z3,'*')4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y为纵坐标的曲线，要求有图形标注。
clc clearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2)); plot(x,y,'b*-'); title('绘图'); xlabel('x 坐标'); ylabel('y 坐标'); legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解 81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y xclc cleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];b=[4,-3,9,-8]; c=b/a; x=c(1,1) y=c(1,2) z=c(1,3) w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1，采用不同插值方法进行插值，并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。
第一次上机:基础准备及入门2.1 指令窗运行入门2.1.1 最简单的计算器使用法【例2.5.1-1】求23)]47(212[÷-⨯+的算术运算结果。
（1）用键盘在MATLAB 指令窗中输入以下内容 >> (12+2*(7-4))/3^2（2）在上述表达式输入完成后，按【Enter 】键，该就指令被执行。
（3）在指令执行后，MATLAB 指令窗中将显示以下结果。
ans = 2【例2.5.1-2】简单矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321A 的输入步骤。
（1）在键盘上输入下列内容A = [1,2,3; 4,5,6; 7,8,9]（2）按【Enter 】键，指令被执行。
（3）在指令执行后，MATLAB 指令窗中将显示以下结果：A =1 2 3 4 5 6 7 8 9【例2.5.1-3】矩阵的分行输入A=[1,2,3 4,5,6 7,8,9]（以下是显示结果）A =1 2 3 4 5 6 7 8 9【例2.5.1-4】指令的续行输入（以下格式在除Notebook 外的MATLAB 环境中可运行）S = 1 – 1/2 + 1/3 –1/4 + 1/5 – 1/6 ... + 1/7 – 1/8 S =0.63452.1.2数值、变量和表达式2.1.2.1复数和复数矩阵【例2.5.2.5-1】复数ieziziz63212,21,43π=+=+=表达，及计算321zzzz=。
（1）经典教科书的直角坐标表示法z1= 3 + 4iz1 =3.0000 +4.0000i（2）采用运算符构成的直角坐标表示法和极坐标表示法z2 = 1 + 2 * i %运算符构成的直角坐标表示法z3=2*exp(i*pi/6) %运算符构成的极坐标表示法z=z1*z2/z3z2 =1.0000 +2.0000iz3 =1.7321 + 1.0000iz =0.3349 + 5.5801i【例2.5.2.5-2】复数矩阵的生成及运算A=[1,3;2,4]-[5,8;6,9]*iB=[1+5i,2+6i;3+8*i,4+9*i]C=A*BA =1.0000 - 5.0000i 3.0000 - 8.0000i2.0000 - 6.0000i 4.0000 - 9.0000iB =1.0000 + 5.0000i2.0000 + 6.0000i3.0000 + 8.0000i4.0000 + 9.0000iC =1.0e+002 *0.9900 1.1600 - 0.0900i1.1600 + 0.0900i 1.3700【例2.5.2.5-3】求上例复数矩阵C的实部、虚部、模和相角。
P第一次实验答案要求以0.01秒为间隔，求出y的151个点, 并求出其导数的值和曲线。
clcclearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*si n(4*sqrt(3)*x+pi/3) y1=diff(y)subplot(2,1,1)Plot(x,y)subplot(2,1,2) plot(x(1:150),y1)2绘制极坐标系下曲线(a,b,n自定数据)— a cos b n vclccleara=10;b=pi/2;n=5;theta=0:pi/100:2*pi; rho=a*cos(b+n*theta);polar(theta,rho) z2=X.*2-Y.*3;xlabel( 'x')ylabel( 'y')zlabel( 'z') surf(X,Y,z1)hold on surf(X ,Y, z2)k=fi nd(abs(z1-z2)<0.5);x1=X(k)y1=Y(k) z3=x1.A2-2*y1.A2 hold onplot3(x1,y1,z3, '*')4、设y cos x 0.53sin x(1 x2)把x=0~2 n间分为101点，画出以x为横坐标，y 为纵坐标的曲线，要求有图形标注clcclearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+si n(x)*3./(1+x.A2));plot(x,y, 'b*-');title('绘图’)；xlabel( 'x 坐标');ylabel( 'y 坐标');legend('原函数')gtext( 'y=cos(x)(0.5+3*sin(x)/(1+xA2))' )3.列出求下列空间曲面交线的程序乙=x2 _2y2z2 = 2x _ 3yclcclearx=[-5:0.5:5];[X,Y]=meshgrid(x);z1=X.A2-2*Y.A2; 5、求下列联立方程的解3x 4y - 7z - 12w 二45x - 7y 4z 2w - -3x 8z - 5w = 9-6x 5y - 2z 10w = -8clccleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10]; b=[4,-3,9,-8];第二次试验答案1、编制m文件，等待键盘输入，输入密码20120520 ,密码正确，显示输入密码正确，程序结束；否则提示，重新输入。
Matlab 上机练习二班级 学号 姓名按要求完成题目，并写下指令和运行结果。
（不需要画图）1、 求⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x 的共轭转置。
>> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i];>> x’ans =- -- ++ -- ++ -2、计算⎥⎦⎤⎢⎣⎡=572396a 与⎥⎦⎤⎢⎣⎡=864142b 的数组乘积。
>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8];>> a.*bans =12 36 38 42 403、 对于B AX =，如果⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=753467294A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=282637B ，求解X 。
>> A=[4 9 2;7 6 4;3 5 7];>> B=[37 26 28]’;>> X=A\BX =4、 ⎥⎦⎤⎢⎣⎡-=463521a ，⎥⎦⎤⎢⎣⎡-=263478b ，观察a 与b 之间的六种关系运算的结果。
>> a=[1 2 3;4 5 6];>> b=[8 –7 4;3 6 2];>> a>bans =0 1 01 0 1 >> a>=bans =0 1 01 0 1 >> a<bans =1 0 1 0 1 0 >> a<=bans =1 0 1 0 1 0 >> a==bans =0 0 0 0 0 0>> a~=bans =1 1 11 1 15、[]7.0=-a，在进行逻辑运算时，a相当于什么样的逻辑量。
82.05-相当于a=[1 1 0 1 1]。
第4章图形处理功能1 内容简介基本内容主要包括：（1）二维图形（2）三维图形（3）图形处理的基本技术2 达到的目标（1）掌握二维图形的绘制。
（2）掌握三维图形的绘制。
（3）掌握图形处理的基本技术3 具体内容３．1 二维图形３．1．1 基本绘图命令（1）当plot函数仅有一个输入变量例4-1y=[5 2 3 8 5]; %y 行矩阵plot(y) %一条线例4-2y=[5 2 3 8 5;2 4 3 1 5;1 1 1 1 1]; %y 矩阵plot(y) %5条线,等于矩阵的列数（2）当plot函数有两个输人变量例4-3x=0:0.01*pi:pi;y=sin(x).*cos(x);plot(x,y)例4-3x=0:0.01*pi:pi;y=[sin(x);cos(x); sin(x).*cos(x)];plot(x,y)例4-4x1=0:0.01*pi:pi;x2=pi:0.01*pi:2*pi;x=[x1' x2'];y=[sin(x1') cos(x2')];plot(x,y)例4-5x1=1:5;x2=6:10;y1=x1;y2=2*x2;plot([x1;x2],[y1;y2])%plot([x1' x2'],[y1' y2'])（3）当plot函数有三个输入变量时MATLAB语言中提供的对曲线的线型、颜色以及标识的控制符如表4．l所示。
例4-6 绘制带有显示属性设置的二维图形。
x=0.5*pi: 0.1*pi:2*pi;y=sin(x);z=cos(x);plot (x, y, '--ko', x, z, '-. r*')３．1．2 特殊的二维图形函数（1）特殊坐标系的二维图形函数（ａ）对数坐标例4-7 绘制X坐标为对数坐标的二维图形。
x=0.5*pi: 0.1*pi:2*pi;y=sin(x);semilogx (x, y, '-ro')（ｂ）极坐标例4-8绘制极坐标下的二维图形。
实验一 Matlab使用方法和程序设计一、实验目的1、熟悉MATLAB上机环境，练习MATLAB命令、m文件；2、熟悉Matlab的数据表示、基本运算和程序控制语句；3、熟悉Matlab图形绘制、数据处理；4、熟悉Matlab程序设计的基本方法。
二MA TLAB简介MATLAB环境是一种为数值计算、数据分析和图形显示服务的交互式的环境。
MATLAB 有3种窗口，即：命令窗口（The Command Window）、m-文件编辑窗口（The Edit Window）和图形窗口（The Figure Window），而Simulink另外又有Simulink模型编辑窗口。
1．命令窗口（The Command Window）当MA TLAB启动后，出现的最大的窗口就是命令窗口。
用户可以在提示符“>>”后面输入交互的命令，这些命令就立即被执行。
在MA TLAB中，一连串命令可以放置在一个文件中，不必把它们直接在命令窗口内输入。
在命令窗口中输入该文件名，这一连串命令就被执行了。
因为这样的文件都是以“.m”为后缀，所以称为m-文件。
2．m-文件编辑窗口（The Edit Window）我们可以用m-文件编辑窗口来产生新的m-文件，或者编辑已经存在的m-文件。
在MATLAB主界面上选择菜单“File/New/M-file”就打开了一个新的m-文件编辑窗口；选择菜单“File/Open”就可以打开一个已经存在的m-文件，并且可以在这个窗口中编辑这个m-文件。
3．图形窗口（The Figure Window）图形窗口用来显示MA TLAB程序产生的图形。
图形可以是2维的、3维的数据图形，也可以是照片等。
二、实验内容1、帮助命令2、矩阵运算（1）矩阵的乘法和乘方已知A=[1 2;3 4]:B=[5 5;7 8]:求A^2*B( 2 )矩阵除法已知A=[1 2 3:4 5 6:7 8 9]:B=[1 0 0:0 2 0:0 0 3],求矩阵左除A\B,右除A/B。
13sin 25457.0145.175.375.0222-++=---y e x z x y x练习一：1) 以两种方式打开 MA TLAB 工作窗口，进入MATLAB 6.0 的工作环境，并尝试用不同的方式退出。
2) 尝试、熟悉 MA TLAB 6.0 的各栏菜单以及各个工具栏的功能。
3) 重新启动 MATLAB 6.0，进入 MATLAB 工作窗口，用 who 命令查看当前工作空间内有无变量及其值。
4) 绘制函数y=cos(5x+2)/sin(3x+1) 的图像，并求解当 x=2 时的函数值。
5) 此时再次用 who 命令查看工作空间内的变量名及其值，与（ 3）比较，同时用whos 命令查看变量，比较与 who 命令的不同。
6) 熟练掌握 MATLAB 的通用命令。
7) 练习并熟练掌握 MA TLAB 的帮助命令，学会利用 MATLAB 的帮助信息。
8) 用 lookfor 命令查找函数 cos 的信息，并与 help 命令查找的结果相比较，注意采用两种命令之间的差别。
练习二： 1) 创建一个表达式式并求当x =1, y =2时的z 值。
2) 计算i x y 108++=的值，其中451-+-=x 。
3) 创建一个包含“ 她是我们班最好的学生”字符串。
练习三：1) 将如下矩阵 A 进行转置和求逆。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----=412303214A 2) 在MATLAB 环境下，用下面三条指令创建矩阵 C ，看输出怎样的结果。
a=2.7358;b=33/79; C =[1,2*a+i*b,b*sqrt(a);sin(pi/4),a+5*b,3.5+i]练习三：3) 用 rand 函数建立一个 3×3 维随机数值矩阵，并将该数值矩阵转换为符号矩阵，比较它们的不同。
4) 求下述矩阵的逆、特征根和矩阵的行列式值。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333231232221131211a a a a a a a a a A ，⎥⎦⎤⎢⎣⎡=22211211a a a a A 练习四： 1) 用曲面图命令 surf 表现函数 22y x z +=的图像。
实验1 MATLAB7.0基本编程实验目的：① 掌握脚本和函数；② 掌握matlab 中的变量和M 文件的流控制语句；③ 掌握函数的设计和实现。
实验要求:给出程序和实验结果。
实验内容：一、一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。
输出全部水仙花数。
二、从键盘输入若干个数，当输入10时结束输入，求这些数的平均值和它们之和。
三、求[1000，2000]之间第一个能被17整除的整数。
四、若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。
求[1,10000]之间的全部完数。
五、Fibonacci 数列定义如下：1f =12f =112,(2)n n n f f f n --=+>求Fibonacci 数列的第20项。
六、设计一个猜数游戏。
首先由计算机产生[1,100]之间的随机整数，然后由用户猜测所产生的随机数。
根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示“High”，小于则显示“Low”，等于则显示“You w o n”，同时退出游戏。
用户最多可以猜7次。
实验结果：一、 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。
输出全部水仙花数。
程序如下：for m=100:999m1=fix(m/100); %求m 的百位数字m2=rem(fix(m/10),10); %求m 的十位数字m3=rem(m,10); %求m 的个位数字if m==m1*m1*m1+m2*m2*m2+m3*m3*m3disp(m)endend153370371407二、从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。
程序如下：sum=0;cnt=0;val=input('Enter a number (end in 0):');while (val~=0)sum=sum+val;cnt=cnt+1;val=input('Enter a number (end in 0):');endif (cnt > 0)summean=sum/cntendEnter a number (end in 0):45Enter a number (end in 0):65Enter a number (end in 0):342Enter a number (end in 0):45Enter a number (end in 0):67Enter a number (end in 0):12Enter a number (end in 0):0sum =576mean =96三、求[100，200]之间第一个能被21整除的整数。
1 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。
输出全部水仙花数。
for m=100:999m1=fix(m/100); %求m的百位数字m2=rem(fix(m/10),10); %求m的十位数字m3=rem(m,10); %求m的个位数字if m==m1*m1*m1+m2*m2*m2+m3*m3*m3disp(m)endend2.从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和。
sum=0;n=0;val=input('Enter a number (end in 0):');while (val~=0)sum=sum+val;n=n+1;val=input('Enter a number (end in 0):');endif (n > 0)summean=sum/nend3. 若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。
求[1,500]之间的全部完数。
for m=1:500s=0;for k=1:m/2if rem(m,k)==0s=s+k;endendif m==sdisp(m);endend4. 从键盘上输入数字星期，在屏幕上显示对应英文星期的单词。
function weekn=input('input the number:');if isempty(n)errror('please input !!')endif n>7|n<1error('n between 1 and 7')endswitch ncase 1disp('Monday')case 2disp('Tuesday')case 3disp('Wednesday')case 4disp('Thursday')case 5disp('Friday')case 6disp('Saturday')case 7disp('Sunday')end5. 某公司销售电脑打印机的价格方案如下：（）如果顾客只买一台打印机，则一台的基本价格为$150。
要求：实验前提前预习，实验时独立上机完成每个练习题的编程，执行出正确结果。
把上机执行后的程序理解 并完整的写在课堂作业本上，以备教务处检查与自己考试复习用。
2012年下学期MATLAB 语言上机实验内容 使用班级：1112201、202、203姓名： 班级： 学号：实验一: matlab 基础与符号计算1、指令>>A= [2 2 ; 2 2];>>A .^2= 。
>>A^2= 。
2、对3、数组>>x=10:-2:2; y=1:5;>>(x<5)&(y>2)，ans= 。
>> y==2|x~=2， ans= 。
3、（1）对于0>x ，求12011122+∞=∑⎪⎭⎫⎝⎛+-+k k x x k 。
（提示：理论结果为x ln ）（2）求∑∞=0k kx 的符号解，并进而用该符号解求∑∞=-0)31(k k，∑∞=0)1(k kπ，∑∞=03k k的准 确值。
（提示：注意subs 的使用）4、（1）通过符号计算求t t y sin )(=的导数dtdy。
（2）然后根据此结果，求-=0t dt dy 和2π=t dtdy。
5、指令>>A= [2 1 ; 1 3];>>A .^2-A= 。
>>A^2+A= 。
6、数组>>x=10:-2:2; y=5:-1:1;>>(x>=4)&(y>=2)，ans= 。
>> y==3|x~=3， ans= 。
7计算二重积分⎰⎰+211222)(x dydx y x 。
8在]2,0[π区间，画出dt ttx y x⎰=0sin )(曲线，并计算)5.4(y 。
9设系统的冲激响应为tet h 3)(-=，求该系统在输入t t u cos )(=，0≥t 作用下的输出。
（提示：运用卷积进行计算）10求⎪⎩⎪⎨⎧>≤⎪⎪⎭⎫⎝⎛-=τττt t t A t f 01)(的Fourier 变换，并画出2,2==τA 时的幅频谱。
11、方程求方程的解初始条件为，,4)0(r 3,r(0),423'2''===++-te r r r12、求Tk ke k f )(λ-=的Z 变换表达式。
13、求方程2,122==+xy y x 的解。
14求边值问题1)0(,0)0(,34,43==+-=+=g f g f dxdg g f dx df 的解。
15求4633)(23++++=s s s s s F 的Laplace 反变换。
实验2: 数组（矩阵）及数值计算1、指令>>A=ones(4,4); >>B=A(3,[2 3]); >>C=A(3,:)<1; 写出输出结果：>>B=_________________________。
>>C=_________________________。
2、 指令>>B=magic(3)；， 写出求矩阵B 的特征值D 和特征向量V 的语句：。
3、>>A=rand(100,100)；，用向量化运算的方法，用1条语句实现求矩阵A 中大于0.4的所有元素个数： 4要求在闭区间]2,0[π上产生具有10个等距采样点的一维数组。
试用两种不同的指令实现。
5由指令rand('state',0),A=rand(3,5)生成二维数组A ，试求该数组中所有大于0.5的元素的位置，分别求出它们的“全下标”和“单下标”。
6在时间区间 [0,10]中，绘制t e y t 2cos 15.0--=曲线。
要求分别采取“标量循环运算法”和“数组运算法”编写两段程序绘图。
7先运行指令x=-3*pi:pi/15:3*pi; y=x; [X,Y]=meshgrid(x,y); warning off; Z=sin(X).*sin(Y)./X./Y; 产生矩阵Z 。
（1）请问矩阵Z 中有多少个“非数”数据？（2）用指令surf(X,Y,Z); shading interp 观察所绘的图形。
（3）请写出绘制相应的“无裂缝”图形的全部指令。
8、指令>>A=3*ones(5,5); >>B=A(3,[1 2]); >>C=A(:，3)>1; 写出输出结果：>>B=_________________________。
>>C=_________________________。
9、 指令>>B=magic(4), 按要求分别写出指令：求矩阵B 的行列式：_______________________。
求矩阵B 的逆矩阵：_______________________。
求矩阵B 的迹： 。
求矩阵B 的特征值D 和特征向量V ： 。
10采用数值计算方法，画出dt ttx y x⎰=0sin )(在]10 ,0[区间曲线，并计算)5.4(y 。
11用quad 求取dx x exsin 7.15⎰--ππ的数值积分，并保证积分的绝对精度为910-。
12求函数5.08.12cos 5.1)5(sin )(206.02++-=t t t et t f t 在区间]5,5[-中的最小值点。
13设0)0(,1)0(,1)(2)(3)(22===+-dtdy y t y dt t dy dt t y d ，用数值法和符号法求5.0)(=t t y 。
14求矩阵b Ax =的解，A 为4阶魔方阵，b 是)14(⨯的全1列向量。
15求矩阵b Ax =的解，A 为4阶魔方阵，⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=4321b 。
16求解二元函数方程组⎩⎨⎧=+=-0)cos(0)sin(y x y x 的解。
（用符号法解）17已知有理分式)()()(x D x N x R =，其中)5.0)(3()(33++=x x x x N ，)125)(22()(232++-+=x x x x x D 。
（1）求该分式的商多项式)(x Q 和余多项式)(x r 。
（2）用程序验算)()()()(x N x r x Q x D =+是否成立。
18已知系统冲激响应为h(n)=[0.05,0.24,0.40,0.24,0.15,-0.1,0.1] ，系统输入u(n)由指令randn('state',1);u=2*(randn(1,100)>0.5)-1产生，该输入信号的起始作用时刻为0。
试用直杆图（提示：用stem 指令）画出分别显示该系统输入、输出信号的两张子图。
实验3: 图形绘制及m 文件编程1． 已知椭圆的长、短轴2,4==b a ，用“小红点线”画椭圆⎩⎨⎧==t b y ta x sin cos 。
2． 根据表达式θρcos 1-=绘制如图p5-2的心脏线。
（提示：采用极坐标绘线指令polar ） 3． A,B,C 三个城市上半年每个月的国民生产总值见表p5.1。
试画出如图p5-3所示的三城市上半年每月生产总值的累计直方图。
（提示：使用指令bar ）表p5.1 各城市生产总值数据（单位：亿元） 城市 1月 2月 3月 4月 5月 6月 A 170 120 180 200 190 220 B 120 100 110 180 170 180 C705080 100951204． 用绿实线绘制)sin(t x =，)cos(t y =，t z =的三维曲线，曲线如图p5-5所示。
（提示：使用plot3指令）-1-0.50.51-11510155． 在区域]3,3[,-∈y x ，绘制224y x xez --=的如图p5-6的三维（透视）网格曲面。
（不得使用ezmesh ）6． 在]4,4[,ππ-∈y x 区间里，根据表达式yx y x z ++=)sin(，绘制如图p5-7所示的曲面。
7． 试用图解法回答：（1）方程组⎪⎩⎪⎨⎧=+=++0))cos(sin(1.0)1(22y x y x y 有多少个实数解？（2）求出离0,0==y x 最近、且满足该方程组的一个近似解。
8.一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。
输出全部水仙花数。
9.Fibonacci 数组的元素满足Fibonacci 规则：12+++=k k k a a a ，),2,1( =k ；且121==a a 。
现要求该数组中第一个大于10000的元素。
10.请分别写出用for 和while 循环语句计算1000000210000002.02.02.012.0+++==∑= i iK 的程序。
此外，还请写出避免循环的数值、符号计算程序。
11.若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。
求[1,500]之间的全部完数。
12.随机输入20个数，按从小到大的顺序输出实验4: 连续、离散时间系统的simulink 仿真1利用SIMULINK 求解dxet I tx ⎰-= 02)(在区间]1 ,0[∈t 积分，并求出积分值)1(I 。
（参见例4.1-5）2、已知系统方程为4)0(r 3,r(0),423'''===++初始条件为，r r r ，请建立simulink 模型，并仿真执行。
（参见例4.1-9）3、设计一个数字低通滤波器)(z F ，从受噪声干扰的多频率混合信号)(t x 中获取10Hz 的信号。
)()1002cos(5.1)102sin()(t n t t t x +⋅⋅+⋅⋅=ππ （7.3） 在此，)2.0,0(~)(2N t n ，s sT k f k t ⋅=⋅=1。
采样频率取1000=s f （Hz ），即采样周期001.0=s T （秒）。
4、在图7.3-1所示的电路中，已知)(3.0H L =，)(3.0F C =，)(21Ω=R ，)(01.02Ω=R ，)(53Ω=R ，)(1)0(V V C -=-， )(1)0(A i L =-， )(10V V s =，开关K 在0=t 时闭合。
试采用SIMULINK 的SimPowerSystems 模块库器件进行元件级仿真，求 L i 和c V 。
图7.3-1 二阶RLC 电路。
