用matlab仿真轴承故障信号

求问用matlab仿真滚动轴承故障信号的具体方法~

在MATLAB的COMMAND框里直接敲入 x=0:pi/依5:贰*pi; y依=sin(x)； y贰=cos(x)； plot(x,y依,x,y贰) 这个是仿真正弦和余弦的 不知道符合你要求不，至于连接到示波器，你自己看连线

% 将以下程序copy到matlab编辑器中运行，或直接在工作区运行即可

fs = 20e3; % 采样频率
fn = 3e3; % 固有频率
y0 = 5; % 位移常数
g = 0.1; % 阻尼系数
T = 0.01; % 重复周期
N = 4096; % 采样点数
NT = round(fs*T); % 单周期采样点数
t = 0:1/fs:(N-1)/fs; % 采样时刻
t0 = 0:1/fs:(NT-1)/fs; % 单周期采样时刻
K = ceil(N/NT)+1; % 重复次数
y = [];
for i = 1:K
y = [y,y0*exp(-g*2*pi*fn*t0).*sin(2*pi*fn*sqrt(1-g^2)*t0)];
end
y = y(1:N);
Yf = fft(y); % 频谱
figure(1)
plot(t,y);
axis([0,inf,-4,5])
title('轴承故障仿真信号时域波形图')
xlabel('Time(s)')
ylabel('Amplitude')
figure(2)
f = 0:fs/N:fs-fs/N;
plot(f/1e3,abs(Yf));
xlabel('Frequency(KHz)');
ylabel('\itYm(\itfm)')
title('轴承故障仿真信号幅度谱图')

% 将以下程序copy到matlab编辑器中运行，或直接在工作区运行即可

fs = 20e3; % 采样频率
fn = 3e3; % 固有频率
y0 = 5; % 位移常数
g = 0.1; % 阻尼系数
T = 0.01; % 重复周期
N = 4096; % 采样点数
NT = round(fs*T); % 单周期采样点数
t = 0:1/fs:(N-1)/fs; % 采样时刻
t0 = 0:1/fs:(NT-1)/fs; % 单周期采样时刻
K = ceil(N/NT)+1; % 重复次数
y = [];
for i = 1:K
y = [y,y0*exp(-g*2*pi*fn*t0).*sin(2*pi*fn*sqrt(1-g^2)*t0)];
end
y = y(1:N);
Yf = fft(y); % 频谱
figure(1)
plot(t,y);
axis([0,inf,-4,5])
title('轴承故障仿真信号时域波形图')
xlabel('Time(s)')
ylabel('Amplitude')
figure(2)
f = 0:fs/N:fs-fs/N;
plot(f/1e3,abs(Yf));
xlabel('Frequency(KHz)');
ylabel('\itY\rm(\itf\rm)')
title('轴承故障仿真信号幅度谱图')

---

凤阳县13872808045： 使用matlab对故障轴承振动信号去噪 - ？
党舒丁香： % 将以下程序copy到matlab编辑器中运行,或直接在工作区运行即可fs = 20e3; % 采样频率 fn = 3e3; % 固有频率 y0 = 5; % 位移常数 g = 0.1; % 阻尼系数 T = 0.01; % 重复周期 N = 4096; % 采样点数 NT = round(fs*T); % 单周期采样点...

凤阳县13872808045： matlab在轴承故障诊断中在哪儿写程序 - ？
党舒丁香： 可以用MATLAB的M编辑器,一个M文件里写一个函数,调用和其他语言环境一样的,矩阵嘛,随便找本MATLAB基础的书上都有的.

凤阳县13872808045： 请问针对轴承的故障检测用振动法可以吗?有没有使用了轴承的设备的从正常到故障时用振动法检测时的相关数 - ？
党舒丁香： 单纯的轴承检测可以使用仪器,一般有两种,即:速度型和加速度型,速度型是测高中低三频的振动值,加速度型是测振动分贝值.装了轴承的设备要测出轴承的故障就比较麻烦,要使用专用的设备,这类设备由专业厂家开发,检测结果不太准确,难以作为判断依据.“振动烈度”在轴承专业中没有听说过,这应该不是专业的说法.衡量轴承振动的大小可以用“分贝”,或者用三频数值来描述.

凤阳县13872808045： matlab信号仿真 - ？
党舒丁香： 例子上都说了, 采用db5小波,选择一维离散小波变幻,在Matlab中设置fs=4096hz t=1s=x0(t)+x1(t)+xn(t)具体到高斯白噪声要有函数,xn(t)只是一个缩写的表达式,Matlab需要像x=sint这样的具体表达,之后用dwt(s,db5 4)这样分解重构,用plot作图希望能解决您的问题.

凤阳县13872808045： 利用matlab工具箱,输入L=Link([ pi/2 10 20 pi/2]),以这种方式建立的模...？
党舒丁香： 6206属于深沟球轴承.实际的深沟球轴承运行时是正游隙(即内外圈和滚珠之间有间隙)滚道曲率半径是滚子直径的0.52~0.54倍(轴承厂家不同,数值不同).建议检查各处尺寸设置是否符合实际,保持架是否有干涉.

凤阳县13872808045： MATLAB仿真信号 - ？
党舒丁香： 老乡好,请问你的仿真做出来了吗?我是这样想的,fs=4096; N=4096; Ts=1/fs; t=0:Ts:(N-1)*Ts; f0=16; % 频率 T=1/f0; A = 1; % 位移常数 NT = round(fs*T); % 单周期采样点数 t0 = 0:1/fs:(NT-1)/fs; % 单周期采样时刻 K = ceil(N/NT)+1; % 重复次数 p0 = []; for i = 1:K p0 = [p0,A*exp(-1000*t0).*sin(2*pi*600*t0)]; end p0 = p0(1:N); figure; plot(t,p0); 你看对吗?你是怎么仿真的?

凤阳县13872808045： matlab仿真显示 model error 怎么解决 - ？
党舒丁香： matlab不能仿真连续信号,只能在离散时间点上进行仿真.sample time就是取样周期,-1是和前一个元件相同 另外,虚机团上产品团购,超级便宜

凤阳县13872808045： Matlab与TruckSim联合仿真时,报错,报错如下图所示. - ？
党舒丁香： 兄弟,首先您这图实在是看不清.但是看这个报错的信息,应该是出现代数环了,用unit delay模块应该就可以解决问题.就是在你的反馈信号前加一个unit delay模块,在discrete模块库中.看小图,我猜你没用goto/from模块,所以我猜你这不是顶层的模块,反馈应该是从这个模块外部输入的,所以在sum的输入前加个unit delay试试吧.