MATLAB环境下轴心轨迹的绘制及降噪前处理方法(附代码示例)

MATLAB环境下轴心轨迹的绘制(包含降噪前处理)
程序运行环境为MATLAB R2018A,执行轴心轨迹的绘制。
轴心轨迹显示了轴心线的路径,可将来自于两个正交、共面传感器的一维时域信息组合成转子轴横向运动的二维图。
轴心轨迹是否需滤波处理视具体情况而定。
轴心轨迹/时域图是通过将一个轨迹与用于创建轨迹/时域图的两个XY时域图合并后得出的。
时域图显示在轨迹图的右侧,Y图位于X图的上方。
轨迹/时域图可用于建立轴心轨迹上特性的时间关系。
压缩包=数据+程序+参考
subplot(2,2,1),
plot(data.noisy, 'b');
grid on
title('Noisy Record ','Rotation',0,'FontSize',14);
xlabel({'Sample'},'FontSize',12);
ylabel('Amplitude (count)','FontSize',12)
subplot(2,2,2),
imagesc(data.t, dn.as, abs(dn.wnoisy));
title('Time-Frequency ','Rotation',0,'FontSize',14);
xlabel({'Time (s)'},'FontSize',12);
ylabel('Scale (a)','FontSize',12)

ID:5290722879802698

哥廷根数学学派2023

  • 7
    点赞
  • 11
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是一个基于 Kalman 滤波的 MATLAB 音频降噪示例代码,同时绘制了原始音频信号、降噪后的音频信号以及它们的时域图和频域图: ```matlab % 读取音频文件并进行预处理 [x, fs] = audioread('noisy_audio.wav'); x = x(:, 1); % 取左声道 x = x / max(abs(x)); % 归一化处理 % 定义卡尔曼滤波模型 A = 1; H = 1; Q = 0.0001; R = 0.1; P = 1; x0 = 0; % 对音频信号进行卡尔曼滤波处理 y = zeros(size(x)); for k = 2:length(x) [x0, P] = kalman_filter(x(k), x0, P, A, H, Q, R); y(k) = x0; end % 绘制原始音频信号和降噪后的音频信号的时域图和频域图 t = (0:length(x)-1) / fs; f = (0:length(x)-1) / length(x) * fs; X = abs(fft(x)); Y = abs(fft(y)); subplot(2, 2, 1), plot(t, x), xlabel('Time (s)'), ylabel('Amplitude'), title('Original Audio Signal (Time Domain)'); subplot(2, 2, 2), plot(f, X), xlabel('Frequency (Hz)'), ylabel('Magnitude'), title('Original Audio Signal (Frequency Domain)'); subplot(2, 2, 3), plot(t, y), xlabel('Time (s)'), ylabel('Amplitude'), title('Denoised Audio Signal (Time Domain)'); subplot(2, 2, 4), plot(f, Y), xlabel('Frequency (Hz)'), ylabel('Magnitude'), title('Denoised Audio Signal (Frequency Domain)'); % 反归一化和反向量化处理,并保存结果到新的音频文件中 y = y * max(abs(x)); audiowrite('denoised_audio.wav', y, fs); % 定义 Kalman 滤波函数 function [xk, Pk] = kalman_filter(zk, xk_1, Pk_1, A, H, Q, R) xk_ = A * xk_1; Pk_ = A * Pk_1 * A' + Q; Kk = Pk_ * H' / (H * Pk_ * H' + R); xk = xk_ + Kk * (zk - H * xk_); Pk = (eye(size(Kk, 1)) - Kk * H) * Pk_; end ``` 这段代码实现了对一个名为 `noisy_audio.wav` 的音频文件进行降噪处理,并绘制了原始音频信号和降噪后的音频信号的时域图和频域图。其中,时域图表示信号在时间轴上的变化情况,频域图表示信号在频率轴上的变化情况,可以用于观察信号的频率成分和噪声的频谱特征等信息。 需要注意的是,本示例代码绘制的时域图和频域图仅限于演示用途,实际应用中可能需要更加精细的绘图方式和参数设置,以得到更准确的结果。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值