MATLAB练习汇总Word格式.docx
《MATLAB练习汇总Word格式.docx》由会员分享，可在线阅读，更多相关《MATLAB练习汇总Word格式.docx（16页珍藏版）》请在冰豆网上搜索。
-10，10]范围内画出函数的三维图形。
练习三6
1、求代数方程
的根。
6
2．把1开五次方，并求其全部五个根。
（提示：
解x5-1=0）6
3.设方程的根为x=[-3,-5,-8,-9],求它们对应的x多项式的系数。
4.求方程f（x）=x^3+1.1x^2+0.55x+0.125=0的根6
5.求
的“商”及“余”多项式。
6.求多项式f（x）=x^3+4x^2+5x+6的导数和f'
（2）6
7.计算矩阵
与
之和。
8求解在x=8时多项式（x-1）（x-2）（x-3）（x-4）的值6
y=（8-1）*（8-2）*（8-3）*（8-4）6
9、计算多项式
的微分和积分。
10、将（x-6）（x-3）（x-8）展开为系数多项式的形式。
7
练习四7
1、绘制
冲激信号7
t=-1:
0.01:
6;
x=[t==0];
plot（t,x）7
axis（[-1,6,-0.1,1]）7
2、绘制
阶跃信号7
3、绘制
指数信号，将a=1和a=-1时的图形分别画于同一界面中，便于形成对比7
4、绘制
门函数7
5、已知某LTI系统的微分方程为7
y’’（t）+2y’（t）+100y（t）=f（t）7
其中，
，求系统的输出y（t）.7
6、绘制阶跃序列
8
7、已知某LTI系统的微分方程为8
y’’（t）+2y’（t）+100y（t）=10f（t）8
求系统的冲激响应和阶跃响应的波形.8
8、画出x（n）、h（n）、y（n）的波形8
axis（[-20,20,-0.05,1.05]）;
subplot（2,2,4）
m=linspace（-20,20,8001）;
plot（m,y）axis（[-20,20,-100,680]）;
第五次练习8
1、用MATLAB画出信号
的波形8
2、若某连续系统的输入为e（t），输出为r（t），系统的微分方程为：
9
若
求出系统的零状态响应y（t）9
3、求门函数
的傅里叶变换，并画出幅度频谱图9
4、已知系统函数
，求其零极点图。
5、已知系统的传递函数为
，求其频率特性。
6、计算图所示的系统的传递函数。
练习六10
7、设
10
（1）画出它在x=[0,4]区间内的曲线。
求出它的过零点的值。
（2）求此曲线在x轴上方及下方的第一块所围的面积的大小。
8、已知两信号
，
，求卷积
MATLAB程序如下：
9、设x=sint,y=sin（Nt+α）,10
若α=常数，令N=1,2,3,4,在四个子图中分别画出其曲线11
若N=2,取α=0，π/3，π/2，及π，四个子图中分别画出其曲线11
10、将一幅图片缩小1倍。
11
练习一
1、求下列联立方程的解
程序
A=[3,4,-7,-12;
5,-7,4,2;
1,0,8,-5;
-6,5,-2,10];
B=[4;
-3;
9;
-8];
X=A\B
程序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];
C1=A*B'
C2=A'
*B,C3=A.*B
inv（C1）,inv（C2）,inv（C3）
x=input（'
x='
）;
ifx>
=10
y=x^2+3;
elseifx>
=0
y=x^3+4*x;
else
y=x^5+x;
end
y
4、求x向量的最大值
x=[6787569943]
Max=max（x）
5、求x向量的最小值
Min=min（x）
6、求和
Sum=sum（x）
练习二
x=0:
pi*3;
y=sin（x）;
plot（x,y）
把x=0~2π间分为101点，画出以x为横坐标，y为纵坐标的曲线
解：
程序
x=linspace（0,2*pi,101）
y=cos（x）.*（0.5+3*sin（x）./（1+x.^2））;
plot（x,y）,grid
x=（1:
0.1:
10）;
%给出自变量x的定义域
y=2*sin（1+x）;
%写出函数形式
plot（x，y）
a=0:
pi/50:
2*pi;
X=cos（a）*3;
%参数方程
Y=sin（a）*2;
plot（X,Y）;
xlabel（'
x'
）,ylabel（'
y'
title（'
椭圆'
）
t=0:
10;
>
z=t;
x=t.*sin（3*t）;
y=t.*cos（3*t）;
plot3（x,y,z）
%产生单位球面的三维坐标
X=2*X0;
Y=2*Y0;
Z=2*Z0;
%产生半径为2的球面的三维坐标
surf（X0,Y0,Z0）;
%画单位球面
shadinginterp%采用插补明暗处理
holdon;
mesh（X,Y,Z）;
colormap（hot）;
holdoff%采用hot色图
hiddenoff%产生透视效果
axisequal,axisoff%不显示坐标轴
x=-10:
0.5:
y=x;
[X,Y]=meshgrid（x,y）;
s=sqrt（X.^2+Y.^2）+eps;
z=sin（s）./（s）;
surf（X,Y,z）
练习三
p=[3472912];
r=roots（p）
2．把1开五次方，并求其全部五个根。
解x5-1=0）
p=[10000-1];
程序>
a=poly（[-3,-5,-8,-9]）
a=1252238311080
即a（x）=x4+25x3+223x2+831x+1080
4.求方程f（x）=x^3+1.1x^2+0.55x+0.125=0的根
p=[11.10.550.125];
p1=conv（[1,0,2],conv（[1,4],[1,1]））;
p2=[1011];
[q,r]=deconv（p1,p2）;
cq=‘商多项式为’;
cr='
余多项式为'
;
disp（[cq,poly2str（q,'
s'
）]）,disp（[cr,poly2str（r,'
）]）
（2）
a=[1456];
e=polyder（a）
fx=poly2str（e,'
）
f=polyval（e,2）
a=[535;
374;
798];
b=[242;
679;
836];
c=a+b
8求解在x=8时多项式（x-1）（x-2）（x-3）（x-4）的值
y=（8-1）*（8-2）*（8-3）*（8-4）
>
p=[4–12–145];
pder=polyder（p）;
pders=poly2sym（pder）
pint=polyint（p）;
pints=poly2sym（pint）
pders=
12*x^2-24*x-14
pints=
x^4-4*x^3-7*x^2+5*x
x=[638];
y=poly（x）;
Y=poly2str（y）
练习四
冲激信号
t=-1:
plot（t,x）
axis（[-1,6,-0.1,1]）
阶跃信号
x=[t>
=0];
axis（[-1,6,-0.1,1.1]）
指数信号，将a=1和a=-1时的图形分别画于同一界面中，便于形成对比
a1=1;
a2=-1;
t=0:
1;
y1=exp（a1*t）;
y2=exp（a2*t）;
plot（t,y1,'
r'
holdon;
plot（t,y2,'
b'
门函数
t=-5:
5;
y=rectpuls（t,2）;
plot（t,y）;
ft'
axis（[-55-0.51.5]）;
5、已知某LTI系统的微分方程为
y’’（t）+2y’（t）+100y（t）=f（t）
，求系统的输出y（t）.
Y=dsolve（'
D2y+2*Dy+100*y-10*sin（2*pi*t）=0'
'
Dy（0）=0,y（0）=0'
t'
ezplot（Y,[0,6]）
n=[-5:
5];
x=[n>
stem（x）
7、已知某LTI系统的微分方程为
y’’（t）+2y’（t）+100y（t）=10f（t）
求系统的冲激响应和阶跃响应的波形.
ts=0;
te=5;
dt=0.01;
sys=tf（[10],[1,2,100]）;
n=[-5:
升级会员 