n=100 ; %取樣的資料點
dt=0.01; %取樣的時間間隔
n=0:n-1;
t=n*dt; %時間序列
frequency=n/(n*dt); %頻率序列
h=[1 zeros(1,n-1)]; %脈衝訊號，僅第乙個值為1，其他值為0
fo=2; fl=10; %頻率漸增的函式，將用於設定訊號從2hz增加到10hz
signal=chirp(t,fo,1,fl); %在時間序列上產生頻率逐漸增強的訊號--鳥鳴訊號
xh=conv(signal,h); %將鳥鳴訊號和脈衝訊號卷積
%卷積之後再fourier變換
xh=dcf(xh(1:n),n);
figure
subplot(2,2,1);
plot(frequency,real(xh)*2/n); %顯示卷積後fourier變換的實部
xlabel('頻率/hz');
ylabel('振幅');
title('訊號卷積後fourier變換的實部');
subplot(2,2,2)
plot(frequency,imag(xh)*2/n); %顯示卷積後fourier變換的虛部
xlabel('頻率/hz');
ylabel('振幅');
title('訊號卷積後fourier變換的虛部');
%先進行fourier變換，然後再作乘積
x=dcf(signal,n); %對鳥鳴訊號作fourier變換
h=dcf(h,n); %對脈衝訊號作fourier變換
xh1=x.*h;
subplot(2,2,3);
plot(frequency,real(xh1)*2/n); %顯示卷積後fourier變換的實部
xlabel('頻率/hz');
ylabel('振幅');
title('頻率域乘積xh1的實部');
subplot(2,2,4)
plot(frequency,imag(xh1)*2/n); %顯示卷積後fourier變換的虛部
xlabel('頻率/hz');
ylabel('振幅');
title('頻率域乘積xh1的虛部');
仿射變換和射影變換 等距變換 相似變換
參考 射影變換組成了乙個群，這個群被稱為射影變換群，n n可逆實矩陣稱為一般線性群gl n 當把相差非零純量因子的矩陣都視為等同時，便得到射影對映群，記為pl n 在平面，射影變換為pl 3 射影變換在平面的變換矩陣形式如下，也就是乙個3 3的矩陣。其中當上面矩陣的最後一行為 0，0，1 時的變換就...
基變換與座標變換
近來正在看 3d math primer for graphics and game development 3d數學基礎 圖形與遊戲開發 清華大學出版社，當看到第八章 矩陣與線性變換時，不禁產生疑問 當繞x軸旋轉時，求出旋轉後的基向量矩陣rx 1 0 0 rx p q r 0 cos sin 0 ...
Z變換和反變換
在 matlab語言中有專門對訊號進行正反 z變換的函式 ztrans 和itrans 其呼叫格式分別如下 lf ztrans f 對f n 進行z 變換，其結果為 f z lf ztrans f,v 對f n 進行z 變換，其結果為 f v lf ztrans f,u,v 對f u 進行z 變換，...