例2

例2

%2017/8/21日
%工程北629
%功能：信号相乘：低通滤波，采用通带：0-5Hz的低通滤波器对原始信号的振幅谱进行滤波。
%思想：使用一个滤波的例子，说明信号相乘。  信号点乘：sigmult函数clear all; clc;%输入信号
dt = 0.02;            %采样间隔0.02s
df = 1/(6000*dt);     %信号长度为120秒的信号频率分辨率(1/信号时间长度以s为单位)
n = 0:2999;           %折叠频率之前取出3000个数据点进行操作
f = n * df;           %给出频率序列%信号处理
sig = rand(1,length(n));                             %利用随机序列产生信号振幅谱  rand函数在[0,1]均匀分布的随机序列
filt = [ones(1, 5/df),zeros(1,length(n) - 5/df)];    %理想滤波器幅频响应函数
[y,n1] = sigmult(filt, n, sig, n);                   %信号相乘(点乘) == 对sig进行滤波 %输出信号
subplot(3,1,1),plot(n*df,sig),title('原始信号振幅谱');xlabel('频率f /Hz');
subplot(3,1,2),plot(n*df,filt,'LineWidth',3),title('滤波器幅度谱');xlabel('频率f /Hz');    %绘制滤波器的幅频响应
subplot(3,1,3),plot(n1*df,y),title('滤波后幅度谱'),xlabel('频率f /Hz');%总结：该程序不是非常懂？？？？ 主要在于：信号的频率此处是自己定义的频率可以吗？？？