响应式网站开发的设计原则,国际新闻联播,简单网站制作代码,外贸社交网站排名一、离散序列傅里叶变化——DTFT 1、DTFT公式 2、Matlab算法实现 function[X]dtft(x,n,w,flag)%计算离散时间付里叶变换
%[X]dtft(x,n,w)
%X在w频率点上的DTFT数组
%xn点有限长度序列
%n样本位置向量
%w频率点位置向量X x * (exp(-j).^(n * w));3、DTFT一些画图代码 function …一、离散序列傅里叶变化——DTFT 1、DTFT公式 2、Matlab算法实现 function[X]dtft(x,n,w,flag)%计算离散时间付里叶变换
%[X]dtft(x,n,w)
%X在w频率点上的DTFT数组
%xn点有限长度序列
%n样本位置向量
%w频率点位置向量X x * (exp(-j).^(n * w)); 3、DTFT一些画图代码 function [] signal_write(X,w,flag)
% X:数据
% w:频率向量magXabs(X);angXangle(X);realXreal(X);imagXimag(X);if(flag 1)figure();magXabs(X);angXangle(X);realXreal(X);imagXimag(X);subplot(2,2,1);plot(w/pi,magX);gridxlabel(以pi为单位的频率);title(幅度部分);ylabel(幅度)subplot(2,2,3);plot(w/pi,angX);gridxlabel(以pi为单位的频率);title(相角部分);ylabel(弧度)subplot(2,2,2);plot(w/pi,realX);gridxlabel(以pi为单位的频率);title(实部);ylabel(实部)subplot(2,2,4);plot(w/pi,imagX);gridxlabel(以pi为单位的频率);title(虚部);ylabel(虚部)
endif(flag 2)plot(w/pi,magX);gridxlabel(以pi为单位的频率);title(幅度部分);ylabel(幅度)
endif(flag 3)plot(w/pi,angX);gridxlabel(以pi为单位的频率);title(相角部分);ylabel(弧度)
endif(flag 4)plot(w/pi,realX);gridxlabel(以pi为单位的频率);title(实部);ylabel(实部)
endif(flag 5)plot(w/pi,imagX);gridxlabel(以pi为单位的频率);title(虚部);ylabel(虚部)
end 二、求LTI系统的频率响应H %example2.4
clear all;close all;b[1];
a[1 -0.8];m0;length(b)-1;
l0:length(a)-1;%频率分点
K500;
k-2*K:1:2*K;
wpi*k/K;%构建分子和分母的傅里叶变换
numb*exp(-j*m*w); %分母
dena*exp(-j*l*w); %分子
hnum./den;magHabs(h);
angHangle(h);figure(1)
subplot(2,1,1);plot(w/pi,magH),grid,title(幅度部分)
subplot(2,1,2);plot(w/pi,angH),grid,title(相角部分)
n0:100;
xcos(0.05*pi*n);
yfilter(b,a,x);
figure(2)
subplot(2,1,1);plot(n,x),grid,title(输入信号)
subplot(2,1,2);plot(n,y),grid,title(输出信号) 三、采样与重构 Matlab代码 function [ ] caiyang(Fs,N,jt,flag)
%UNTITLED3 此处显示有关此函数的摘要
% 此处显示详细说明
% Dt 模拟时间间隔在特定精度下信号为模拟的
% t 模拟时刻序列
% n 离散时间索引
% Ts 采样周期
% Fs 采样频率
% xa 在特定精度下由离散信号逼近模拟信号
% jt 是否需要重构
% flag 5 第五题
% 6 第六题Dt0.00005; % 模拟时间间隔在特定精度下信号为模拟的
Ts1/Fs; % 采样周期
n-N:1:N; % 离散时间索引
nTsn*Ts; % 序列时刻索引
t-N*Ts:Dt:N*Ts; % 模拟时刻序列 %% 只是对应相应的作业、、
if flag 5xaexp(-1000*abs(t)); % 在特定精度下由离散信号逼近模拟信号x1exp(-1000*abs(nTs)); % Fs5000 样本/s:x1为采样后的离散时间序列
endif flag 6xasin(1000*pi*t); % 在特定精度下由离散信号逼近模拟信号x1sin(1000*pi*nTs); % Fs5000 样本/s:x1为采样后的离散时间序列 endif flag 7xa sin(20*pi*t pi/4);x1 sin(20*pi*nTs pi/4);
end%%K 500; % 对pi进行K等分:相当于对单位园2pi进行1000等分
dk pi/K; % pi 的等分步进索引
w 0 : dk : pi; % 角度步进索引
X x1 * exp(-j* n*w);% 对x1序列做离散傅立叶变换
Xr real(X);
w [-fliplr(w),w(2:end)]; % 定义w负半轴
Xr [fliplr(Xr),Xr(2:end)]; % 由于实部偶对称得到Xr的负半轴 %% 决定是否重构
if jt 1figure();% 绘出xasubplot(3,1,1)plot(t*1000,xa);hold on % 绘出x(jw)stem(n*Ts*1000,x1,r.),hold off,title(时域波形) % 绘出以pi归一化的数字频率对应的频域实部波形subplot(3,1,2);plot(w/pi,Xr);title(频域波形) subplot(3,1,3)chonggou(x1,Fs,N);
endif jt 0figure();% 绘出xasubplot(2,1,1);plot(t*1000,xa);hold on % 绘出x(jw)stem(n*Ts*1000,x1,r.),hold off,title(时域波形) % 绘出以pi归一化的数字频率对应的频域实部波形subplot(2,1,2);plot(w/pi,Xr);title(频域波形)
endif jt 2% 绘出以pi归一化的数字频率对应的频域实部波形plot(w/pi,Xr);title(频域波形) endif jt 3 figure();subplot(2,1,1);% 绘出xaplot(t*1000,xa);hold on % 绘出x(jw)plot(n*Ts*1000,x1,r.),hold off,title(时域波形) xa x1 * sinc(Fs*(ones(length(nTs),1) * t-nTs*ones(1,length(t)))); % 内插重构subplot(2,1,2);plot(t*1000,xa, k ),hold onplot(n*Ts*1000,x1,r.),hold off ,title(重构波形 )axis([-N/Fs*1000,N/Fs*1000,min(x1),max(x1)]);end 重构代码 function [ ] chonggou(x1,Fs,N)
%UNTITLED4 此处显示有关此函数的摘要
% 此处显示详细说明
% x1 抽样序列
% Fs 采样率
% t 时间轴
% Dt 离散间隔模拟信号Dt 0.00005; % 模拟时间间隔在特定精度下信号为模拟的
n -N:N;
nTs n/Fs;
t -N/Fs:Dt:N/Fs; % 模拟时刻序列
xa x1 * sinc(Fs*(ones(length(nTs),1) * t-nTs*ones(1,length(t)))); % 内插重构plot(t*1000,xa, k ),hold on
stem(nTs*1000,x1, r. ),hold off ,title(重构波形 )
axis([-N/Fs*1000,N/Fs*1000,min(x1),max(x1)]);end转载于:https://www.cnblogs.com/BlueMountain-HaggenDazs/p/4969034.html
