matlab习题与答案
Matlab及工程应用
----作业
胡章芳
电气工程与自动化（1）班
0513080109
燕如龙
重庆邮电大学电工程学院
习题一MATLAB运算基础
一、目的
1、熟悉启动和退出MATLAB的方法；
2、熟悉MATLAB命令窗口的组成；
3、掌握建立矩阵的方法；
4、掌握MATLAB各种表达式的书写规则以及常用函数的使用。
二、内容
1．先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。
提示：利用冒号表达式生成a向量，求各点的函数值时用点乘运算。>> z1=2*sin(85*pi/180)/(1+exp(2))
z1 =
0.2375
>> x=[2 1+2i;-0.45 5];
>> z2=1/2*log(x+sqrt(1+x^2))
z2 =
0.7114 - 0.0253i 0.8968 + 0.3658i
0.2139 + 0.9343i 1.1541 - 0.0044i
>> a=-3.0:0.1:3.0;
z3=1/2*(exp(0.3*a)-exp(0.3*a)).*sin(a+0.3)
z3 =
Columns 1 through 6
0 0 0 0 0 0 Columns 7 through 12
0 0 0 0 0 0 Columns 13 through 18
0 0 0 0 0 0 Columns 19 through 24
0 0 0 0 0 0 Columns 25 through 30
0 0 0 0 0 0 Columns 31 through 36
0 0 0 0 0 0 Columns 37 through 42
0 0 0 0 0 0 Columns 43 through 48
0 0 0 0 0 0 Columns 49 through 54
0 0 0 0 0 0 Columns 55 through 60
0 0 0 0 0 0 Column 61
>> A=[12 34 -4;34 7 87;3 65 7]; >> B=[1 3 -1;2 0 3;3 -2 7]; >> A+6*B
ans =
18 52 -10
46 7 105
21 53 49
>> A-B+eye(3)
ans =
12 31 -3
32 8 84
0 67 1
>> A*B
ans =
68 44 62
309 -72 596
154 -5 241
>> A.*B
ans =
68 0 261
9 -130 49
>> A^3
ans =
37226 233824 48604
247370 149188 600766
78688 454142 118820 >> A.^3
ans =
1728 39304 -64
39304 343 658503
27 274625 343 >> A/B
ans =
16.4000 -13.6000 7.6000
35.8000 -76.2000 50.2000
67.0000 -134.0000 68.0000
>> B\A
ans =
109.4000 -131.2000 322.8000
-53.0000 85.0000 -171.0000
-61.6000 89.8000 -186.2000
>> [A,B]
ans =
12 34 -4 1 3 -1
34 7 87 2 0 3
3 65 7 3 -2 7 >> [A([1,3],:);B^2]
ans =
3 65 7
4 5 1
11 0 19
20 -5 40
（1）求它们的乘积C；
（2）将矩阵C的右下角3x 2子矩阵赋给D；
（3）查看MATLAB工作空间的使用情况。
>> A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20;21 22 23 24 25 ]; >> B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 13 11];
>> C=A*B
C =
93 150 77
258 335 237
423 520 397
588 705 557
753 890 717
>> D=A(3:5,4:end)
D =
14 15
19 20
24 25
>> who
Your variables are:
A B C D
>> whos
Name Size Bytes Class
A 5x5 200 double array
B 5x3 120 double array
C 5x3 120 double array
D 3x2 48 double array
Grand total is 61 elements using 488 bytes
4．完成下列操作：
（1）求[100，999]之间能被21整除的数的个数。
提示；先利用冒号表达式，再利用find和Length函数。
>> a=100:999;
c=rem(a/21,1);
d=find(c==0);
n=length(d)
n =
43
习题二程序设计及函数和数据的可视化
一、目的
1、掌握建立和执行M文件的方法；
2、掌握利用if、switch、for、While语句实现的方法；
3、掌握定义和调用MATLAB函数的方法；
4、掌握绘制二维、三维图形的常用函数以及绘制图形的辅助操作。
二、内容
1．求下列分段函数的值
要求：（1）用if语句实现，分别输出x＝-5.0，-3.0，1.0，2.0，2.5，3.0，5.0时的y值。
>> x=input('please input x:');
if x<0&x~=-3
y=x^2+x-6
elseif 0<=x<10,x~=2&x~=3
y=x^2-5*x+6
else
y=x^2-x-1
end
please input x:-5.0
y =
14
please input x:-3.0
y =
30
please input x:1.0
y =
2
please input x:2.0
y =
please input x:2.5
y =
-0.2500
please input x:3.0 y = 0
please input x:5.0 y = 6
2．输入一个百分制成绩，要求输出成绩等级A,B,C,D,E.其中90~99分为A ，80~89 分为B ，70~79分为C ，60~69分为D ，60分以下为E 。 要求： 分别用if 语句和switch 语句实现。 function switch_case(number);
g=input('please input number:','s'); switch(g) case('90~99') disp('A'); case('80~89') disp('B'); case('70~79') disp('C'); case('60~69') disp('D'); case('<60') disp('E'); otherwise disp('输入错误!') end 3．根据
2
2
2
2
1111...
6
1
2
3
n
π
=
+
+
+，求π的近似值，当n 分别取100、1000、10000
时，结果是多少？ n=input('please input n:'); s=0;
for m=1:n
s=s+1.0/(m*m);
t=sqrt(6*s);
end
t
please input n:100
t =
3.1321
please input n:1000
t =
3.1406
please input n:10000
t =
3.1415
4．一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。试输出全部水仙花数。
for a=1:9;
for b=0:9;
for c=0:9;
i=100*a+10*b+c;
if a^3+b^3+c^3==i;
disp(i);
end
end end end 153 370 371 407
5．定义三个函数文件，分别求给定双曲正弦、余弦和正切函数（如下），并在命令文件中调用该函数文件。
s i n h ()2
x x e e x --= c o s h ()2x x e e x -+= t a n h
()x x
x x e e x e e ---=+ 并用你定义的函数分别绘制上述三条曲线。 >> function [imean stdev]=stat(x) mean=1/2*(exp(x)-exp(-x)); stdev=sinh(x); mean=stdev;
function [imean stdev]=stat(x) mean=1/2*(exp(x)+exp(-x)); stdev=sinh(x); mean=stdev;
>> function [imean stdev]=stat(x)
mean=(exp(x)-exp(-x))/(exp(x)+exp(-x)); stdev=sinh(x); mean=stdev;
6．已知y1＝x2，y2＝cos(2x)，y3=y1*y2，完成下列操作：(1)在同——坐标系下用不同的颜色和线型绘制3条曲线； x=0:2*pi;
y1=x.^2;
plot(x,y1,'r--');
hold on;
y2=cos(2*x);
plot(x,y2,'g-.');
hold on;
y3=y1.*y2;
plot(x,y3);
hold off;
(2)以子图形式绘制3条曲线。x=0:2*pi;
y1=x.^2;
plot(x,y1) ;
x=0:2*pi;
y2=cos(2*x); plot(x,y2)
x=0:2*pi;
y1=x.^2;
y2=cos(2*x); y3=y1.*y2;
plot(x,y3);
(3)分别用条形图、阶梯图、杆图和填充图绘制3条曲线。 x=0:2*pi;
y1=x.^2;
stairs(x,y1,'r--');
hold on;
y2=cos(2*x);
stairs(x,y2,'g-.');
hold on;
y3=y1.*y2;
stairs(x,y3);
hold off;
x=0:2*pi;
y1=x.^2;
stem(x,y1,'r--');
hold on;
y2=cos(2*x);
stem(x,y2,'g-.');
hold on;
y3=y1.*y2;
stem(x,y3);
hold off;
x=0:2*pi;
y1=x.^2;
hist(x,y1,'r--'); hold on;
y2=cos(2*x);
hist(x,y2,'g-.'); hold on;
y3=y1.*y2;
hist(x,y3);
hold off;
7．绘制极坐标曲线P＝asin(b+nθ)，并分析参数a，b，n对曲线形状的影响. t = 0 : .01 : 2 * pi;
polar(t,2*sin(4+5*t),'--r')
t = 0 : .01 : 2 * pi;
polar(t,22*sin(44+55*t),'--r')
8．分别用plot和fplot函数绘制函数y＝sin(1/x)的曲线，分析两曲线的差别。alpha=1./x;
y=sin(alpha);
plot(x,y)
x=0:1:50;
y=sin(1./x);
figure(2)
fplot(@(x) sin(1./x),[0,50])
习题三Matlab在工程中的应用
一、习题目的
1、掌握生成特殊矩阵和矩阵分析的方法；
2、掌握线性方程组的求解方法；
3、了解统计和分析的方法以及多项式的常用运算；
4、掌握求数值导数、数值积分、常微分方程数值求解、非线性代数方程数组求
解的；
5、掌握符号表达式的运算法则、符号矩阵运算、符号函数极限及导数、符号函
数定积分和不定积分的方法。
二、习题内容
1．分别用3种不同的数值方法解线性方程组。
>> A=[6 5 -2 5;9 -1 4 -1;3 4 2 -2;3 -9 0 2];
>> B=[-4 13 1 11]';
>> x=A\B
x =
0.6667
-1.0000
1.5000
-0.0000
>> A=[6 5 -2 5;9 -1 4 -1;3 4 2 -2;3 -9 0 2];
B=[-4 13 1 11]';
C=inv(A);
x=C*B
x =
0.6667
-1.0000
1.5000
-0.0000
>>
[x,y,z,u]=solve('6*x+5*y-2*z+5*u+4','9*x-y+4*z-u-13','3*x+4*y+2*z-2*u-1','3*x-9* y+2*u-11')
x =
y =
2/3
z =
-1
u =
3/2
2．将20个学生5门功课的成绩存入矩阵P中，进行如下处理：
(1) 分别求每门课的最高分、最低分及相应学生序号。
(2) 分别求每门课的平均分和标推方差。
(3) 5门课总分的最高分、最低分及相应学生序号。
(4) 将5门课总分按从大到小顺序存入zcj中，相应学生序号存入xsxh.
p=[1,67,78,65,87,89;2,98,78,67,89,78;...
3,99,78,98,67,78;4,78,98,45,67,89;...
5,64,65,78,99,86;6,70,71,73,87,88;...
7,82,73,67,56,89;8,61,47,87,56,0;...
9,46,97,67,87,56;10,34,68,56,87,88;...
11,56,78,78,74,75;12,56,67,86,76,98;...
13,78,44,56,43,56;14,98,32,56,65,16;...
15,67,76,57,87,78;16,45,56,56,67,87;...
17,67,67,56,83,84;18,87,89,89,90,69;...
19,9,78,78,66,87;20,56,78,78,78,89];
a=p(:,2);b=p(:,3);c=p(:,4);d=p(:,5);e=p(:,6);
a1=max(a)
a2=min(a)
matlab基础练习题带答案
Matlab 基础练习题 常量、变量、表达式 1、 MATLAB 中，下面哪些变量名是合法的？（ ） （A ）_num （B ）num_ （C ）num- （D ）-num 2、 在MA TLAB 中，要给出一个复数z 的模，应该使用（ ）函数。 （A ）mod(z) （B ）abs(z) （C ）double(z) （D ）angle(z) 3、 下面属于MATLAB 的预定义特殊变量的是？（ ） （A ）eps （B ）none （C ）zero （D ）exp 4、 判断：在MA TLAB 的内存工作区中，存放一个英文字符 'a' 需要占用1个字节，存放 一个中文字符‘啊’需要占用2个字节。（ 错，都是2个字节 ） 5、 判断：MA TLAB 中，i 和j （ 对 ） 6、 判断：MA TLAB 中，pi 代表圆周率，它等于3.14。（ 错，后面还有很多位小数 ） 7、 在MA TLAB 中，若想计算的5 1)3.0sin(21+= πy 值，那么应该在MA TLAB 的指令窗中 输入的MA TLAB 指令是__y1=2*sin(0.3*pi)/(1+sqrt(5))_。 8、 在MA TLAB 中，a = 1，b = i ，则a 占_8__个字节，b 占_16_个字节，c 占________字 节。 9、 在MA TLAB 中，inf 的含义是__无穷大__，nan 的含义是__非数（结果不定）___。 数组 1、 在MA TLAB 中，X 是一个一维数值数组，现在要把数组X 中的所有元素按原来次序 的逆序排列输出，应该使用下面的（ ）指令。 （A ）X[end:1] （B ）X[end:-1:1] （C ）X (end:-1:1) （D ）X(end:1) 2、 在MA TLAB 中，A 是一个字二维数组，要获取A 的行数和列数，应该使用的MATLAB 的命令是（ ）。 （A ）class(A) （B ）sizeof(A) （C ）size(A) （D ）isa(A) 3、 在MATLAB 中，用指令x=1:9生成数组x 。现在要把x 数组的第二和第七个元素都 赋值为0，应该在指令窗中输入（ ） （A ）x([2 7])=(0 0) （B ）x([2,7])=[0,0] （C ）x[(2,7)]=[0 0] （D ）x[(2 7)]=(0 0) 4、 在MA TLAB 中，依次执行以下指令：clear;A=ones(3,4); A(:)=[-6:5];这时， 若在指令窗中输入指令b=A(:,2)'，那么，MATLAB 输出的结果应该是（ ） （A ）b = -3 -2 -1 （B ）b = -2 -1 0 1 （C ）b = -5 -1 3 （D ）b = -5 -2 1 4 5、 在MA TLAB 中，A = 1:9，现在执行如下指令L1 = ~(A>5)，则MATLAB 的执行结果应 该是L1 =___ 1 1 1 1 1 0 0 0 0___。
matlab课后习题解答第二章
第2章符号运算 习题2及解答 1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度” 对象，还是“符号”符号对象 3/7+; sym(3/7+; sym('3/7+'); vpa(sym(3/7+) 〖目的〗 不能从显示形式判断数据类型，而必须依靠class指令。 〖解答〗 c1=3/7+ c2=sym(3/7+ c3=sym('3/7+') c4=vpa(sym(3/7+) Cs1=class(c1) Cs2=class(c2) Cs3=class(c3) Cs4=class(c4) c1 = c2 = 37/70 c3 = c4 = Cs1 = double Cs2 = sym Cs3 = sym Cs4 = sym 2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认 为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') 〖目的〗 理解自由符号变量的确认规则。 〖解答〗 symvar(sym('sin(w*t)'),1) ans = w symvar(sym('a*exp(-X)'),1)
ans = a symvar(sym('z*exp(j*th)'),1) ans = z 5求符号矩阵???? ??????=3332 31 232221 131211 a a a a a a a a a A 的行列式值和逆，所得结果应采用“子表达式置换”简洁化。 〖目的〗 理解subexpr 指令。 〖解答〗 A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]') DA=det(A) IA=inv(A); [IAs,d]=subexpr(IA,d) A = [ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA = a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs = [ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d*(a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d = 1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31) 8（1）通过符号计算求t t y sin )(=的导数 dt dy 。（2）然后根据此结果，求- =0t dt dy 和2 π = t dt dy 。 〖目的〗 diff, limit 指令的应用。 如何理解运行结果。 〖解答〗 syms t
matlab考试题及答案
%1、编制一个解数论问题的函数文件：取任意整数，若是偶数，则用2除，否则乘3加1，重复此过程，直到整数变为1。 function f=NO_1(X); Y(1)=X;k=1; while (X~=1) k=k+1; if (mod(X,2)==0) X=X/2; else X=3*X+1; end Y(k)=X; end plot(Y,'b.') end % 2、编制程序产生一个数组，满足：a1=1，a2=1，从第三个元素开始，每个元素等于前两个元素的和，直到数组的前后两个元素的比值比小于1e-4，并且以红色点线的形式画出这个数组。 clear; A(1)=1; A(2)=1; i=3; Z=1; while (abs(Z)>=1e-4) A(i)=A(i-1)+A(i-2); Z=A(i-1)/A(i); i=i+1; end plot (A,'r.') % 3、编写一个函数，能够产生分段函数。function y=test_3_1(X) if (X<=2) y=*X; elseif (X>6) y=; else y=调用分段函数，绘制曲线。 clear; i=1; for j=0::2; x(i)=j;y(i)=test_3_1(j)*test_3_1(j+2); i=i+1; end plot(x,y) % 4、在2pi周期内画正弦函数曲线，并加注坐标轴标识和标题，然后在3pi/4，pi，5pi/4处分别加入带箭头的说明性文本，最后加注图例。 clear; t = 0:pi/50:2*pi; n = length(t); y = sin(t); plot(t,y,'-bo','linewidth',1) xlabel('X');ylabel('Y'); title('正弦函数曲线'); text,, ' \leftarrow 3pi/4','FontSize',18); text,, ' \leftarrow 3pi/4','FontSize',18); text,,' \leftarrow 3pi/4','FontSize',18); hleg1 = legend('sin(x)'); % 5、A为任意一个n*m矩阵，写程序来计算A 中有多少个零元素，并输出个数。 A=input('输入一个矩阵 A = ') n=length(find(A==0)) % 6、A为任意一个向量，写程序找出A中的最小元素，并且输出这个最小元素。 A=input('输入一个向量 A = '); x=length(A); i=1; y=A(i); while (i
matlab课后习题及答案详解
第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 在安装MATLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？
利用MATLAB进行时域分析
自动控制原理与系统课程实验报告 实验题目：利用MATLAB进行时域分析 班级：机电1131班姓名：刘润学号：38号 一、实验目的及内容 时域分析法是一种直接在时间域中对系统进行分析的方法，具有直观、准确的优点，并且可以提供系统时间响应的全部信息。在此实验中，主要介绍时域法进行系统分析，包括一阶系统、二阶系统以及高阶系统,以及系统的性能指标。通过实验，能够快速掌握、并利用MATLAB及控制系统箱对各种复杂控制系统进行时域分析。 二、实验设备 三、实验原理 典型的二阶系统在不同的阻尼比的情况下，它们的阶跃响应输出特性的差异是很大的。若阻尼比过小，则系统的振荡加剧，超调量大幅度增加；若阻尼比过大，则系统的响应过慢，又大大增加了调整时间，下面通过此实验课题分析输出响应变化规律： 已知二阶振荡环节的传递函数为：G（s）=ωn*ωn/（s*s+2*ζ*ωn*s+ωn*ωn）， 其中ωn=0.4，ζ从0变化到2，求此系统的单位阶跃响应曲线，并分析当ζ发生变化时，二阶系统的响应有什么样的变化规律。
四、实验步骤编出程序如下图： 五、实验结果画出图表如下图：
六、结果分析 （1）当ξ=0（无阻尼）（零阻尼）时： 无阻尼时的阶跃响应为等幅振荡曲线。如图ξ=0曲线。 （2）当0<ξ<1（欠阻尼）时： 对应不同的ξ，可画出一系列阻尼振荡曲线，且ξ越小，振荡的最大振幅愈大。如图ξ=0.4曲线。 (3)当ξ=1（临界阻尼）时： 临界阻尼时的阶跃响应为单调上升曲线。如图ξ=1曲线。 （4）当ξ>1（过阻尼）时： 过阻尼时的阶跃响应也为单调上升曲线。不过其上升的斜率较临界阻尼更慢。如图ξ=1.6曲线 七、教师评语
matlab习题及答案
2. 用MATLAB 语句输入矩阵A 和B 3．假设已知矩阵A ，试给出相应的MATLAB 命令，将其全部偶数行提取出来， 赋给B 矩阵，用magic(8)A =命令生成A 矩阵，用上述命令检验一下结果是不是正确。 4．用数值方法可以求出∑=++++++==63 63622284212i i S Λ，试不采用循环的 形式求出和式的数值解。由于数值方法是采用double 形式进行计算的，难以保证有效位数字，所以结果不一定精确。试采用运算的方法求该和式的精确值。 5．选择合适的步距绘制出下面的图形。 （1）)/1sin(t ，其中)1,1(-∈t ； （2）)tan(sin )sin(tan t t -，其中),(ππ-∈t 6. 试绘制出二元函数2 2 2 2 )1(1)1(1),(y x y x y x f z +++ +-= =的三维图和三 视图 7. 试求出如下极限。 （1）x x x x 1)93(lim +∞ →； （2）1 1lim 0-+→→xy xy y x ； （3）2 2)()cos(1lim 2 2 220 0y x y x e y x y x +→→++- 8. 已知参数方程? ??-==t t t y t x sin cos cos ln ，试求出x y d d 和3 /2 2d d π=t x y 9. 假设?-=xy t t e y x f 0 d ),(2 ，试求2222 22y f y x f x f y x ??+???-?? 10. 试求出下面的极限。 （1）??????-++-+-+-∞→1)2(1 161141121lim 2222n n Λ； （2）)131211( lim 2 222π πππn n n n n n n ++++++++∞ →Λ 11. 试求出以下的曲线积分。 （1）?+l s y x d )(22，l 为曲线)sin (cos t t t a x +=，)cos (sin t t t a y -=， )20(π≤≤t 。
matlab课程设计题目
课题一： 连续时间信号和系统时域分析及MATLAB实现 课题要求： 深入研究连续时间信号和系统时域分析的理论知识。利用MATLAB强大的图形处理功能、符号运算功能以及数值计算功能，实现连续时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB实现常用连续时间信号的时域波形（通过改变参数，分析其时域特性）。 1、单位阶跃信号， 2、单位冲激信号， 3、正弦信号， 4、实指数信号， 5、虚指数信号， 6、复指数信号。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘， 4、微分， 5、积分 三、用MATLAB实现信号的时域变换（参数变化，分析波形变化） 1、反转， 2、使移（超时，延时）， 3、展缩， 4、倒相， 5、综合变化 四、用MATLAB实现信号简单的时域分解 1、信号的交直流分解， 2、信号的奇偶分解 五、用MATLAB实现连续时间系统的卷积积分的仿真波形 给出几个典型例子，对每个例子，要求画出对应波形。 六、用MATLAB实现连续时间系统的冲激响应、阶跃响应的仿真波形。 给出几个典型例子，四种调用格式。 七、利用MATLAB实现连续时间系统对正弦信号、实指数信号的零状态响应的仿真波形。 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题二： 离散时间信号和系统时域分析及MATLAB实现。 课题要求： 深入研究离散时间信号和系统时域分析的理论知识。利用MATLAB强大的图
形处理功能、符号运算功能以及数值计算功能，实现离散时间信号和系统时域分析的仿真波形。 课题内容： 一、用MATLAB绘制常用信号的时域波形（通过改变参数分析其时域特性） 1、单位序列， 2、单位阶跃序列， 3、正弦序列， 4、离散时间实指数序列， 5、离散时间虚指数序列， 6、离散时间复指数序列。 二、用MATLAB实现信号的时域运算 1、相加， 2、相乘， 3、数乘。 三、用MATLAB实现信号的时域变换（参数变化，分析波形的变化） 1、反转， 2、时移（超时，延时）， 3、展缩， 4、倒相。 四、用MATLAB实现离散时间系统卷积和仿真波形 给出几个典型例子，对每个例子要求画出e(k)，h(k)，e(i)，h(i)，h(-i)，Rzs(k)波形。 五、用MATLAB实现离散时间系统的单位响应，阶跃响应的仿真波形 给出几个典型例子，四中调用格式。 六、用MATLAB实现离散时间系统对实指数序列信号的零状态响应的仿真波形 给出几个典型例子，要求可以改变激励的参数，分析波形的变化。 课题三： 连续时间信号傅里叶级数分析及MATLAB实现。 课题要求： 深入研究连续时间信号傅里叶级数分析的理论知识，利用MATLAB强大的图形处理功能，符号运算功能以及数值计算功能，实现连续时间周期信号频域分析的仿真波形。 课题内容： 一、用MATLAB实现周期信号的傅里叶级数分解与综合 以周期矩形波信号为例，绘出包含不同谐波次数的合成波形，观察合成波形与原矩形 波形之间的关系及吉布斯现象。
matlab程序设计第三章课后习题答案
1. p138 第6题在同一坐标轴中绘制下列两条曲线并标注两曲线交叉点。 >> t=0:0.01:pi; >> x1=t; >> y1=2*x1-0.5; >> x2=sin(3*t).*cos(t); >> y2=sin(3*t).*sin(t); >> plot(x1,y1,'r-',x2,y2,'g-') >> axis([-1,2,-1.5,1]) >> hold on >> s=solve('y=2*x-0.5','x=sin(3*t)*cos(t)','y=sin(3*t)*sin(t)'); >> plot(double(s.x),double(s.y),'*'); 截图：
p366 第4题绘制极坐标曲线，并分析对曲线形状的影响。 function [ output_args ] = Untitled2( input_args ) %UNTITLED2 Summary of this function goes here % Detailed explanation goes here theta=0:0.01:2*pi; a=input('请输入a的值：'); b=input('请输入b的值：'); n=input('请输入n的值：'); rho=a*sin(b+n*theta); polar(theta,rho,'k'); end 下面以a=1，b=1，n=1的极坐标图形为基础来分析a、b、n的影响。
对a的值进行改变：对比发现a只影响半径值的整倍变化 对b的值进行改变：对比发现b的值使这个圆转换了一定的角度
对n的值进行改变：对比发现当n>=2时有如下规律 1、当n为整数时，图形变为2n个花瓣状的图形 2、当n为奇数时，图形变为n个花瓣状的图形 分别让n为2、3、4、5
matlab习题及答案2
MATLAB 基本运算 1．在MATLAB 中如何建立矩阵?? ?? ??194375，并将其赋予变量a ？>>a=[573;491] 2．在进行算术运算时，数组运算和矩阵运算各有什么要求？ 进行数组运算的两个数组必须有相同的尺寸。进行矩阵运算的两个矩阵必须满足矩阵运算规则，如矩阵a 与b 相乘（a*b ）时必须满足a 的列数等于b 的行数。 3．数组运算和矩阵运算的运算符有什么区别？ 在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b 为矩阵乘，a.*b 为数组乘。 4．计算矩阵??????????897473535与??????????638976242之和。>>a=[535;374;798]; >>b=[242;679;836]; >>a+b ans = 7 779 1413151214 5．计算??????=572396a 与?? ????=864142b 的数组乘积。>>a=[693;275]; >>b=[241;468]; >>a.*b ans = 12 36384240 6．“左除”与“右除”有什么区别？ 在通常情况下，左除x=a\b 是a*x=b 的解，右除x=b/a 是x*a=b 的解，一般情况下，a\b ≠b/a 。
7．对于B AX =，如果??????????=753467294A ，???? ??????=282637B ，求解X 。>>A=[492;764;357]; >>B=[372628]’; >>X=A\B X = -0.5118 4.0427 1.3318 8．已知：???? ??????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 >>a=[123;456;789]; >>a.^2 ans = 1 4916 253649 6481 >>a^2 ans = 30 364266 81961021261509．[]7.0802.05--=a ，在进行逻辑运算时，a 相当于什么样的逻辑量。 相当于a=[11011]。 10．在sin(x )运算中，x 是角度还是弧度？ 在sin(x)运算中，x 是弧度，MATLAB 规定所有的三角函数运算都是按弧度进行运算。
matlab课后答案完整版
ones表示1矩阵 zeros表示0矩阵 ones(4)表示4x4的1矩阵 zeros(4)表示4x4的0矩阵 zeros(4,5)表示4x5的矩阵 eye(10,10)表示10x10的单位矩阵rand(4,5)表示4x5的伴随矩阵 det(a)表示计算a的行列式 inv(a)表示计算a的逆矩阵 Jordan(a)表示求a矩阵的约当标准块rank(a)表示求矩阵a的秩 [v,d]=eig(a)对角矩阵 b=a’表示求a矩阵的转置矩阵 sqrt表示求平方根 exp表示自然指数函数 log自然对数函数 abs绝对值 第一章 一、5（1） b=[97 67 34 10;-78 75 65 5;32 5 -23 -59]; >> c=[97 67;-78 75;32 5;0 -12]; >> d=[65 5;-23 -59;54 7]; >> e=b*c e = 5271 11574 -11336 664 1978 3112 （2）a=50:1:100 二、1 、x=-74; y=-27; z=(sin(x.^2+y.^2))/(sqrt(tan(abs(x+y)))+pi) z = 2、a=::; >> b=exp*a).*sin(a+ 3、x=[2 4; 5]; y=log(x+sqrt(1+x.^2))/2 y =4、a*b表示a矩阵和b矩阵相乘 a.*b表示a矩阵和b矩阵单个元素相乘A(m,n)表示取a矩阵第m行，第n列 A(m,:)表示取a矩阵第m行的全部元素 A(:,n)表示取a矩阵的第n列全部元素 A./B表示a矩阵除以b矩阵的对应元素， B.\A等价于A./B A.^B表示两个矩阵对应元素进行乘方运算A.^2表示a中的每个元素的平方 A^2表示A*A 例：x=[1,2,3]; y=[4,5,6]; z=x.^y z= 1 3 2 729 指数可以是标量（如y=2）.底数也可以是标量（如x=2） 5、a=1+2i; >> b=3+4i; >> c=exp((pi*i)/6) c = + d=c+a*b/(a+b) d = + 第二章 二、4、（1） y=0;k=0; >> while y<3 k=k+1; y=y+1/(2*k-1); end >> display([k-1,y-1/(2*k-1)]) ans = 第三章 二1（1）
MATLAB习题参考答案
第1章 MATLAB简介 1、MA TLAB的主要特点有：①语言简洁，编程效率高。②人机界面友善，交互性好。③绘图功能强大，便于数据可视化。④学科众多、领域广泛的MATLAB工具箱。⑤源程序的开放性。 MATLAB的典型应用领域有：①自动控制②汽车③电子④仪器仪表⑤生物医学⑥信号处理⑦通信等。 2、填空题 ⑴命令窗口、命令历史窗口、当前目录窗口 ⑵查阅、保存、编辑 ⑶清除图形窗、清除命令窗口中显示内容、清除MATLAB工作空间中保存的变量。 3、如果想查看某一变量具体内容或者对其修改操作，可以在工作空间中双击该变量名称，可以打开数组编辑器，在数组编辑器中可以查看变量的具体内容，也可以对其修改。如果想删除MATLAB内存中的变量，可以在工作空间中选中该变量，然后利用工作空间窗口的菜单命令或工具条中的快捷图标进行删除。 4、1+2+3+4+5+... （+ 后面可以直接跟...，也可以在+和...中加一个空格。） 6+7+8+9 1+2+3+4+5 ... （5后面必须跟一个空格，不能直接跟...，否则报错，这在预置一个+6+7+8+9 大数组时很重要。） 第2章矩阵与数值数组 1、填空题： ⑴非数、无穷大、机器零阈值，浮点数相对精度,eps= 2.2204e-016。 ⑵全下标、单下标。 2、阅读程序题： （本题主要考察数组的寻访、赋值和简单运算，提示：带；的语句不显示结果） ⑴ans = 2 3 7 Sa = 10 20 30 A = 1 20 5 30 9 10 4 6 8 10 ⑵ A = 1 3 5 7 9 11 13 15
2 4 6 8 10 12 14 16 ans = 1 5 9 13 2 6 10 14 3 7 11 15 4 8 12 16 A = 0 0 5 7 0 0 13 15 2 4 0 0 10 12 0 0 ⑶ ans = -1 -4 6 4 ans = 3 0 5 -2 ans = 3 6 9 12 ans = 3 6 9 12 3、A=magic(4)； L=A<10 L = 0 1 1 0 1 0 0 1 1 1 1 0 1 0 0 1 4、省略。 第3章字符串、元胞和构架数组 1. 直接创建法：S1=['Where there is life,' 'there is hope. '] %注意第2行要加入6个空格使其与第一 行字符数相等，否则报错。 S2=str2mat('Where there is life,','there is hope.') S3=strvcat('Where there is life,','there is hope.') %后两种方法则不用考虑两行 字符的数目 2. 填空题：A（2，3）；A{2，3} 3. 阅读程序题： ⑴subch = ABc12
MATLAB习题及答案 (入门)
习题 1， 计算?? ?? ??=572396a 与??? ???=864142b 的数组乘积。 a.*b ans = 12 36 3 8 42 40 2， 对于B AX =，如果??????????=753467294A ，???? ??????=282637B ，求解X 。 inv(a)*b ans = -0.5118 4.0427 1.3318 3， 已知：??? ? ? ?????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 a.*a ans = 1 4 9 16 25 36 49 64 81 >> a^2 ans = 30 36 42 66 81 96 102 126 150
4， 角度[]604530=x ，求x 的正弦、余弦、正切和余切。(应用sin,cos,tan.cot) >> sin(x./pi) ans = -0.1242 0.9826 0.2465 >> cos(x./pi) ans = -0.9923 -0.1857 0.9692 >> tan(x./pi) ans = 0.1252 -5.2915 0.2543 >> cot(x./pi) ans = 7.9894 -0.1890 3.9321 5， 将矩阵?? ????=7524a 、??????=3817b 和??? ???=2695c 组合成两个新矩阵： （1）组合成一个4?3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵 元素，第三列为按列顺序排列的c 矩阵元素，即 ?? ??? ???? ???237 912685 574 [reshape(a,4,1),reshape(b,4,1),reshape(c,4,1)] ans = 4 7 5 5 8 6 2 1 9 7 3 2 （2）按照a 、b 、c 的列顺序组合成一个行矢量，即
MATLAB习题及参考答案经典.doc
习题： 1， 计算?? ????=572396a 与??? ???=864142b 的数组乘积。 2， 对于B AX =，如果???? ? ?????=753467294A ，??????????=282637B ，求解X 。 3， 已知：?? ?? ? ?????=987654321a ，分别计算a 的数组平方和矩阵平方，并观察其结果。 4， 角度[]604530=x ，求x 的正弦、余弦、正切和余切。(应用sin,cos,tan.cot) 5， 将矩阵?? ?? ??=7524a 、??????=3817b 和??? ???=2695c 组合成两个新矩阵： （1）组合成一个4?3的矩阵，第一列为按列顺序排列的a 矩阵元素，第二列为按列顺序排列的b 矩阵元素，第三列为按列顺序排列的c 矩阵元素，即 ?? ??? ???? ???237 912685 574 （2）按照a 、b 、c 的列顺序组合成一个行矢量，即 []296531877254 6， 将(x -6)(x -3)(x -8)展开为系数多项式的形式。(应用poly,polyvalm) 7， 求解多项式x 3-7x 2+2x +40的根。(应用roots) 8， 求解在x =8时多项式(x -1)(x -2) (x -3)(x -4)的值。(应用poly,polyvalm) 9， 计算多项式9514124234++--x x x x 的微分和积分。(应用polyder,polyint ，poly2sym)
10， 解方程组???? ? ?????=??????????66136221143092x 。(应用x=a\b) 11， 求欠定方程组?? ? ???=???? ??5865394742x 的最小范数解。(应用pinv) 12， 矩阵???? ? ?????-=943457624a ，计算a 的行列式和逆矩阵。(应用det,inv) 13， y =sin(x )，x 从0到2π，?x =0.02π，求y 的最大值、最小值、均值和标准差。(应用max,min,mean,std) 14， 参照课件中例题的方法，计算表达式() 2 2 e 1053y x y x z ---=的梯度并绘图。(应用meshgrid, gradient, contour, hold on, quiver) 15， 用符号函数法求解方程a t 2+b*t +c=0。(应用solve) 16， 用符号计算验证三角等式：(应用syms,simple) 17， 求矩阵?? ? ? ??=2221 1211a a a a A 的行列式值、逆和特征根。(应用syms,det,inv,eig) 18， 因式分解：6555234-++-x x x x (应用syms, factor) 19， ? ??? ?? ?? =)sin()log(12 x x e x x a f ax ，用符号微分求df/dx 。(应用syms,diff) 20， 符号函数绘图法绘制函数x=sin(3t)cos(t)，y=sin(3t)sin(t)的图形，t 的变化范围为[0,2π]。(应用syms,ezplot) 21， 绘制曲线13++=x x y ，x 的取值范围为[-5,5]。(应用plot) 22， 有一组测量数据满足-at e =y ，t 的变化范围为0~10，用不同的线型和标记点画出a=0.1、a=0.2和a=0.5三种情况下的曲线,在图中添加标题-at e =y ，并用箭头线标识出各曲线a 的取值，并添加标题-at e =y 和图例框。(应用plot,title,text,legend) 23 24， x= [66 49 71 56 38]，绘制饼图，并将第五个切块分离出来。
最新MATLAB基础教程-薛山第二版-课后习题答案
《MATLAB及应用》实验指导书《MATLAB及应用》实验指导书 班级：T1243-7 姓名：柏元强 学号：20120430724 总评成绩： 汽车工程学院 电测与汽车数字应用中心
目录 实验04051001 MATLAB语言基础 (1) 实验04051002 MATLAB科学计算及绘图 (18) 实验04051003 MATLAB综合实例编程 (31)
实验04051001 MATLAB语言基础 1实验目的 1)熟悉MATLAB的运行环境 2)掌握MATLAB的矩阵和数组的运算 3)掌握MATLAB符号表达式的创建 4)熟悉符号方程的求解 2实验内容 第二章 1.创建double的变量，并进行计算。 (1)a=87，b=190，计算 a+b、a-b、a*b。 clear,clc a=double(87); b=double(190); a+b,a-b,a*b (2)创建 uint8 类型的变量，数值与(1)中相同，进行相同的计算。 clear,clc a=uint8(87); b=uint8(190); a+b,a-b,a*b 2．计算： (1) () sin 60
(2) e3 (3) 3cos 4??π ??? clear,clc a=sind(60) b=exp(3) c=cos(3*pi/4) 3．设2u =，3v =，计算： (1) 4 log uv v (2) () 2 2 e u v v u +- (3) clear,clc u=2;v=3; a=(4*u*v)/log(v) b=((exp(u)+v)^2)/(v^2-u) c=(sqrt(u-3*v))/(u*v) 4．计算如下表达式： (1) ()() 3542i i -+ (2) () sin 28i - clear,clc (3-5*i)*(4+2*i) sin(2-8*i) 5．判断下面语句的运算结果。 (1) 4 < 20 (2) 4 <= 20
matlab程序设计例题及答案
1.编写程序：计算1/3+2/5+3/7+……+10/21 法一： s=0; for i=1:10 s=s+i/(2*i+1); end s s = 4.4096 法二： sum((1:10)./(3:2:21)) ans = 4.4096 2.编写程序：计算1~100中即能被3整除，又能被7整除的所有数之和。 s=0; for i=1:100 if mod(i,3)==0&&mod(i,7)==0 s=s+i; end,end s s = 210 3.画出y=n！的图（1<=n<=10）,阶乘的函数自己编写，禁用MATLAB自带的阶乘函数。 x=1:10; for i=1:10 try y(i)=y(i-1)*i; catch y(i)=1; end,end plot(x,y)
12345678910 0.511.522.533.54x 10 6 4.一个数恰好等于它的因子之和，这个数就称为完数。例如，6的因子为1，2,3，而6=1+2+3，因此6就是一个完数。编程找出2000以内的所有完数。 g=[]; for n=2:2000 s=0; for r=1:n-1 if mod(n,r)==0 s=s+r; end end if s==n g=[g n]; end end g g =6 28 496
5.编写一个函数，模拟numel函数的功能，函数中调用size函数。 function y=numelnumel(x) m=size(x); y=m(1)*m(2); numelnumel([1 2 3;4 5 6]) ans = 6 6. 编写一个函数，模拟length函数的功能，函数中调用size函数。 function y=lengthlength(x) m=size(x); y=max(m(1),m(2)); lengthlength([1 2 3;4 5 6]) ans = 3 7.求矩阵rand（5）的所有元素和及各行平均值，各列平均值。 s=rand(5); sum=sum(sum(s)) mean2=mean(s,2) mean1=mean(s) sum = 13.8469
matlab例题
五、某公司投资2000万元建成一条生产线。投产后，在时刻t 的追加成本和追加收益分别为3/225)(t t t G ++=（百万元/年），3/218)(t t H -=（百万元/年）。试确定该生产线在何时停产可获最大利润？最大利润是多少？ 提示：利用函数?=T G H t R 0t 20-d ))t (-)t (()(（百万元），由于H （t ）-G （t ）单调 下降，所以H （t ）=G （t ）时，R （t ）取得最大利润。 5.解：构造函数f(t)=H(t)-G(t)=13-t-3t 2/3=0 ; 令t 1/3=x,则f(t)=-t 3-3t 2+13 可得矩阵P=[-1,-3,0,13] 求最佳生产时间的源程序如下： p=[-1,-3,0,13]; x=roots(p); t=x.^3 运行结果如下： t = 3.6768 +21.4316i 3.6768 -21.4316i 4.6465 再分别将t 的三个值带入函数f(t)，比较大小后，得到最大利润与最佳生产时间。 求最大利润的程序代码如下： ① t=3.6768 +21.4316i; x=0:0.01:t; y=13-x-3*x.^(2/3); trapz(x,y) 运行结果： ans = 25.2583 ② t=3.6768 -21.4316i; x=0:0.01:t; y=13-x-3*x.^(2/3); trapz(x,y) 运行结果： ans = 25.2583 ③ t=4.6465; x=0:0.01:t; y=13-x-3*x.^(2/3); trapz(x,y) 运行结果： ans = 26.3208 比较以上三组数据，可知最佳生产时间t=4.6465年，可获得的最大利润 26.3208（百万元/年）。 clear; close; fplot('18-t^(2/3)',[0,20]);grid on;hold on; fplot('5+t+2*t^(2/3)',[0,20],'r');hold off; %发现t 约为4
matlab课后习题答案
第4章数值运算 习题 4 及解答 1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分diff 或数值梯度gradient指令计算)(t y'曲线绘制 y'，然后把)(t y和)(t 在同一张图上，观察数值求导的后果。（模拟数据从prob_data401.mat获得） 〖目的〗 ●强调：要非常慎用数值导数计算。 ●练习mat数据文件中数据的获取。 ●实验数据求导的后果 ●把两条曲线绘制在同一图上的一种方法。 〖解答〗 （1）从数据文件获得数据的指令 假如prob_data401.mat文件在当前目录或搜索路径上 clear load prob_data401.mat （2）用diff求导的指令 dt=t(2)-t(1); yc=diff(y)/dt; %注意yc的长度将比y短1 plot(t,y,'b',t(2:end),yc,'r') （3）用gradent求导的指令（图形与上相似） dt=t(2)-t(1); yc=gradient(y)/dt;
plot(t,y,'b',t,yc,'r') grid on 〖说明〗 ● 不到万不得已，不要进行数值求导。 ● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级 以上。 ● 求导会使数据中原有的噪声放大。 2 采用数值计算方法，画出dt t t x y x ? =0 sin )(在]10 ,0[区间曲线，并计算)5.4(y 。 〖提示〗 ● 指定区间内的积分函数可用cumtrapz 指令给出。 ● )5.4(y 在计算要求不太高的地方可用find 指令算得。 〖目的〗 ● 指定区间内的积分函数的数值计算法和cumtrapz 指令。 ● find 指令的应用。 〖解答〗 dt=1e-4; t=0:dt:10; t=t+(t==0)*eps; f=sin(t)./t; s=cumtrapz(f)*dt; plot(t,s,'LineWidth',3) ii=find(t==4.5); s45=s(ii) s45 =
matlab练习题和答案
matlab练习题和答案 控制系统仿真实验 Matlab部分实验结果 目录 实验一 MATLAB基本操 作 ............................................................................................ 1 实验二 Matlab编程 .................................................................................................... 5 实验三Matlab底层图形控制 (6) 实验四控制系统古典分析.............................................................................................12 实验五控制系统现代分析 . (15) 实验六 PID控制器的设 计 ...........................................................................................19 实验七系统状态空间设计.............................................................................................23 实验九直流双闭环调速系统仿真 . (25) 实验一 MATLAB基本操作 1 用MATLAB可以识别的格式输入下面两个矩阵 1233，，1443678，i，，,,2357,,,,2335542，i,,,, A,1357B,,,2675342， i,,3239,,,,189543，，,,1894，， 再求出它们的乘积矩阵C，并将C矩阵的右下角2×3子矩阵赋给D矩阵。赋值完成后，调 用相应的命令查看MATLAB工作空间的占用情况。 A=[1,2,3,3;2,3,5,7;1,3,5,7;3,2,3,9;1,8,9,4];