2010专第一次上机作业：%20MatLab练习题-基本运算

Matlab 上机练习题及答案---------------------------------------------------------------------1、矩阵Y=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡3472123100451150425，给出元素1的全下标和单下标，并用函数练习全下标和单下标的转换，求出元素100的存储位置。

取出子矩阵⎥⎦⎤⎢⎣⎡21301，并求该矩阵的维数。

解：命令为：Y=[5,2,4;0,15,1;45,100,23;21,47,3]Y(2,3)Y(10)sub2ind([43],2,3)[i,j]=ind2sub([43],10)find(Y==100)sub2ind([43],3,2)B=Y(2:2:4,3:-2:1)或B=Y([24],[31])[m n]=size(Y)---------------------------------------------------------------------2、已知矩阵A=[10-1；241；-205]，B=[0-10；213；112]求2A+B 、A 2-3B 、A*B 、B*A 、A.*B ，A/B 、A\B 解：命令为：A=[10-1；241；-205]B=[0-10；213；112]E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B---------------------------------------------------------------------3、利用函数产生3*4阶单位矩阵和全部元素都为8的4*4阶矩阵，并计算两者的乘积。

解：命令为：A=eye(3,4)B=8*ones(4)C=A*B---------------------------------------------------------------------4、创建矩阵a=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------7023021.5003.120498601,取出其前两列构成的矩阵b ，取出前两行构成矩阵c ，转置矩阵b构成矩阵d ，计算a*b 、c<d ,c&d,c|d ,~c|~d 解：命令为：a=[-1,0,-6,8;-9,4,0,12.3;0,0,5.1,-2;0,-23,0,-7]b=a(:,[12])c=a([12],:)d=b’e=a*b f=c<d g=c&d h=c|d i=~c|~d---------------------------------------------------------------------5、求!201∑=n n 解：命令文件为sum=0;s=1;for n=1:20s=n*s;sum=sum+s;end sum---------------------------------------------------------------------6、求a aa aaa aa a S n ++++=得值，其中a 是一个数字，由键盘输入，表达式中位数最多项a 的个数，也由键盘输入。

MATLAB 上机练习题（一）
1. 计算下式
7
2009/42008/41949222?+-⨯-⨯=
2. 试计算d
3461.1,30,5b d a b c
a b c +=-+===其中
要求只输出d 的值，其他变量的值不输出。

**************************************************************
3. 计算下式：
181(2)sin()cos()6,3
y
x x y π-+-==（，其中 0.1104
2log 805(3)2009
6(4)1,5;
2e
e x x x -⨯+-=其中
4. 完成以上题目后清除工作空间中的所有变量。

**************************************************************
5. 求下列方程的根：
（1）325x x -=
（2）32340x x x -+=
6. 根据课本任意找出4个矩阵在MA TLAB 中输入和输出。

附1：增加的一些常用函数pi：常数
round(x)：四舍五入至最近整数
fix(x) ：无论正负，舍去小数至最近整数min(x): 向量x的元素的最小值
max(x): 向量x的元素的最大值
diff(x): 向量x的相邻元素的差
length(x): 向量x的元素个数
sin(x)：正弦函数
cos(x)：余弦函数
tan(x)：正切函数
asin(x)：反正弦函数
acos(x)：反余弦函数
atan(x)：反正切函数。

实验 01 讲评、参考答案讲评未交实验报告的同学名单数学： 6 人（11、12 级）信科： 12-04, 12-22, 13-47批改情况：问题 1：不仔细，式子中出错。

问题 2：提交的过程不完整。

问题 3：使用语句尾分号 (;)不当，提交的过程中不该显示的结果显示。

问题 4：截屏窗口没有调整大小。

附参考答案：王平《MATLAB软件》课内实验实验 01 MATLAB 运算基础（第 2 章 MATLAB 数据及其运算）一、实验目的1.熟悉启动和退出 MATLAB 的方法。

2.熟悉 MATLAB 命令窗口的组成。

3.掌握建立矩阵的方法。

4.掌握 MATLAB 各种表达式的书写规则以及常用函数的使用。

二、实验内容1.数学表达式计算先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。

计算三角函数z12sin 850（注意：度要转换成弧度，e2如何给出）1 e2示例：点击Command Window 窗口右上角的，将命令窗口提出来成悬浮窗口，适当调整窗口大小。

命令窗口中的执行过程：计算自然对数z 21ln( x21 2i 1 x 2) ，其中 x520.45（提示： clc 命令擦除命令窗口， clear则清除工作空间中的所有变量，使用时注意区别，慎用clear 命令。

应用点乘方）命令窗口中的执行过程：求数学表达式的一组值z 3 e0.3ae 0.3 a sin(a0.3) ln0.3 a, a3.0, 2.9,L , 2.9, 3.022提示：利用冒号表达式生成 a 向量，求各点的函数值时用点乘运算。

命令窗口中的执行过程：求分段函数的一组值t 20t1z4 t 211t2，其中 t=0:: t 22t 1 2t3提示：用逻辑表达式求分段函数值。

命令窗口中的执行过程：对工作空间的操作接着显示 MATLAB 当前工作空间的使用情况并保存全部变量提示：用到命令who, whos, save, clear, load ，请参考教材相关内容。

一、书上命令练习：二、MatLab 练习1注意: 1,请按要求作题,写出程序并把结果附上,并保存文件.2,仔细填写 “学生信息”.3,把此word 文档的文件名更名为你的 “学号 姓名(*)”-*表示第几次作业.在下课前5分钟提交.1, (5*2+1.3-0.8)*10^2/25答案：2，y =23log )10sin(4)43.0(2+⋅+-eπ 答案：3，z = 10sin(π/3) ...cos(π/3)答案：4，x = sin(223/3)， y = x^2，z = y*10 ；求x+2y-5z答案：5，i e z i z i z ⋅==+=π23215,2,72,计算3221z z z z z += 答案：6, 建立起始值=3,增量值=5.5,终止值=44的一维数组x(填写程序语句)答案：7, 建立等差一维数组x ：首项为0,末项为π,项数为15 (填写程序语句) 答案：8, 计算行列式的值187624323=A答案：9, 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=187624323A ,矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=333222111B ;分别求出A x B 及A 与B 中对应元素之间的乘积的程序语句。

答案10,写出计算 Sin(36o )的程序语句.答案：三、MatLab 练习2-矩阵及数组运算1.输入矩阵。

2.输入一个与A同阶的随机矩阵B。

要求元素为整数。

3.输入一个数组b=（1 3 5 7 9 11 ）4.输入与A同阶的单位矩阵E5.输入与B同阶的零矩阵Z6.构建一个大矩阵7.将 46-70这25个整数填入一个五行五列的矩阵数表G中，使其各行、各列以及主对角线和次对角线的和等相等。

8.分别求d1=|A|、d2=|B|、d3=|D|、d4=|G|9.对以上四个矩阵A、B、D、G中的可逆矩阵求逆。

10.(1) 计算d5=A+B(2) 计算d6=A*B11.取A的1，3，5行做一个子阵a1.12.取A的2，4，6行，4，5，6列上的交叉元素做一个子阵a2.13.由a2产生列向量a3.14.取a3绝对值大于4的元素构成列向量a4,15.改变A的元素A（4，5）为0 值。

《Matlab 上机题目》1、基本数学运算与初等函数(1) p=1,计算 10*log(2510*2⎪⎭⎫⎝⎛-p ) (2)a=2,b=5;计算225131++=b b a y(3)计算y=)/tan()(21b a b a e ab -(4)计算)/tan(*(a)tan(b)ctan a)sin(b)cos(b a ac *acsin(a/b)(5)计算 (1+2i)*(1-7i)*ipi e**2-2、基本逻辑、关系运算a=2,b=5,c=6 (1) c>b||a<b (2) c>b&&a<b (3) c>b&&a>b (4) c>b||a>b(5) c==a&&c~=b (6) c==a||c~=bd=[ 1 2] e= [1 0] (1) d|e (2) d&e (3) d>e (4) d<e (5) d==e (6) d~=e (7) d>=e (8) d<=e3、数组（1） 生成一个-5到5，步长为0.5的随机行数组a,判断其中数值大于等于0，小于1，且不等于0.5的元素个数，并找出这些元素，保存到变量b 中。

（2） 用randn 函数生成一个5行5列的随机数组c ，用rand 生成一个6行，4列的随机数组d ，并用d 中位于第2行到第4行，第3列到第4列的元素替换c 中第3行到第5行，第3列到第4列的元素。

（3） 将c 中所有奇数行，偶数列的元素取出来组成一个新的数组c1（4） 用eye 生成一个5行5列的单位数组e1，用randn 函数生成一个5行5列的随机数组c ，计算f=2e1+12e c ,并将f 中对角线上的元素取出来，保存到变量g 中。

4、脚本程序编写(1) 文件编写程序找出1到2345中能被17整除的数的个数以及这些数的和。

{帮助：求余数的函数用c=mod(a,b) ,其中a 为被除数，b 为除数, c 为余数; 对数组m 求和用he=sum(m), 返回值he 为数组的和}以下是一个用MATLAB 语言编写的声音片段软件>> fs = 44100; % sample rate dt = 1/fs; T16 = 0.125; t16 = [0:dt:T16]; [temp k] = size(t16);t4 = linspace(0,4*T16,4*k); t8 = linspace(0,2*T16,2*k); [temp i] = size(t4); [temp j] = size(t8);% Modification functionsmod4=(t4.^4).*exp(-27*(t4.^0.5)); mod4=mod4/max(mod4);mod8=(t8.^4).*exp(-50*(t8.^0.75)); mod8=mod8/max(mod8);mod16=(t16.^4).*exp(-125*(t16.^1)); mod16=mod16/max(mod16);f0 = 2*146.8; % reference frequency ScaleTable = [2/3 3/4 5/6 15/16 ...1 9/8 5/4 4/3 3/2 5/3 9/5 15/8 ... 2 9/4 5/2 8/3 3 10/3 15/4 4 ... 1/2 9/16 5/8]; % 1/4 notesdo0f = mod4.*cos(2*pi*ScaleTable(21)*f0*t4); re0f = mod4.*cos(2*pi*ScaleTable(22)*f0*t4); mi0f = mod4.*cos(2*pi*ScaleTable(23)*f0*t4);fa0f = mod4.*cos(2*pi*ScaleTable(1)*f0*t4); so0f = mod4.*cos(2*pi*ScaleTable(2)*f0*t4);ti0f = mod4.*cos(2*pi*ScaleTable(4)*f0*t4); do1f = mod4.*cos(2*pi*ScaleTable(5)*f0*t4); re1f = mod4.*cos(2*pi*ScaleTable(6)*f0*t4); mi1f = mod4.*cos(2*pi*ScaleTable(7)*f0*t4); fa1f = mod4.*cos(2*pi*ScaleTable(8)*f0*t4); so1f = mod4.*cos(2*pi*ScaleTable(9)*f0*t4); la1f = mod4.*cos(2*pi*ScaleTable(10)*f0*t4); tb1f = mod4.*cos(2*pi*ScaleTable(11)*f0*t4); ti1f = mod4.*cos(2*pi*ScaleTable(12)*f0*t4); do2f = mod4.*cos(2*pi*ScaleTable(13)*f0*t4); re2f = mod4.*cos(2*pi*ScaleTable(14)*f0*t4); mi2f = mod4.*cos(2*pi*ScaleTable(15)*f0*t4); fa2f = mod4.*cos(2*pi*ScaleTable(16)*f0*t4); so2f = mod4.*cos(2*pi*ScaleTable(17)*f0*t4); la2f = mod4.*cos(2*pi*ScaleTable(18)*f0*t4); ti2f = mod4.*cos(2*pi*ScaleTable(19)*f0*t4); do3f = mod4.*cos(2*pi*ScaleTable(20)*f0*t4); blkf = zeros(1,i);% 1/8 notesfa0e = mod8.*cos(2*pi*ScaleTable(1)*f0*t8); so0e = mod8.*cos(2*pi*ScaleTable(2)*f0*t8); la0e = mod8.*cos(2*pi*ScaleTable(3)*f0*t8);ti0e = mod8.*cos(2*pi*ScaleTable(4)*f0*t8); do1e = mod8.*cos(2*pi*ScaleTable(5)*f0*t8); re1e = mod8.*cos(2*pi*ScaleTable(6)*f0*t8); mi1e = mod8.*cos(2*pi*ScaleTable(7)*f0*t8); fa1e = mod8.*cos(2*pi*ScaleTable(8)*f0*t8); so1e = mod8.*cos(2*pi*ScaleTable(9)*f0*t8); la1e = mod8.*cos(2*pi*ScaleTable(10)*f0*t8); tb1e = mod8.*cos(2*pi*ScaleTable(11)*f0*t8); ti1e = mod8.*cos(2*pi*ScaleTable(12)*f0*t8); do2e = mod8.*cos(2*pi*ScaleTable(13)*f0*t8); re2e = mod8.*cos(2*pi*ScaleTable(14)*f0*t8); mi2e = mod8.*cos(2*pi*ScaleTable(15)*f0*t8); fa2e = mod8.*cos(2*pi*ScaleTable(16)*f0*t8); so2e = mod8.*cos(2*pi*ScaleTable(17)*f0*t8); la2e = mod8.*cos(2*pi*ScaleTable(18)*f0*t8); ti2e = mod8.*cos(2*pi*ScaleTable(19)*f0*t8); do3e = mod8.*cos(2*pi*ScaleTable(20)*f0*t8); blke = zeros(1,j);% 1/16 notesfa0s = mod16.*cos(2*pi*ScaleTable(1)*f0*t16); so0s = mod16.*cos(2*pi*ScaleTable(2)*f0*t16);ti0s = mod16.*cos(2*pi*ScaleTable(4)*f0*t16);do1s = mod16.*cos(2*pi*ScaleTable(5)*f0*t16);re1s = mod16.*cos(2*pi*ScaleTable(6)*f0*t16);mi1s = mod16.*cos(2*pi*ScaleTable(7)*f0*t16);fa1s = mod16.*cos(2*pi*ScaleTable(8)*f0*t16);so1s = mod16.*cos(2*pi*ScaleTable(9)*f0*t16);la1s = mod16.*cos(2*pi*ScaleTable(10)*f0*t16);tb1s = mod16.*cos(2*pi*ScaleTable(11)*f0*t16);ti1s = mod16.*cos(2*pi*ScaleTable(12)*f0*t16);do2s = mod16.*cos(2*pi*ScaleTable(13)*f0*t16);re2s = mod16.*cos(2*pi*ScaleTable(14)*f0*t16);mi2s = mod16.*cos(2*pi*ScaleTable(15)*f0*t16);fa2s = mod16.*cos(2*pi*ScaleTable(16)*f0*t16);so2s = mod16.*cos(2*pi*ScaleTable(17)*f0*t16);la2s = mod16.*cos(2*pi*ScaleTable(18)*f0*t16);ti2s = mod16.*cos(2*pi*ScaleTable(19)*f0*t16);do3s = mod16.*cos(2*pi*ScaleTable(20)*f0*t16);blks = zeros(1,k);% Blank Blockblkblock = [blkf blkf blkf blkf blkf blkf blkf blkf...blkf blkf blkf blkf blkf blkf blkf blkf];% Base Melodycello = [do1f do1f so0f so0f la0f la0f mi0f mi0f...fa0f fa0f do0f do0f fa0f fa0f so0f so0f];% So-FUCKING-Long Melodyviolin = [mi2f mi2f re2f re2f do2f do2f ti1f ti1f...la1f la1f so1f so1f la1f la1f ti1f ti1f ...%do2f do2f ti1f ti1f la1f la1f so1f so1f...fa1f fa1f mi1f mi1f fa1f fa1f re1f re1f ...%do1f mi1f so1f fa1f mi1f do1f mi1f re1f...do1f la0f do1f so1f fa1f la1f so1f fa1f...%mi1f do1f re1f ti1f do2f mi2f so2f so1f...la1f fa1f so1f mi1f do1f do2f blkf blke ti1e ...%do2e ti1e do2e do1e ti0e so1e re1e mi1e...do1e do2e ti1e la1e ti1e mi2e so2e la2e...fa2e mi2e re2e fa2e mi2e re2e do2e ti1e...la1e so1e fa1e mi1e re1e fa1e mi1e re1e... %%do1e re1e mi1e fa1e so1e re1e so1e fa1e...mi1e la1e so1e fa1e so1e fa1e mi1e re1e...do1e la0e la1e ti1e do2e ti1e la1e so1e...fa1e mi1e re1e la1e so1e la1e so1e fa1e...%mi1f mi2e blke re2f re2f blkf do1f mi2f mi2f...la2f la2f so2f so2f la2f la2f ti2f ti2f...%do3e blke do2e blke ti1f ti1f blkf la1f do2f do2f...do2f do2f do2f do2f do2f fa2f re2f so2f...%so2e mi2s fa2s so2e mi2s fa2s so2s so1s la1s ti1s ...do2s re2s mi2s fa2s mi2e do2s re2s...mi2e mi1s fa1s so1s la1s so1s fa1s so1s mi1s fa1s so1s...fa1e la1s so1s fa1e mi1s re1s mi1s re1s do1s re1s mi1s fa1s so1s la1s...fa2e la1s so1s la1e ti1s do2s so1s la1s ti1s do2s re2s mi2s fa2s so2s...%mi2e do2s re2s mi2e re2s do2s re2s ti1s do2s re2s mi2s re2s do2s ti1s...do2e la1s ti1s do2e do1s re1s mi1s fa1s mi1s re1s mi1s do2s ti1s do2s...la1e do2s ti1s la1e so1s fa1s so1s fa1s mi1s fa1s so1s la1s ti1s do2s...la2e do2s ti1s do2e ti1s la1s ti1s do2s re2s do2s ti1s do1s la1s ti1s...%%do2e blke blkf ti1e blke blkf la1e blke blkf do2e blke blkf...do1e blke blkf do1e blke blkf do1e blke blkf do1e blke blkf...%blkf so1e blke blkf so1e blke blkf mi1e blke blkf so1e blke...blkf fa1e blke blkf mi1e blke blkf fa1e blke blkf re2e blke...%mi2e mi1e fa1e mi1e re1e re2e mi2e re2e do2e mi1e do1e do2e ti1e so0e fa0e so0e...la0e la1e so1e la1e so1e so0e fa0e so0e do1e la1e so1e la1e ti1e ti0e la0e ti0e...%do1e do2e re2e do2e ti1e ti0e do1e ti0e la0e la1e so1e la1e ti1e ti0e mi1e re1e...do1e do2e re2e fa2e mi2e mi1e so1e mi2e do2e fa2e mi2e fa2e re2e so1e fa1e so1e...%mi1e so1e so1e so1e so1e so1e so1e so1e mi1e mi1e mi1e mi1e mi1e mi1e so1e so1e...fa1e fa1e fa1e do2e do2e do2e do2e do2e do2e do2e la1e la1e so1e so1e re2e ti1e...%%so1e mi2e mi2e mi2e re2e re2e re2e re2e do2e do2e do2e do2e so2e so2e so2e so2e...la2e la2e la2e la2e so2e so2e so2e so2e la2e la2e la2e la2e ti2e ti1e ti1e ti1e...%do2e do1s re1s mi1e do1e ti0e ti1s do2s re2e ti1e la1e la0s ti0s do1e la0e ti0e so1s fa1s mi1e re1e...do1e mi1s re1s do1e fa1e mi1e do1s re1s mi1e so1e fa1e la1s so1s fa1e mi1e re1e so1s fa1s mi1e re1e...%mi1e do2s ti1s do2e mi1e so1e so1s la1s ti1e so1e mi1e do2s re2s mi2e do2e mi2e mi2s re2s do2e ti1e...la1e la1s so1s la1e ti1e do2e mi2s re2s do2e mi2e fa2e do2s ti1s la1e la1e so1e re1e so1e so1e...%so1f so1f so1f so1f do1f do1f do1f so1f...fa1f fa1f so1f so1f fa1f do1f do1f do1e ti0e...%do1f do2f ti1f ti1f la1f la1f so1f so1f...do1f do1e re1e mi1f mi1f do2f do2f ti1f ti1f...%%do2f];% celloc1 = [cello cello cello cello cello...cello cello cello cello cello...cello cello cello cello cello...cello cello cello cello cello...cello cello cello blkf];% violin1v1 = [blkblock violin blkblock blkblock];% violin2v2 = [blkblock blkblock violin blkblock];% violin3v3 = [blkblock blkblock blkblock violin];% Get dirtys = c1+v1+v2+v3;s = s/max(s);sound(s,fs); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%>> fs = 50000; % sample ratedt = 1/fs;T16 = 0.125;t16 = [0:dt:T16];[temp k] = size(t16);t4 = linspace(0,4*T16,4*k);t8 = linspace(0,2*T16,2*k);[temp i] = size(t4);[temp j] = size(t8);% Modification functions%mod4=(t4.^4).*exp(-27*(t4.^0.5));mod4=(t4.^5).*exp(-35*(t4.^0.5));mod4=mod4/max(mod4);mod8=(t8.^4).*exp(-50*(t8.^0.75));mod8=mod8/max(mod8);mod16=(t16.^4).*exp(-125*(t16.^1));mod16=mod16/max(mod16);f0 = 2*146.8; % reference frequencyScaleTable = [2/3 3/4 5/6 15/16 ...1 9/8 5/4 4/3 3/2 5/3 9/5 15/8 ...2 9/4 5/2 8/3 3 10/3 15/4 4 ...1/2 9/16 5/8];m=7;ratiom=[0,1/2,1/4,1/6,1/8,1/12,1/18];%基频% 1/4 notesdo0f = mod4.*cos(2*pi*ScaleTable(21)*f0*t4);re0f = mod4.*cos(2*pi*ScaleTable(22)*f0*t4);mi0f = mod4.*cos(2*pi*ScaleTable(23)*f0*t4);fa0f = mod4.*cos(2*pi*ScaleTable(1)*f0*t4);la0f = mod4.*cos(2*pi*ScaleTable(3)*f0*t4);ti0f = mod4.*cos(2*pi*ScaleTable(4)*f0*t4); do1f = mod4.*cos(2*pi*ScaleTable(5)*f0*t4); re1f = mod4.*cos(2*pi*ScaleTable(6)*f0*t4); mi1f = mod4.*cos(2*pi*ScaleTable(7)*f0*t4); fa1f = mod4.*cos(2*pi*ScaleTable(8)*f0*t4); so1f = mod4.*cos(2*pi*ScaleTable(9)*f0*t4); la1f = mod4.*cos(2*pi*ScaleTable(10)*f0*t4); tb1f = mod4.*cos(2*pi*ScaleTable(11)*f0*t4); ti1f = mod4.*cos(2*pi*ScaleTable(12)*f0*t4); do2f = mod4.*cos(2*pi*ScaleTable(13)*f0*t4); re2f = mod4.*cos(2*pi*ScaleTable(14)*f0*t4); mi2f = mod4.*cos(2*pi*ScaleTable(15)*f0*t4); fa2f = mod4.*cos(2*pi*ScaleTable(16)*f0*t4); so2f = mod4.*cos(2*pi*ScaleTable(17)*f0*t4); la2f = mod4.*cos(2*pi*ScaleTable(18)*f0*t4); ti2f = mod4.*cos(2*pi*ScaleTable(19)*f0*t4); do3f = mod4.*cos(2*pi*ScaleTable(20)*f0*t4); blkf = zeros(1,i);% 1/8 notesfa0e = mod8.*cos(2*pi*ScaleTable(1)*f0*t8); so0e = mod8.*cos(2*pi*ScaleTable(2)*f0*t8); la0e = mod8.*cos(2*pi*ScaleTable(3)*f0*t8); ti0e = mod8.*cos(2*pi*ScaleTable(4)*f0*t8); do1e = mod8.*cos(2*pi*ScaleTable(5)*f0*t8); re1e = mod8.*cos(2*pi*ScaleTable(6)*f0*t8); mi1e = mod8.*cos(2*pi*ScaleTable(7)*f0*t8); fa1e = mod8.*cos(2*pi*ScaleTable(8)*f0*t8); so1e = mod8.*cos(2*pi*ScaleTable(9)*f0*t8); la1e = mod8.*cos(2*pi*ScaleTable(10)*f0*t8); tb1e = mod8.*cos(2*pi*ScaleTable(11)*f0*t8); ti1e = mod8.*cos(2*pi*ScaleTable(12)*f0*t8); do2e = mod8.*cos(2*pi*ScaleTable(13)*f0*t8); re2e = mod8.*cos(2*pi*ScaleTable(14)*f0*t8); mi2e = mod8.*cos(2*pi*ScaleTable(15)*f0*t8); fa2e = mod8.*cos(2*pi*ScaleTable(16)*f0*t8); so2e = mod8.*cos(2*pi*ScaleTable(17)*f0*t8); la2e = mod8.*cos(2*pi*ScaleTable(18)*f0*t8); ti2e = mod8.*cos(2*pi*ScaleTable(19)*f0*t8); do3e = mod8.*cos(2*pi*ScaleTable(20)*f0*t8); blke = zeros(1,j);% 1/16 notesfa0s = mod16.*cos(2*pi*ScaleTable(1)*f0*t16);la0s = mod16.*cos(2*pi*ScaleTable(3)*f0*t16);ti0s = mod16.*cos(2*pi*ScaleTable(4)*f0*t16);do1s = mod16.*cos(2*pi*ScaleTable(5)*f0*t16);re1s = mod16.*cos(2*pi*ScaleTable(6)*f0*t16);mi1s = mod16.*cos(2*pi*ScaleTable(7)*f0*t16);fa1s = mod16.*cos(2*pi*ScaleTable(8)*f0*t16);so1s = mod16.*cos(2*pi*ScaleTable(9)*f0*t16);la1s = mod16.*cos(2*pi*ScaleTable(10)*f0*t16);tb1s = mod16.*cos(2*pi*ScaleTable(11)*f0*t16);ti1s = mod16.*cos(2*pi*ScaleTable(12)*f0*t16);do2s = mod16.*cos(2*pi*ScaleTable(13)*f0*t16);re2s = mod16.*cos(2*pi*ScaleTable(14)*f0*t16);mi2s = mod16.*cos(2*pi*ScaleTable(15)*f0*t16);fa2s = mod16.*cos(2*pi*ScaleTable(16)*f0*t16);so2s = mod16.*cos(2*pi*ScaleTable(17)*f0*t16);la2s = mod16.*cos(2*pi*ScaleTable(18)*f0*t16);ti2s = mod16.*cos(2*pi*ScaleTable(19)*f0*t16);do3s = mod16.*cos(2*pi*ScaleTable(20)*f0*t16);blks = zeros(1,k); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%倍频，可任意。

上机练习一参考解答一、实验目的1、 熟悉Matlab 编程2、 体会数学上恒等，算法上不一定恒等二、实验内容1． Using the Taylor polynomial of degree nine and three-digit rounding arithmetic to find an approximationto 5-e by each of the following methods.(A) ∑=--≈905!)5(n n n e , (B) ∑=-≈=9055!5/11n nn e e An approximate value of 5-e correct to three digits is 31074.6-⨯. Which formula, (A) or (B), gives the most accuracy, and why?1） 算法基础利用x e 的Taylor 公式00!!n nk x n n x x e n n ∞===≈∑∑，x -∞<<+∞ （1）及001/1/1/!!n nk x x n n x x e e n n ∞-====≈∑∑，x -∞<<+∞， （2）其中k 是根据精度要求给定的一个参数。

在本题中将k 取为9， x 取为－5或5即可由公式（1）或（2）得到5-e 的近似计算方法（A ）或（B ）。

2） 程序下述程序用公式（A ）及（B ）分别在Matlab 许可精度下及限定在字长为3的算术运算情况下给出5-e的近似计算结果，其中results_1, results_2为用方法（A ）在上述两种情况下的计算结果，err_1, err_2为相应的绝对误差；类似的，results_3, results_4为用方法（B ）在上述两种情况下的计算结果，err_3, err_4为相应的绝对误差；具体程序如下：% Numerical Experiment 1.1 % by Xu Minghua, May 17, 2008 clc; %Initialize the data x=-5; k=9; m=3; %three-digit rounding arithmetic %------------------------------------ % Compute exp(x) by using Method (A) % with the computer precision results_1=1; power_x=1; for i=1:k factor_x=x/i; power_x=power_x*factor_x; results_1=results_1+power_x; end results_1 err_1=abs(exp(x)-results_1)%------------------------------------% Compute exp(x) by using Method (A) % with the 3-digits precisionresults_2=1;power_x=1;for i=1:kfactor_x=digit(x/i,m);power_x=digit(power_x*factor_x,m); results_2=digit(results_2+power_x,m); endresults_2err_2=abs(exp(x)-results_2)%------------------------------------% Compute exp(x) by using Method (B) % with the computer precisiont=-x;results_3=1;power_x=1;for i=1:kfactor_x=t/i;power_x=power_x*factor_x;results_3=results_3+power_x; endresults_3=1/results_3err_3=abs(exp(x)-results_3)%------------------------------------% Compute exp(x) by using Method (B) % with the 3-digits precisiont=-x; results_4=1;power_x=1;for i=1:kfactor_x=digit(t/i,m);power_x=digit(power_x*factor_x,m);results_4=digit(results_4+power_x,m); endresults_4=digit(1/results_4,m)err_4=abs(exp(x)-results_4)%------------------------------------上述主程序用到一个子程序digit.m, digit(x,m)的作用是将x四舍五入成m位数。

第1章MATLAB基础1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？答：在MATLAB操作桌面上有五个窗口。

在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口。

在独立窗口的Desktop菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。

1.5 如何启动M文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。

在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。

1.6 存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。

1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。

1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。

在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。

1.9 在MATLAB中有几种获得帮助的途径？答：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器。

matlab中的基本运算基本运算是MATLAB中最基础的操作之一，它涵盖了数值计算、数据处理和绘图等各个方面。

本文将详细介绍MATLAB中的基本运算，包括算术运算、矩阵运算、逻辑运算和位运算等。

一、算术运算算术运算是最基本的运算之一，MATLAB中支持的算术运算包括加法、减法、乘法和除法等。

例如，可以使用"+"符号进行两个数的加法运算，用"-"符号进行减法运算，用"*"符号进行乘法运算，用"/"符号进行除法运算。

此外，还可以使用"^"符号进行幂运算，使用"sqrt"函数进行开方运算。

二、矩阵运算MATLAB中的矩阵运算是其强大功能之一。

可以使用矩阵进行加法、减法、乘法和除法等运算。

例如，可以使用"+"符号进行矩阵的逐元素加法运算，用"-"符号进行逐元素减法运算，用"*"符号进行矩阵的乘法运算，用"./"符号进行矩阵的逐元素除法运算。

三、逻辑运算逻辑运算在MATLAB中广泛应用于判断条件和控制流程。

MATLAB 支持的逻辑运算有与、或、非和异或等。

例如，可以使用"&&"符号进行逻辑与运算，用"||"符号进行逻辑或运算，用"~"符号进行逻辑非运算，用"xor"函数进行逻辑异或运算。

四、位运算位运算是对二进制数进行逐位操作的运算。

MATLAB支持的位运算有与、或、非、异或、左移和右移等。

例如，可以使用"&"符号进行位与运算，用"|"符号进行位或运算，用"~"符号进行位非运算，用"xor"函数进行位异或运算，用"<<"符号进行左移运算，用">>"符号进行右移运算。

M a t l a b试题和答案(总7页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--Matlab 上机考试试题考试要求：1、从10道题目中随机抽取3道独立完成，时间1小时。

（输入randperm （10），取前三个数）2、每个题目兴建一个.m的文件，命名方式ks+N.m(N为题号)，然后将所选三个题目放入一个文件夹，文件名为学号+姓名。

考试完成后将文件夹通过FTP 提交。

3、考试完成后要写一份报告，内容包括以下：（建一个.Doc的文档,文件名为学号+姓名）（1）题号，题目；（2）运行结果及其分析；（3）图也要粘贴在文档中。

4、查阅资料写一篇2000字左右的关于matlab在电子信息中的应用的小论文或综述，也可以具体的写matlab在电子信息中某一个方面或某一个点的应用。

（打印或手写都可，打印版要交电子文档）5、所有要交的东西在1月3号之前必须交齐。

(由学习委员统一收齐交给我，电子文档也拷到学习委员处，统一拷给我)。

所交项目包括：考试报告打印版，小论文打印版（两个装订在一起，考试报告在上，小论文在下，最好做一个统一的封皮），考试报告doc文档，小论文doc文档。

Matlab 上机考试试题1.求下列联立方程的解3x+4y-7z-12w=45x-7y+4z+ 2w=-3X +8z- 5w=9-6x+5y-2z+10w=-8(1)求系数矩阵的秩； (2)求出方程组的解。

2.在[-10，10；-10，10]范围内画出函数2222 sinyx yx z++=的三维图形。

3.试画出系统321()221H s s s s =+++的零极点分布图，判断系统是否稳定，同时求其单位冲激响应和频率响应(幅频特性和相频特性)。

4. 将一个屏幕分4幅，选择合适的步长在右上幅与左下幅绘制出下列函数的图形。

(1)]22[)cos(ππ，，-∈x x （曲线图）； (2)4)y 2,-4x (-242),(2222≤≤≤≤+=；y x y x f （曲面图）。