求信号y的功率谱密度,画出图形
来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/11/08 15:12:43
求信号y的功率谱密度,画出图形
clear;
clc;
fs=30*10^6; %采样率1k
sweeptime=10^(-4); %扫描时间
t=0:1/fs:sweeptime;%扫描时间为0~12s
f0=5*10^6; %初始频率
f1=10*10^6; %结束频率
y=chirp(t,f0,sweeptime,f1); %产生扫频信号 频率在0~12s内从10Hz增加到120Hz
figure(1) %画第一幅图
subplot(2,1,1) %在第一个位置绘图
%plot(t(1:2000),y(1:2000)); %显示扫描信号波形
plot(y);
title('signal chirp figure') %标题
xlabel('Time(s)') %横坐标标签
ylabel('Amplitude') %纵坐标标签
clear;
clc;
fs=30*10^6; %采样率1k
sweeptime=10^(-4); %扫描时间
t=0:1/fs:sweeptime;%扫描时间为0~12s
f0=5*10^6; %初始频率
f1=10*10^6; %结束频率
y=chirp(t,f0,sweeptime,f1); %产生扫频信号 频率在0~12s内从10Hz增加到120Hz
figure(1) %画第一幅图
subplot(2,1,1) %在第一个位置绘图
%plot(t(1:2000),y(1:2000)); %显示扫描信号波形
plot(y);
title('signal chirp figure') %标题
xlabel('Time(s)') %横坐标标签
ylabel('Amplitude') %纵坐标标签
在你的程序后面加上这个
p=(abs(fft(y,4096))/4096).^2; %求功率谱
subplot(2,1,2);
plot(p);
title('功率谱');
我也是菜鸟,不对的话见笑了.
再问: 你这求出来的是功率谱吗?应该是简单的做傅里叶变换得到的频谱而已吧
再答: 我觉得有平方就是功率谱,没有平方是频谱。另外你要的是功率谱密度还是功率谱?我看了另一端程序,改成画功率谱密度的话应该是如下程序 yy=fft(y,4096); p=yy.*conj(yy)/4096; %计算功率谱密度 ff=10000*(0:2047)/4096; %计算变换后不同点对应的频率值 subplot(2,1,2); plot(ff,p(1:2048)); ylabel('幅值'); xlabel('频率(Hz)'); title('功率谱密度(周期图法)');
再问: 我要的是通过间接法得到的功率谱密度,就是对信号先求自相关,然后求傅里叶变换得到的PSD
再答: 你再看看这段函数,也是看了别人的给你改的,不知道是不是你所需要 nfft=4096; yy=xcorr(y,'unbiased'); %计算序列的自相关函数 y_fft=fft(yy,nfft); P=abs(y_fft); index=0:round(nfft/2-1); k=index*fs/nfft; PSD=10*log10(P(index+1)); subplot(2,1,2); plot(k,PSD); ylabel('幅值'); xlabel('频率(Hz)'); title('功率谱密度(间接法)'); 这个思想是对的,横轴你可以根据需要稍加改动
再问: 结果画出来的图是这样的,我感觉应该不对啊!你可以看下我的百度知道的另外一个问题
再答: 没太理解你的意思。你最终目的是想要另一个问题中图对对应的时域信号么?这两个问题的联系我没有搞清楚,这个问题中y的频谱也不是你问的另个问题中的频谱啊
再问: 另一个问题你的理解是对的。我这个问题其实也就是想找到这样一个时域信号,而提出来的。说白了,我自己也画出了功率谱,但是就是觉得图不应该这样。我的最终问题就是要这样一个时域信号
p=(abs(fft(y,4096))/4096).^2; %求功率谱
subplot(2,1,2);
plot(p);
title('功率谱');
我也是菜鸟,不对的话见笑了.
再问: 你这求出来的是功率谱吗?应该是简单的做傅里叶变换得到的频谱而已吧
再答: 我觉得有平方就是功率谱,没有平方是频谱。另外你要的是功率谱密度还是功率谱?我看了另一端程序,改成画功率谱密度的话应该是如下程序 yy=fft(y,4096); p=yy.*conj(yy)/4096; %计算功率谱密度 ff=10000*(0:2047)/4096; %计算变换后不同点对应的频率值 subplot(2,1,2); plot(ff,p(1:2048)); ylabel('幅值'); xlabel('频率(Hz)'); title('功率谱密度(周期图法)');
再问: 我要的是通过间接法得到的功率谱密度,就是对信号先求自相关,然后求傅里叶变换得到的PSD
再答: 你再看看这段函数,也是看了别人的给你改的,不知道是不是你所需要 nfft=4096; yy=xcorr(y,'unbiased'); %计算序列的自相关函数 y_fft=fft(yy,nfft); P=abs(y_fft); index=0:round(nfft/2-1); k=index*fs/nfft; PSD=10*log10(P(index+1)); subplot(2,1,2); plot(k,PSD); ylabel('幅值'); xlabel('频率(Hz)'); title('功率谱密度(间接法)'); 这个思想是对的,横轴你可以根据需要稍加改动
再问: 结果画出来的图是这样的,我感觉应该不对啊!你可以看下我的百度知道的另外一个问题
再答: 没太理解你的意思。你最终目的是想要另一个问题中图对对应的时域信号么?这两个问题的联系我没有搞清楚,这个问题中y的频谱也不是你问的另个问题中的频谱啊
再问: 另一个问题你的理解是对的。我这个问题其实也就是想找到这样一个时域信号,而提出来的。说白了,我自己也画出了功率谱,但是就是觉得图不应该这样。我的最终问题就是要这样一个时域信号
若知道一个信号的功率谱密度,此信号乘以一余弦波后的功率谱密度如何求
运用matlab进行随机信号的功率谱密度估计仿真【急求】
正弦信号的功率谱密度是什么以及物理意义
求正弦信号x(t)=Asin(wt+φ)的自相关函数和功率谱密度函数
请问一道题,具有6.5MHZ带宽的某高斯信道,若信道中信号功率与噪声功率谱密度之比为45.5MHZ,求信道容量?
关于信号的功率谱密度信号的功率谱密度与信号的幅度值有关系吗?比如说,两个类型相同的信号,只不过最大的峰值幅度不同,那么,
在MATLAB中信号分析的频谱与功率谱密度相似是怎么回事?
能量信号与功率信号为什么说能量信号的频谱密度是连续的,而功率信号的频谱密度是离散的.为什么说周期信号都是功率信号.
求matlab周期三角波信号频谱分析的代码,能画出三角波信号、幅度谱和相位谱.
为什么用matlab画出的高斯白噪声功率谱密度不是直线
rt为什么matlab画出的 高斯白噪声功率谱密度不是直线.而是如下图
求下列直线的斜率以及在y轴上的截距,并画出图形