matlab 时间序列降采样

在 MATLAB 中，可以使用 resample 函数来进行时间序列的降采样。resample 函数可以将原始时间序列的采样率降低到较低的频率。

下面是一个简单的示例代码，演示如何使用 resample 函数进行时间序列的降采样：

% 创建原始时间序列
Fs = 1000; % 原始采样率（Hz）
t = 0:1/Fs:1-1/Fs; % 时间向量
x = sin(2*pi*100*t) + sin(2*pi*200*t); % 原始信号

% 降采样为新的采样率
Fs_new = 200; % 新的采样率（Hz）
[P,Q] = rat(Fs_new/Fs); % 计算降采样的比率
x_resampled = resample(x, P, Q); % 降采样

% 绘制原始信号和降采样后的信号
figure;
subplot(2,1,1);
plot(t, x);
title('原始信号');
subplot(2,1,2);
t_new = 0:1/Fs_new:(length(x_resampled)-1)/Fs_new; % 新的时间向量
plot(t_new, x_resampled);
title('降采样后的信号');

在上面的示例中，首先创建了一个包含两个正弦信号的原始时间序列。然后，使用 resample 函数将原始信号的采样率从 1000 Hz 降采样为 200 Hz。最后，绘制了原始信号和降采样后的信号。

请注意，降采样可能会导致信号的失真或信息损失，因此在降采样之前应该仔细考虑采样率的选择和信号的特性。

相关问题

较长时间序列降采样 matlab代码

下面是一个简单的 Matlab 代码示例，可以将一个较长的时间序列进行降采样：

% 设置采样率和时间间隔
fs = 1000; % 原始采样率
fs_new = 100; % 新的采样率
dt = 1/fs; % 原始时间间隔
dt_new = 1/fs_new; % 新的时间间隔

% 生成一个示例时间序列
t = 0:dt:10; % 时间向量
x = sin(2*pi*10*t) + sin(2*pi*50*t) + randn(size(t));

% 计算降采样后的时间序列
t_new = 0:dt_new:10; % 新时间向量
x_new = interp1(t, x, t_new); % 线性插值

% 绘制原始和降采样后的时间序列
figure;
subplot(211);
plot(t, x);
xlabel('Time (s)');
ylabel('Amplitude');
title('Original Signal');
subplot(212);
plot(t_new, x_new);
xlabel('Time (s)');
ylabel('Amplitude');
title('Downsampled Signal');

上述代码中，我们首先设置了原始采样率和时间间隔，然后生成了一个示例时间序列。接着，我们计算了降采样后的时间序列，其中使用了 Matlab 中的 `interp1` 函数进行线性插值。最后，我们绘制了原始和降采样后的时间序列。

matlab 5点降采样代码

在MATLAB中，实现5点线性降采样的代码通常涉及对信号数据进行等间距的抽取，使得新的样本间隔是原始间隔的两倍。这是一个简单的示例：

% 假设你有一个长度为N的连续信号，例如时间序列 signal
% 采样频率Fs

original_samples = length(signal);
new_sampling_rate = Fs / 2; % 新的采样率是原始的一半

% 计算新的采样点索引
new_sample_indices = (0:(new_sampling_rate):length(signal)-new_sampling_rate)';

% 进行5点降采样，实际上这是平均值采样
downsampled_signal = interp1(1:length(signal), signal, new_sample_indices, 'mean', 'extrap');

% 如果你想清楚地看到降采样的效果，可以绘制原始信号和降采样后的信号
figure;
plot(original_samples, signal, 'b', 'LineWidth', 2, 'DisplayName', 'Original Signal');
hold on;
plot(new_sample_indices, downsampled_signal, 'r', 'LineWidth', 2, 'DisplayName', 'Downsampled Signal');
legend('show');
xlabel('Sample Index');
ylabel('Amplitude');