医学信号处理作业.pdfVIP

先天下之忧而忧，后天下之乐而乐。——范仲淹

1.设x(n)和y(n)是有限长的序列，x(n)[1,0.1,1,0.1]，y(n)[0.1,1,0.1,1]，箭头



所指位置表示n＝0的序列值，箭头右边依次是n＝1、2、3┉，箭头左边依次是n＝－

1、－2、－3┉。求这两个序列的线性相关函数。（分布使用直接计算法和表格法求解）

2.试采用傅里叶变换对一段医学信号进行处理。

要求：有原信号波形、源代码和处理结果。

肌电原始信号：

做傅里叶变换：

N=10000;

M=1;

y1=fft(x1,N);

subplot(4,1,1);

plot(f(1:N/2),y1(1:N/2));

axis([0500020]);

gridon;

先天下之忧而忧，后天下之乐而乐。——范仲淹

3.试采用频谱分析对一段医学信号进行处理。

要求：有原信号波形、源代码和处理结果。

肌电原始信号：

做频谱分析：

clear;

closeall;

%fft频率分析

a=load(EMG.txt);

y=fft(a,10000);%做10000点福利叶变换

fs=1000;

N=length(y);

mag=abs(y);

f=(0:N-1)/N*fs;

figure;

plot(f,mag);%做幅频谱

xlabel(频率);

ylabel(幅值);

title(肌电幅频N=10000);

figure;

plot(f,angle(y));%做相频谱

xlabel(频率);

ylabel(相位);

title(肌电相频N=10000);

gridon;

先天下之忧而忧，后天下之乐而乐。——范仲淹

肌电幅频谱N=10000

250

200

150

值

幅