MATLAB画功率谱密度

基于Matlab功率谱密度估计方法

摘要

在实际情况下, 许多平稳信号无法导出数学表达式, 要准确获取这些信号的功率谱密度存在一定的困难。根据维纳-辛钦 (Wiener Khintchine)定理,提出一种基于Matlab编程实现这类信号的功率谱密度的估计方法。通过仿真实验表明, 该方法简单易行,准确性较高。

关键词：平稳信号；功率谱密度；估计方法；

Estimation Method for Power Spectral Density of Stationary

Signals

Abstract

In practice m any mathematical expressions can no t be derived from stationary signals as there are: some difficulties in getting the power spectral density of the signals According to Wiener Khintchine theorem, the paper suggested an estimation method for power spectral d

振动台在使用中经常运用的公式

振动台在使用中经常运用的公式

1、 求推力（F）的公式

F=（m0+m1+m2+ ??）A??????????公式（1） 式中：F—推力（激振力）（N）

m0—振动台运动部分有效质量（kg） m1—辅助台面质量（kg）

m2—试件（包括夹具、安装螺钉）质量（kg）

A— 试验加速度（m/s2）

2、 加速度（A）、速度（V）、位移（D）三个振动参数的互换运算公式 2.1 A=ωv ????????????????????公式（2） 式中：A—试验加速度（m/s2）

V—试验速度（m/s） ω=2πf（角速度） 其中f为试验频率（Hz）

2.2 V=ωD×10-3

??????????????????公式（3） 式中：V和ω与“2.1”中同义

D—位移（mm0-p）单峰值

2.3 A=ω2

D×10-3 ??????????????????公式（4） 式中：A、D和ω与“2.1”，“2.2”中同义 公式（4）亦可简化为：

A=

f2250?D 式中：A和D与“2.3”中同义，但A的单位为g

1g=9.8m/s2

所以： A≈f225

功率谱密度的三种matlab实现方法

一：实验目的：

(1)掌握三种算法的概念、应用及特点； (2)了解谱估计在信号分析中的作用；

(3) 能够利用burg法对信号作谱估计，对信号的特点加以分析。 二;实验内容：

（1）简单说明三种方法的原理。

（2）用三种方法编写程序，在matlab中实现。

（3）将计算结果表示成图形的形式，给出三种情况的功率谱图。 （4）比较三种方法的特性。 （5）写出自己的心得体会。 三：实验原理： 1.周期图法：

周期图法又称直接法。它是从随机信号x(n)中截取N长的一段，把它视为能量有限x(n)真实功率谱Sx(ejw)的估计Sx(ejw)的抽样.

认为随机序列是广义平稳且各态遍历的，可以用其一个样本x(n)中的一段xN(n)来估计该随机序列的功率谱。这当然必然带来误差。由于对xN(n)采用DFT，就默认xN(n)在时域是周期的，以及xN(k)在频域是周期的。这种方法把随机序列样本x(n)看成是截得一段xN(n)的周期延拓，这也就是周期图法这个名字的来历。

2.相关法（间接法):

这种方法以相关函数为媒介来计算功率谱，所以又叫间接法。这种方法的具体步骤是：

第一步：从无限长随机序列x(n)中截取长度N的有限长序列列

通过功率谱密度函数还原出时域函数（公路谱）

频域法的核心是快速傅立叶变换，即功率谱密度函数在离散的采样点上与信号的频谱有着一个确定的关系。如果能够在功率谱密度函数上离散采样，构造出频谱，然后再对其进行傅立叶逆变换，即可得到时域的函数曲线。

下面以公路谱为例的matlab处理程序：

Gx0=256; %参考空间频率n0下的路面功率谱密度 n0=0.1; %参考空间频率n0 L=409.6; l=0.1;

N=L/l; %采样点数

n1=0.01; %空间频率范围n1--nu nu=3;

w=2; %频率指数 no=1/L; %空间频率间隔 Xk=[]; Xm=[];

n=linspace(0.01,3,N/2+1); GxC=Gx0*(n/n0).^(-w); k=0:N/2;

fik=randn(1,N)*2*pi; %产生0到2pi的均匀分布的随机序列 pg=GxC(1:N/2+1);

山西师范大学 本科毕业论文

学院：物理与信息工程学院 专业：电子信息工程

班级：1105班 学号：1152030221 姓名：罗永涛 指导教师：李竹

论文题目：二进制基带信号功率谱密度研究

年 月 日

目录

1.引言................................................................................................................... 1

2.两种基本的基带信号波形............................................................................... 2

2.1单极性波形...................................................................................

山西师范大学 本科毕业论文

学院：物理与信息工程学院 专业：电子信息工程

班级：1105班 学号：1152030221 姓名：罗永涛 指导教师：李竹

论文题目：二进制基带信号功率谱密度研究

年 月 日

目录

1.引言................................................................................................................... 1

2.两种基本的基带信号波形............................................................................... 2

2.1单极性波形...................................................................................

fastsim轮轨蠕滑力计算程序，适用于机车辆动力学计算．-

function F=Get_Fastsim(C,Cp,ab,N) úSTSIM求解蠕滑力程序 f=0.25; %许用粘着系数 a=ab(1);b=ab(2);%接触斑半径

fx=Cp(1);fy=Cp(2);fin=Cp(3); %蠕滑力 C11=C(1);C22=C(2);C23=C(3);%kalker系数 if N<=0

Fx=0;Fy=0;Mz=0;MU=0;

else %法向力小于0，蠕滑力也为0 Nx=10;Ny=10;%接触斑被分割的数目 v=0.25; E=2.0*10^11; G=E/2/(1+v); Tz=pi/2;

L1=8*a/3/C11/G;L2a=8*a/3/C22/G;L2b=pi*a*sqrt(a/b)/4/C23/G; z0=2*N/pi/a/b; dy=2/Ny; Tx=0;Ty=0;Mz=0; n1=a*fx/f/z0/L1; n2=a*fy/f/z0/L2a; fin1=a*b*fin/f/z0/L2b; fin2=a*a*fin/f/z0/L2b; for i=1:Ny

建筑照明设计标准

中华人民共和国国家标准 建筑照明设计标准

Standard for lighting design of buildings GB 50034-2004 前言

本标准系在原国家标准《民用建筑照明设计标准》GBJl33---90和《工业企业照明设计标准》GB 50034---92的基础上，总结了居住、公共和工业建筑照明经验，通过普查和重点实测调查，并参考了国内外建筑照明标准和照明节能标准经修订、合并而成。其中照明节能部分是由国家发展和改革委员会环境和资源综合利用司组织主编单位完成的。

本标准由总则、术语、一般规定、照明数量和质量、照明标准值、照明节能、照明配电及控制、照明管理与监督共八章和二个附录组成。主要规定了居住、公共和工业建筑的照明标准值、照明质量和照明功率密度。 2 术语 3 一般规定

3.1 照明方式和照明种类 3.1.1 按下列要求确定照明方式： 1 工作场所通常应设置一般照明；

2 同一场所内的不同区域有不同照度要求时，应采用分区一般照明；

3 对于部分作业面照度要求较高，只采用一般照明不合理的场所，宜采用混合照明； 4 在一个工作场所内不应只采用局部照明。 3.1.2 按下列要求确定照明种类： 1 工作场所均

%本程序采用振型分解反应谱法计算框架结构水平地震力

%采用KN.M单位

%运行本程序之前请运行CYGD1.M和CYGD2.M求解框架水平侧移刚度

%本程序未考虑扭转耦联振动，只能用于平面框架计算。求解所有振型。

%结构地震影响系数按高规3.3.8选取

%地震作用和作用效应按高规3.3.10计算

clear %清理WORKSPACE

k0=[263770 %各层框架侧移刚度

263770

263770

263770

123582];

m0=[1.904 %各层质量，重力荷载代表值/g

2.677

2.677

2.677

2.677]*1.0e3./9.8;

n1= 0.21712; %单榀框架地震力分配系数

Tg=0.35; %特征周期(按规范选取)

s=0.05; %阻尼比(按规范选取)

r=0.9; %衰减系数(按规范选取)

y1=0.02; %阻尼比调整系数1(按规范选取)

y2=1; %阻尼比调整系数2(按规范选取)

amax=0.08

disp(' ')

disp(' srs.m ver 2.0 July 3, 2006')

disp(' by Tom Irvine Email: tomirvine@aol.com') disp(' ')

disp(' This program calculates the shock response spectrum') disp(' of an acceleration time history, which is pre-loaded into Matlab.') disp(' The time history must have two columns: time(sec) & acceleration') disp(' ') %

clear t; clear y; clear yy; clear n; clear fn; clear a1; clear a2 clear b1; clear b2; clear jnum; clear THM; clear resp; clear x_pos; clear x_neg; %

iunit=input(' Enter acceleration unit: 1= G