急求用matlab画幅频特性和相频特性图。
答案:2 悬赏:30 手机版
解决时间 2021-01-18 20:16
- 提问者网友:niaiwoma
- 2021-01-18 04:24
急求用matlab画幅频特性和相频特性图。
最佳答案
- 五星知识达人网友:狂恋
- 2021-01-18 05:29
%产生峰值为1的三角波,分析其0~63次谐波的幅值谱和相位谱
clf;
Fs =128; %采样频率
T = 1/Fs; % 采样周期
N = 128; % 采样点数
t = (0:N-1)*T; % 时间,单位:S
x=zeros(N);
for n=0:N-1
b=fix((n)/(N/4));
m=n+1;
A=1/(N/4);
if b==0
x(m)=A*n;
elseif b==1||b==2
x(m)=A*(N/2-n);
elseif b==3
x(m)=A*(n-N);
end;
end;
n=0:N-1;
subplot(3,1,1)
plot(t,x);
xlabel('时间/S');
ylabel('振幅');title('时域波形');grid on;
y=fft(x,N); %对信号进行快速Fourier变换
mag=abs(y)*2/N; %求取Fourier变换的振幅;*2/N转变为真实幅值
f=n*Fs/N;
subplot(3,1,2)
plot(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('幅值谱');grid on;
p=mod(angle(y)*180/pi,360);
subplot(3,1,3)
plot(f(1:N/2),p(1:N/2)); %绘出Nyquist频率之前随频率变化的相位
xlabel('频率/Hz');
ylabel('振幅');title('相位谱');grid on;
clf;
Fs =128; %采样频率
T = 1/Fs; % 采样周期
N = 128; % 采样点数
t = (0:N-1)*T; % 时间,单位:S
x=zeros(N);
for n=0:N-1
b=fix((n)/(N/4));
m=n+1;
A=1/(N/4);
if b==0
x(m)=A*n;
elseif b==1||b==2
x(m)=A*(N/2-n);
elseif b==3
x(m)=A*(n-N);
end;
end;
n=0:N-1;
subplot(3,1,1)
plot(t,x);
xlabel('时间/S');
ylabel('振幅');title('时域波形');grid on;
y=fft(x,N); %对信号进行快速Fourier变换
mag=abs(y)*2/N; %求取Fourier变换的振幅;*2/N转变为真实幅值
f=n*Fs/N;
subplot(3,1,2)
plot(f(1:N/2),mag(1:N/2)); %绘出Nyquist频率之前随频率变化的振幅
xlabel('频率/Hz');
ylabel('振幅');title('幅值谱');grid on;
p=mod(angle(y)*180/pi,360);
subplot(3,1,3)
plot(f(1:N/2),p(1:N/2)); %绘出Nyquist频率之前随频率变化的相位
xlabel('频率/Hz');
ylabel('振幅');title('相位谱');grid on;
全部回答
- 1楼网友:三千妖杀
- 2021-01-18 07:02
参考 2.很多时候"过去"总是显得美好的原因是我们追加了太多华丽的片断进去。
我要举报
如以上回答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息,可以点下面链接进行举报!
点此我要举报以上问答信息
大家都在看
推荐资讯