![](https://img-blog.csdnimg.cn/20201014180756922.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
数字信号处理
学习过程积累!
「已注销」
在校大学生,主要做故障诊断,信号处理,深度学习方面,希望指导,公众号:轴承故障诊断与寿命预测,欢迎联系欢迎点赞
展开
-
matlab互相关运算
求两个数组相关性x = input('输入数组1 ');y = input('输入数组2 ');n1 = length(y)-1; n2 = length(x)-1;r = conv(x,fliplr(y));k = (-n1):n2';subplot(211)stem(k,r);xlabel('Lag index'); ylabel('Amplitude');r1=xcorr(...原创 2019-09-21 17:55:18 · 7987 阅读 · 1 评论 -
matlab平均值消噪声
原始信号s=0.9(n(0.9^n)加入随机噪声rand(n)利用多次采样取平均值消除噪声,代码如下:clcclear allclose allr=50;M=0:r-1;s=0.9*M.*(0.9.^M);k=rand(r,1)-0.5;s1=s+k';subplot(311)stem(M,s),title('原始信号')subplot(312)stem(M,s1),t...原创 2019-09-21 10:53:26 · 1683 阅读 · 0 评论 -
MATLAB滑动平均滤波
示例R = 50;d = rand(R,1)-0.5;m = 0:1:R-1;s = 2*m.*(0.9.^m);x = s + d';M = input('输入滤波重复采样数 :');subplot(211)plot(m,d,'r-',m,s,'b--',m,x,'k--')xlabel('Time index n'); ylabel('Amplitude')legend(...原创 2019-09-21 20:23:03 · 2894 阅读 · 0 评论 -
中值滤波
R = 50; a = rand(1,R)-0.4;b = round(a);m = 0:R-1;s = 2*m.*(0.9.^m);x = s + b; y = medfilt1(x,3); subplot(2,1,1)stem(m,x);axis([0 50 -1 8]);xlabel('n');ylabel('Amplitude');title('Impulse Nois...原创 2019-09-21 20:37:19 · 172 阅读 · 0 评论 -
matlab中wextend函数简单示例理解
学习短时傅里叶变换,发现对信号进行延拓时用到wextend函数学习了一下:初始函数COS函数一个周期,-pi到pi进行sym延拓代码如下:clc clear all t=-pi:pi/240:pi t2=-2*pi:pi/240:2*pi+pi/120 z=cos(t) subplot(211) plot(t,z) z2=wextend(1,'sym',z,round(length(z)/...原创 2019-09-18 17:10:23 · 5409 阅读 · 0 评论 -
中值滤波频响
h1 = ones(1,5)/5; h2 = ones(1,14)/14;[H1,w] = freqz(h1, 1, 256);[H2,w] = freqz(h2, 1, 256);m1 = abs(H1); m2 = abs(H2);subplot(211)plot(w/pi,m1,'r-',w/pi,m2,'b--');ylabel('Magnitude'); xlabel('\o...原创 2019-09-21 21:54:09 · 254 阅读 · 0 评论 -
DFT与DTFT区别
k = 0:15;w = 0:511;x = cos(2*pi*k*3/16);X = fft(x);XE = fft(x,512);plot(k/16,abs(X),'o', w/512,abs(XE))xlabel('\omega/\pi'); ylabel('Magnitude')原创 2019-09-22 09:45:07 · 504 阅读 · 0 评论 -
matlab判断IIR传递函数稳定性
num = input('输入分子: ');den = input('输入分母: ');N = max(length(num),length(den));x = 1; y0 = 0; S = 0;zi = zeros(1,N-1);for n = 1:1000 [y] = filter(num,den,x); if abs(y) < 0.000001, break, ...原创 2019-09-22 15:50:20 · 2406 阅读 · 0 评论