英语翻译t = 0:.001:.25;x = sin(2*pi*50*t) + sin(2*pi*120*t);y = x + 2*randn(size(t));plot(y(1:50))title('Noisy time domain signal')Y = fft(y,256);Pyy = Y.*conj(Y)/256;f = 1000/256*(0:127);plot(f,Pyy(1:128))title('Power spectral density')xlabel('F

来源:学生作业帮助网 编辑:作业帮 时间:2024/05/09 13:47:33
英语翻译t = 0:.001:.25;x = sin(2*pi*50*t) + sin(2*pi*120*t);y = x + 2*randn(size(t));plot(y(1:50))title('Noisy time domain signal')Y = fft(y,256);Pyy = Y.*conj(Y)/256;f = 1000/256*(0:127);plot(f,Pyy(1:128))title('Power spectral density')xlabel('F

英语翻译t = 0:.001:.25;x = sin(2*pi*50*t) + sin(2*pi*120*t);y = x + 2*randn(size(t));plot(y(1:50))title('Noisy time domain signal')Y = fft(y,256);Pyy = Y.*conj(Y)/256;f = 1000/256*(0:127);plot(f,Pyy(1:128))title('Power spectral density')xlabel('F
英语翻译
t = 0:.001:.25;
x = sin(2*pi*50*t) + sin(2*pi*120*t);
y = x + 2*randn(size(t));
plot(y(1:50))
title('Noisy time domain signal')
Y = fft(y,256);
Pyy = Y.*conj(Y)/256;
f = 1000/256*(0:127);
plot(f,Pyy(1:128))
title('Power spectral density')
xlabel('Frequency (Hz)')

英语翻译t = 0:.001:.25;x = sin(2*pi*50*t) + sin(2*pi*120*t);y = x + 2*randn(size(t));plot(y(1:50))title('Noisy time domain signal')Y = fft(y,256);Pyy = Y.*conj(Y)/256;f = 1000/256*(0:127);plot(f,Pyy(1:128))title('Power spectral density')xlabel('F
t = 0:.001:.25;%%%产生时间序列t,隐含采样率f=1/0.001=1000Hz
x = sin(2*pi*50*t) + sin(2*pi*120*t);%%产生时域信号,有50Hz和120Hz两个频率
y = x + 2*randn(size(t));%%给x加噪声,randn的用法不用说了吧
plot(y(1:50))%%作图
title('Noisy time domain signal')
Y = fft(y,256);%%%对时域加噪信号做256点傅里叶变换
Pyy = Y.*conj(Y)/256;%%求个频率点处幅值
f = 1000/256*(0:127);%%频率轴
plot(f,Pyy(1:128))
title('Power spectral density')
xlabel('Frequency (Hz)')
你的频率分辨率df=1000/256=3.9Hz,不合适,无法分辨出50Hz和120Hz,望注意.
希望对你有所帮助,有问题欢迎追问,满意请采纳.
给你一个参考程序:
Fs=50;%%采样频率
N=100;%%采样点数
dt=1/Fs;%%时域最小间隔,即时域分辨率
t=(0:N-1)*dt;%%采样时间长度
df=Fs/N;%%频域最小间隔,即频域分辨率
f=(-N/2:N/2-1)*df;
x=sin(2*pi*15*t)+2*sin(2*pi*18*t);
y=fft(x);
Y=fftshift(y);
A=abs(Y);
A=A/(N/2);%%还原真实幅值
figure(1)
subplot(121)
plot(t,x)
xlabel('时间t')
ylabel('幅值x')
grid on
subplot(122)
plot(f,A)
xlabel('频率f')
ylabel('幅值A')
grid on

英语翻译t = 0:.001:.25;x = sin(2*pi*50*t) + sin(2*pi*120*t);y = x + 2*randn(size(t));plot(y(1:50))title('Noisy time domain signal')Y = fft(y,256);Pyy = Y.*conj(Y)/256;f = 1000/256*(0:127);plot(f,Pyy(1:128))title('Power spectral density')xlabel('F 已知sinα,cosα是方程25x² 5(2t 1)x+t²+t的两根且α为锐角,求t25x²- 5(2t+1)x+t²+t=0 1.S=(10-X)X =-X的平方+10X =-(X的平方-10X+25)+25 =-(X-5)得平方+25 X=5时S最大为25 2.h=-5t+30t =-5(t-6t+9)+45 =-5(t-3)+45 t=3时h最大为4 3.BP=12-2t BQ=4t S=4t(12-2t)÷2 =-4(t的平方)+24t (0<t≤6) 因为y=-4(t的平方+6t 设函数f(x)=tx^2+2t^2*x+t^2+t+1/t-1(t>0),求f(x)的最小值h(t) matlab微分方程画图问题已知:x'[t] = 2 (Sin[t] - x[t]) /Sqrt[(Sin[t] - x[t])^2 + (Cos[t] - y[t])^2]y'[t] = 2 (Cos[t] - y[t]) /Sqrt[(Sin[t] - x[t])^2 + (Cos[t] - y[t])^2]初始x[0]=0y[0]=0画出(x,y)的轨迹 一个一阶非线性非齐次微分方程,t'-2t^2+x=0求解t(x) f(x)=xlnx,求f(x)在[t,t+a](t>0)上的最小值! f(x)=x+∫0到1(x+t)f(t)dt 求f(x) 证明F(x)=积分(a到b)(|x-t|T(t)dt)在(a,b)上是凹的T(x)在[a,b]上连续,T(x)>0 微分方程 矩阵 x'(t)=x(t)+y(t)+2t y'(t)=x(t)+y(t)-2t 线性变换相关问题有线性变换T 假设 T^n X=T(T(T(...T(X)..))) 设X为线性空间内一向量,并且T^n-1 X≠0,T^n X=0 求证 X,TX,T^2 X,T^3 X,...,T^n-1 X 线性无关 matlab X=0.5*x(t-1) + 25*x(t-1)/(1+x(t-1)^(2))中的^(2)是什么意思? ∫ 0到x tf(x-t)dt=∫ 0到x (x-t)f(t)dt 为什么? 求解联立微分方程组 x'(t)-x(t)-2y(t)=t ,y'(t)-y(t)-2x(t)=t ,其初始条件为x(0)=2,y(0)=4 matlab问题,>> t=0:0.2:3*pi;x=exp(-t).*sin(t);plot(t,x) >> t=0:0.2:3*pi;x=exp(-t)*sin(t);plot(t,x)>> t=0:0.2:3*pi;x=exp(-t).*sin(t);plot(t,x)>> t=0:0.2:3*pi;x=exp(-t)*sin(t);plot(t,x)这两个为什么一个可以执行,一个却不可以 初三的一元二次方程 1、6x²+x=0 2、0.2t² - 3t=0 3、x(2 - 3x)+(3x - 2)=01、6x²+x=02、0.2t² - 3t=03、x(2 - 3x)+(3x - 2)=0 4、3(x - 5)²=x² - 25 什么地方出错 x=sin(t) y=exp(-t)*cos(t) t=0:0.01:3*pi clc plot(t,x,'r-',t,y,'c--' 求解微分方程x''-1/t*x'+(x')^2=0