五种随机特性噪声分析_图文

1、随机特性噪声分析源程序代码%0、生成5种随机分布函数，长度N=1024x1=unifrnd(-1,1,1,1024; %均匀分布x2=normrnd(0,1,1,1024; %正态分布x3=exprnd(1,1,1024; %指数分布x4=raylrnd(1,1,1024; %瑞利分布x5=chi2rnd(1,1,1024; %卡方分布% 1、时域波形title(时域波形;subplot(3,2,1,plot(1:1024,x1;grid on;title(均匀分布;xlabel(t/s;ylabel(A/v;axis(0 1000 -2 2 ;subplot(3,2,2,plot(1:102

2、4,x2;grid on;title(正态分布;xlabel(t/s;ylabel(A/v;axis(0 1000 -2 5 ;subplot(3,2,3,plot(1:1024,x3;grid on;title(指数分布;xlabel(t/s;ylabel(A/v;axis(0 1000 -2 5 ;subplot(3,2,4,plot(1:1024,x4;grid on;title(瑞利分布;xlabel(t/s;ylabel(A/v;axis(0 1000 -2 5 ;subplot(3,2,5,plot(1:1024,x5;grid on;title(卡方分布;xlabel(t/s;y

3、label(A/v;axis(0 1000 -2 5 ;%2、频域分布figure;f1=fft(x1,1024;f2=fft(x2,1024;f3=fft(x3,1024;f4=fft(x4,1024;f5=fft(x5,1024;subplot(3,2,1,plot(abs(f1,axis(0 1023 0 50;grid on;title(均匀分布 幅频特性;xlabel(f/Hz; ylabel(A/Vsubplot(3,2,2,plot(abs(f2;axis(0 1023 0 50;grid on;title(正态分布;xlabel(f/Hz; ylabel(A/Vsubplot(

4、3,2,3,plot(abs(f3,axis(0 1023 0 100;grid on;title(指数分布;xlabel(f/Hz; ylabel(A/Vsubplot(3,2,4,plot(abs(f4,axis(0 1023 0 50;grid on;title(瑞利分布;xlabel(f/Hz; ylabel(A/Vsubplot(3,2,5,plot(abs(f5,axis(0 1023 0 100;grid on;title(卡方分布;xlabel(f/Hz; ylabel(A/V%3、概率密度函数figure;k1,n1=ksdensity(x1; k2,n2=ksdensity

5、(x2;k3,n3=ksdensity(x3;k4,n4=ksdensity(x4;k5,n5=ksdensity(x5;subplot(3,2,1,plot(n1,k1;grid on;title(均匀分布 概率密度;xlabel(t/s;ylabel(A/v;subplot(3,2,2,plot(n2,k2;grid on;title(正态分布;xlabel(t/s;ylabel(A/v;subplot(3,2,3,plot(n3,k3;grid on;title(指数分布; xlabel(t/s;ylabel(A/v;subplot(3,2,4,plot(n4,k4;grid on;ti

6、tle(瑞利分布; xlabel(t/s;ylabel(A/v;subplot(3,2,5,plot(n5,k5;grid on;title(卡方分布;xlabel(t/s;ylabel(A/v;y1=unifpdf(x1,-1,1;y2=normpdf(x2,0,1;y3=exppdf(x3,1;y4=raylpdf(x4,1;y5=chi2pdf(x5,1;%4、均值figure;m1=mean(x1,m2=mean(x2,m3=mean(x3,m4=mean(x4,m5=mean(x5;subplot(3,2,1,plot(1:1024,m1;grid on;title(均匀分布;xla

7、bel(t/s;ylabel(（A/v）;subplot(3,2,2,plot(1:1024,m2;grid on;title(高斯分布;xlabel(时间（t/s）;ylabel(幅度（A/v）;subplot(3,2,3,plot(1:1024,m3;grid on;title(指数分布均值;xlabel(t/s;ylabel(（A/v）;subplot(3,2,4,plot(1:1024,m4;grid on;title(瑞利分布;xlabel(t/s;ylabel(（A/v）;subplot(3,2,5,plot(1:1024,m5;grid on;title(卡方分布;xlabel(

8、t/s;ylabel(（A/v）;%5、方差figure;v1=var(x1,v2=var(x2,v3=var(x3,v4=var(x4,v5=var(x5;subplot(3,2,1,plot(1:1024,v1;grid on;title(均匀分布;xlabel(t/s;ylabel(（A/v）;subplot(3,2,2,plot(1:1024,v2;grid on;title(高斯分布方差;xlabel(t/s;ylabel(（A/v）;subplot(3,2,3,plot(1:1024,v3;grid on;title(指数分布;xlabel(t/s;ylabel(（A/v）;sub

9、plot(3,2,4,plot(1:1024,v4;grid on;title(瑞利分布;xlabel(t/s;ylabel(（A/v）;subplot(3,2,5,plot(1:1024,v5;grid on;title(卡方分布;xlabel(t/s;ylabel(（A/v）;%6、自相关函数figure;title(自相关函数图;x_c1,lags=xcorr(x1,200,unbiased;x_c2,lags=xcorr(x2,200,unbiased;x_c3,lags=xcorr(x3,200,unbiased;x_c4,lags=xcorr(x4,200,unbiased;x_c

10、5,lags=xcorr(x5,200,unbiased;subplot(3,2,1,plot(lags,x_c1;grid on;title(均匀分布;subplot(3,2,2,plot(lags,x_c2;grid on;title(正态分布;subplot(3,2,3,plot(lags,x_c3;grid on;title(指数分布;subplot(3,2,4,plot(lags,x_c4;grid on;title(瑞利分布;subplot(3,2,5,plot(lags,x_c5;grid on;title(卡方分布;%7、【功率谱密度】figure;f1=fft(x1,1024

11、;f2=fft(x2,1024;f3=fft(x3,1024;f4=fft(x4,1024;f5=fft(x5,1024;p1=mean(f1.*conj(f1/1024; p2=mean(f2.*conj(f2/1024; p3=mean(f3.*conj(f3/1024; p4=mean(f4.*conj(f4/1024; p5=mean(f5.*conj(f5/1024;subplot(3,2,1,plot(1:1024,abs(p1;grid on;title(均匀分布 功率谱密度;xlabel(频率Hz;ylabel(幅值;%axis(0 511 -5*10-5 10*10-5;su

12、bplot(3,2,2,plot(1:1024,abs(p2;grid on;title(正态分布 功率谱密度;xlabel(频率Hz; ylabel(幅值%axis(0 511 0 10;subplot(3,2,3,plot(1:1024,abs(p3;grid on;title(指数分布 功率谱密度;xlabel(频率Hz; ylabel(幅值%axis(0 511 0 10;pspansubplot(3,2,4,plot(1:1024,abs(p4; grid on;title(瑞利分布 功率谱密度;xlabel(频率Hz; ylabel(幅值%axis(0 511 0 10;subpl

13、ot(3,2,5,plot(1:1024,abs(p5;grid on;title(卡方分布 功率谱密度;xlabel(频率Hz; ylabel(幅值%axis(0 511 0 10;%8、（3）验证N增大，令N=20480figure;x11=normrnd(2,3,1,20480; %正态x_c11,lags=xcorr(x11,20480,unbiased;f_c1=fft(x11,20480;f_c2=fft(x11,10240;f_c3=fft(x11,1024;p11=mean(f_c1.*conj(f_c1/20480; p22=mean(f_c2.*conj(f_c2/1024

14、0;p33=mean(f_c3.*conj(f_c3/1024;subplot(2,2,1;plot(1:20480,abs(p11;grid on;title(N为20840正态分布功率谱密度;xlabel(频率Hz; ylabel(幅值subplot(2,2,2;plot(1:10240,abs(p22;grid on;title(N为10240 正态分布功率谱密度;xlabel(频率Hz; ylabel(幅值subplot(2,2,3;plot(1:1024,abs(p33;grid on;title(N为1024 正态分布功率谱密度;xlabel(频率Hz;ylabel(幅值%figu

15、replot(x_c11;grid on;title(正态分布自相关函数;检验检验程序% 0、生成5种随机分布函数，长度N=1024x1=unifrnd(-1,1,1,1024; %均匀分布x2=normrnd(0,1,1,1024; %正态分布x3=exprnd(1,1,1024; %指数分布x4=raylrnd(1,1,1024; %瑞利分布x5=chi2rnd(1,1,1024; %卡方分布%9、分布叠加R1=zeros(1,1024;for i=0:4z=unifrnd(0,1,1,1024;span s=s+z;pspanendppspanR2=zeros(1,1024;pspanf

16、or j=0:4pspan z1=exprnd(1,1,1024;pspan s1=s1+z1;pspanendpspanfigure;pspanj1,l1=ksdensity(s;p j2,l2=ksdensity(s1; subplot(2,2,1,hist(s;grid on;title(5个均匀分布过程叠加;subplot(2,2,2;plot(l1,j1;grid on;title(曲线拟合;subplot(2,2,3;hist(s1;grid on;title(5个指数分布过程叠加;subplot(2,2,4,plot(l2,j2;grid on;title(曲线拟合;%10、正态

17、分布计算机或不透光烟度计算结果取最后三次自由加速测量结果的算术平均值, 此数据将在仪表板或电子灯牌上显示出来。在计算均值时可以忽略与测量均值相差很大的测量值。必要时，应重新进行检测。11.1汽车车速表指示误差检测操作规范适用于使用 滚筒式汽车车速检验台（下简称检验台）对汽车车速表（下简称车速表）指示误差的检测过程控制。技术要求车速表指示误差根据GB72582004的规定，由车速表指示值V1与实际车速V2之间的下列关系式设定：0V1V检验（V2/ 10）+4 （km/h）也就是当车速表的指示值（V1）为40 km/h时，检验台速度指示仪表的指示值（V2）在32.8 km/h40 km/h范围内为

18、合格。检测操作根据电子灯牌指示，引车员将汽车以不大于5 km/h的速度正直驶入检验台，驱动轮停放于检测滚筒的中间位置，并轻轻刹车停稳车辆。计算机自动控制或操作员控制降下举升台或放松滚筒锁止机构，必要时在非驱动轮前部放入止动块（前轮驱动车使用驻车制动），将车速采样按钮交给引车员。引车员确认汽车前方及驱动轮两旁没有人员后发动汽车，慢慢加速使检验台的两滚筒旋转并观察车速表读数，当汽车车速表的指示值（V1）为40 km/h时，稳定3 s5 s后按下采样按钮，检测值将会在检验台仪表板或电子灯牌上显示出来。然后减速停车（不可猛踩制动器踏板停车），在举升器自动升起后，撤离止动块，被测车可驶离检验台，预备进入

19、下一检测工位。figure; 11.2汽车前照灯发光强度、照射位置偏移量检测操作规范适用于使用前照灯检测仪（下简称灯光仪）对车辆前照灯的光强度、远光光轴及近光转折点偏移量检测过程的控制。技术要求车辆前照灯远光光束发光强度及远、近光光束照射位置偏移量依据GB72582004的规定进行评定。车辆前照灯远光光束发光强度最小值要求见（表5）注：subplot(1,2,1,normplot(s;四灯制是指前照灯具有四个灯光光束，采用四灯制的车辆其中两只对称的灯达到两灯制的要求时视为合格。 2.允许手扶拖拉机运输机组只装用一只前照灯。光束照射位置要求a.前照灯近光照射在距离10 m的屏幕上时，乘用车前照灯

20、近光光束明暗截止线转角或中点的高度应为0.7H0.9H（H为前照灯基准中心高度，下同）；其他机动车（拖拉机运输机组除外）应为0.6H0.8H。机动车（装用一只前照灯的机动车除外）前照灯近光光束水平方向位置向左偏不允许超过170 mm，向右偏不允许超过350 mm。b.前照灯远光照射在距离10 m的屏幕上时，要求在屏幕光束中心离地高度，对乘用车为0.9H1.0H；对其他车辆为0.8Htitle(；前照灯远光光束水平位置要求，左灯向左偏不允许超过170 mm，向右偏不允许超过350 mm，右灯向左或向右偏均不允许超过350 mm。注：前照灯远光照射位置偏移值检验仅对远光光束能单独调整的前照灯进行；

21、远光光束能单独调整的前照灯是指手工或通过使用专用工具能够在不影响近光光束照射角度的情况下调整远光光束照射角度的前照灯，通常情况下远近光束一体的前照灯其远光光束照射角度不能单独进行调整。表5 前照灯远光光束发光强度最小值要求 （cd）均匀分布正态检验 subplot(1,2,2,normplot(s1; title( 指数分布正态检验 ; %11 、参数估计 muhat,sigmahat,muci,sigmaci=normfit(s; % 均值、方差、均值的 muhat1,sigmahat1,muci1,sigmaci1=normfit(s1; % 均值、方差、均值的 、假设检验 h,sig,c

22、i=ttest(s,muhat h1,sig1,ci1=ttest(s1,muhat1 实验结果 结果： 8000;机 动 车 类 型检 测 项 目新注册车在 用 车一灯制两灯制四灯制一灯制两灯制四灯制三 轮 汽 车8000600060000.95置信区间、方差0.95最高设计车速小于70 km/h的汽车置信区间1000080006000其 他 汽 车置信区间、方差0.9518000置信区间150001500012000%11拖拉机运输机组标定功率18kW80006000标定功率18kW60006000一、5000m1 =5000-5.6402e-004 检测操作根据电子灯牌指示，引车员将被检

23、车辆停于灯光仪前方规定的距离处，方向与导轨垂直。（或使被检车的纵向中心线与检测场地的行驶标志线平行，如不平行，车辆应重新停放。）被检车辆应在空载、轮胎气压正常、乘坐一名引车员的条件下进行检测，另被检车辆蓄电池电量应充足，前照灯配光镜上的灰尘应清除干净。引车员置变速器于空档，车辆电源处于充电状态，开启前照灯远光灯，灯光仪的电脑系统自动或手动搜寻前照灯光源，并依次自动检测车辆前照灯左远光灯与右远光灯发光强度及远光光束照射位置偏移值。按规定在测同一侧的两个远光灯（主远光和辅远光灯）时，应遮挡不测量的车灯，为了能有时间遮挡不测量的车灯，自动检测系统需设置遮挡切换延时时间。根据电子灯牌指示，引车员将前照

24、灯远光灯变换为近光灯，灯光仪的电脑系统自动搜寻前照灯光源，并依次自动检测车辆前照灯左、右近光灯光束照射位置偏移值。灯光仪检测完之后应返回到原位，检测结束。11.3汽车喇叭声级检测操作规范m2 = 适用于使用声级计对汽车喇叭声级的检测过程控制。技术要求汽车喇叭声级标准限值（根据GB72582004执行）,汽车喇叭声级: 在距车前2 m，离地高1.2 m处测量时，其值为90dB115dB。检测操作打开声级计的电源开关，将量程开关置于合适档位（一般置于“35130”dB的档位），将计权网络设定为“A”计权，将采样档位设定为“快”档。根据电子灯牌指示，引车员将汽车停放在距声级计2 m的位置（声级计离地

25、高为1.2 m）。根据电子灯牌指示，将喇叭按钮按下2 s3 s的时间进行测量，声级计显示的最大值即为喇叭声级的测量值，其单位为dB。11.4汽车转向轮横向侧滑量检测操作规范适用于使用汽车侧滑试验台（下简称侧滑台）对非独立悬架汽车转向轮横向侧滑量值的检测过程控制。技术要求汽车转向轮横向侧滑量值根据GB72582004的规定，对前轴采用非独立悬架的汽车，其转向轮的横向侧滑量用侧滑台检测时的侧滑量值在5 m/km之间为合格。检测操作根据电子灯牌指示，引车员将被检车辆正直居中驶近侧滑台，操纵转向盘使转向轮处于正中位置，在驱动状态下以不高于5 km/h的车速平稳通过侧滑台，车辆通过滑板过程中不得转动转向

26、盘，不得踩刹车。当转向轮完全通过滑板时，取侧滑台仪表板或电子灯牌上显示出来的最大值，即为转向轮横向侧量的检测值（m/km）。11.5汽车轮荷及整车质量检测操作规范适用于使用汽车轮重试验台（下简称轮重台）对汽车轮荷及整车质量的检测过程控制。检测操作根据电子灯牌指示，引车员将被检车辆的前轴以不大于5 km/h的速度垂直驶入轮重台的中间，并轻轻刹车，待前轴在轮重台上稳定后，轮重台左右仪表板或电子灯牌上显示的数值是该轴的左右轮荷质量。同理，根据电子灯牌指示，引车员将被检车辆的后轴以不大于5 km/h的速度垂直驶入轮重台的中间，并轻轻刹车，待后轴在轮重台上稳定后，轮重台仪表板或电子灯牌上显示的数值是该轴

27、的左右轮荷质量。对于通过式轮荷计应按照生产厂家推荐的速度通过轮荷计。汽车轴荷及质量限值根据GB72582004和GB15892004的规定。1.260911.6汽车制动性能检测操作规范适用于使用滚筒反力式汽车制动试验台（下简称滚筒制动台）、平板汽车制动试验台（以下简称平板制动台）对汽车制动性能进行台试检测过程的控制。技术要求汽车、汽车列车在制动台上测出的制动力应符合表6的要求。对空载检测制动力有质疑时，可用表中规定的满载检测制动力要求进行复检。表6 台 试 检 测 制 动 力 要 求0.95890.9930v4 =0.4342二、参数估计和假设检验其他汽车、汽车列车h =后轴机 动 车 类 型

28、制动力总和与整车重量的百分比轴制动力与轴荷1）的百分比空载满载前轴三轮汽车v3 =45602）乘用车、总质量不大于3500kg的货车6050601）202）6050601）注：）用平板制动检验台检验乘用车时应按动态轴荷计算。2）空载和满载状态下测试均应满足此要求制动力平衡要求在制动力增长全过程中同时测得的左右轮制动力差的最大值，与全过程中测得的该轴左右轮最大制动力中大者之比，对前轴不应大于20%，对后轴（及其他轴）在轴制动力不小于该轴轴荷的60%时，不应大于24%；当后轴（及其他轴）制动力小于该轴轴荷的60%时，在制动力增长全过程中同时测得的左右轮制动力差的最大值不应大于该轴荷的8%。汽车的制动协调时间要