信号与系统课上一些非系统性matlab函数概念

来源:互联网 发布:知乎周刊 epub 编辑:程序博客网 时间:2024/06/06 05:09

这东西还是要看文档,这里仅提供一个简短介绍


subplot(a b p) a表示行 b表示列 p表示图在左数第几个


polyfit(x y n) x y为横纵坐标 n表示要拟合的是几次函数


axis([a b c d]) a-b为横坐标范围 c-d为纵坐标范围


grid on 网格线


plot(x y o’)做点图


hold on 让下一张图打在上一张图上


sys=tf([a,b,c],[d,e])生成sys=as2+bs+c)/(ds+e)


impulse(sys)系统sys初始输入为δ(t)时的输出图像


Heaviside(t-t0) 阶跃函数,t0为移动的距离


stepfun(t,t0)Heaviside一样


rectpuls(t,width)产生一个幅值为1,宽度为width,相对于t=0左右对称的矩形波信号


squarew*t,DUTY)sinwt)类似(t为变量),幅值为+-1的周期性防波信号,DUTY表示占空比(百分比)。


mean(t) t的平均值


conj(a) a的共轭

imag(a) a的虚部

real(a) a的实部

atan(a) arctan


dirac(t) 狄拉克函数


fir1(M,wc)返回6dB截止频率为wcM阶低筒滤波器系数向量。

wcf/(fs/2)的数字频率,fs为抽样频率


Y = filter(B,A,X) 

输入X为滤波前序列,Y为滤波结果序列,B/A提供滤波器系数,B为分子, A为分母 

整个滤波过程是通过下面差分方程实现的: 

 a(1)*y(n) = b(1)*x(n) + b(2)*x(n-1) + ... + b(nb+1)*x(n-nb) - a(2)*y(n-1) - ... - a(na+1)*y(n-na) 


sound(Y,Fs)发声函数 Y是输入信号,Fs是采样频率



syms t; %t为变量

ft=sin(t);

fw=fourier(ft);%fw=fourier(fw,t,5);

ezplot(fw[,xmin,xmax]);%以上是傅里叶变换的示例


n=[a b c d];

ezplot(f,n)%(a,b)和(c,d)分别为xy的上下限



syms t,s;

fs=laplace(ft,t,s);

ft=ilaplace(fs,s,t);


Hs=tf([a b],[c d e]);

pzmap(Hs) %用来画出零极点图


clc:清除面板里的命令

clf:清除图形


MATLAB提供了专门用于求离散系统频响特性的函数freqz()

[H,w]=freqz(B,A,N)



plot的奇淫技巧:plot(x,y,’Color’,[rand(),rand(),rand()])%自动输出随机颜色


阅读全文
0 0
原创粉丝点击