BÁO CÁOLAB 2
Sinh viên: Bùi Minh Sinh
Lớp: 09DT1. Nhóm: 39A.
A- LAB 1A:
I . Exercise 1: Synthesizer
1. Tạo một script file sigsynth.m.
2. Download mẫu âm thanh trumpet.mat . Play âm thanh này với tần số lấy mẫu
Fs=11025.
3. Vẽ một phần nhỏ tín hiệu (khoảng 100 điểm).
* Code Matlab
load trumpet.mat
Fs=11025; soundsc(ed1,Fs)
figure, plot(ed1(200:400)); axis tight % looks very
periodic
wavwrite(ed1,Fs,'trumpet.wav')
* Kết quả:
*Nhận xét : tín hiệu tuần hoàn với chu kì: 40*1/11025= 3.628 ms.
4. Quan sát phổ tần số:
* Code Matlab:
Y = fft(ed1, 512); % take the fft of trumpet
Ymag = abs(Y); % take the mag of Y
f = Fs * (0:256)/512; % get a meaningful axis
plot(f, Ymag(1:257)); % plot Ymag (only half the points are
needed)
xlabel('Frequency (Hz)')
ylabel('Magnitude')
axis tight
*Kết quả:
5. Dùng chức năng "data cursor" để lấy ra các điểm mà tại tần số đó biên độ lớn
Các điểm mà ở các tần số đó biên độ lớn nhất:
freq=[258.4 538.3 796.7 1055 1335 1593 1852 2110 2390
2649 ];
mag=[11.8 37.75 65.11 52.24 52.59 38.99 28.31 12.49 15.72
15.21 ];
6. Tạo hàm addcosines.m để tổng hợp lại tín hiệu từ các điểm lấy được từ bước 5.
Gồm 3 vectơ :thời gian(t), tần số(freq), biên độ(mag). Sử dụng vòng lặp để tổng
hợp các tín hiêu cosin tại các tần số khác nhau thu được ở trên.
Code matlab hàm addcosines:
function [y]=addcosines(t,freq,mag)
sigsyntht=zeros(1,length(t))
for i=1:length(freq)
sigsyntht= sigsyntht + mag(i)*cos(2*pi*freq(i)*t)
end
y=0.999*sigsyntht/max(abs(sigsyntht))
Gọi hàm addcosines:
t=0:1/Fs:3;
freq=[258.4 538.3 796.7 1055 1335 1593 1852 2110 2390
2649 ];
mag=[11.8 37.75 65.11 52.24 52.59 38.99 28.31 12.49 15.72
15.21 ];
sigsynth=addcosines(t,freq,mag)
soundsc(sigsynth,Fs)
>> figure
subplot(2,1,1)
plot(ed1(200:400));
subplot(2,1,2)
plot(sigsynth (200:400));
wavwrite(sigsynth,Fs,'sigsynth1.wav');
*Kết quả:
8.Nhận xét :Hai tín hiệu âm thanh trên nghe không giống nhau.Đồ thị của chúng
cũng không giống nhau.
9.Thực hiện lại với tần số lấy mẫu nhỏ hơn và tần số lấy mẫu lớn hơn:
>> soundsc(y,Fs/2)
>> sound(y,2*Fs)
II. Exercise 2:
1.Tạo scritp phasefun.m .Tạo tín hiệu sig1 bằng cách tổng hợp 2 tín hiệu tuần
hoàn với tần số nhỏ hơn Fs/2.Fs=8000.
2.Tạo tín hiệu sig2 tương tự như sig1 nhưng thành phần thứ 2 đã bị trễ nữa chu kì.
3.Vẽ đồ thị 2 tín hiệu, quan sát kết quả.
*Code Matlab:
Fs=8000;
t=0:1/Fs:1;
sig1=cos(2*pi*500*t)+cos(2*pi*800*t);
sig2=cos(2*pi*500*t)+cos(2*pi*800*t-pi);
figure
subplot(2,1,1);
plot(sig1(200:400));
subplot(2,1,2);
plot(sig2(200:400));
*Kết quả chạy:
>> phasefun;
*Nhận xét: 2 tín hiệu khác nhau.
4.Nghe 2 âm thanh vừa tạo ra.
soundsc(sig1,Fs)
soundsc(sig2,Fs)
*Nhận xét: 2 âm nghe rất giống nhau.
5.Thực hiện lại với thành phần tín hiệu thứ 2 bị trễ 1 khoảng thời gian khác.
sig2=cos(2*pi*500*t)+cos(2*pi*800*(t-1/3200));
soundsc(sig2,Fs);
sig2=cos(2*pi*500*t)+cos(2*pi*800*(t-1/4800));
soundsc(sig2,Fs);
*Nhận xét: Các âm thanh nghe giống nhau.
6.Taoj tín hiệu sig3.Cộng tín hiệu sig3 với chính nó bị làm trễ đi 1 khoảng thời
gian để tạo thành sig4.
*Code matlab:
sig3=cos(2*pi*1000*t);
sig4=sig3+cos(2*pi*1000*(1-1/4000));
subplot(2,1,1);
plot(sig3(200:400));
subplot(2,1,2);
plot(sig4(200:400));
soundsc(sig3,Fs);
soundsc(sig4,Fs);
*Nhận xét:Hai âm thanh nghe giống nhau.
III. Exercise 3 Gibbs phenomena:
1. Plot the magnitude and phase of the coe-cients Ck for k={-10 -9 -8….8 9 10};
*Code matlab:
k=-10:10;
fct=ck(k);
subplot(2,1,1);
plot(abs(fct));
subplot(2,1,2);
plot(phase(fct));
*Kết quả:
. lớn nhất:
freq= [25 8.4 538.3 796.7 1055 1335 1593 18 52 2110 23 90
26 49 ];
mag=[11.8 37.75 65.11 52. 24 52. 59 38.99 28 .31 12. 49 15. 72
15 .21 ];
6. Tạo hàm. addcosines:
t=0:1/Fs:3;
freq= [25 8.4 538.3 796.7 1055 1335 1593 18 52 2110 23 90
26 49 ];
mag=[11.8 37.75 65.11 52. 24 52. 59 38.99 28 .31 12. 49 15. 72
15 .21 ];
sigsynth=addcosines(t,freq,mag)
soundsc(sigsynth,Fs)
>>