基于MATLAB的线性电路频率响应特性分析 上传

课程设计任务书

学生姓名： 专业班级： 指导教师： 工作单位： 题 目: 基于MATLAB 的线性电路频率响应特性分析

初始条件：

MATLAB 软件，微机

主要任务：

利用MATLAB 强大的图形处理功能、符号运算功能和数值计算功能，实现线性电路频率响应特性的仿真波形。

1）绘出RLC 串联电路中AuR=UR/U1的幅频特性及相频特性曲线； 2）绘出RLC 串联电路中AuC=UC/U1的幅频特性及相频特性曲线； 3）绘出RLC 串联电路中AuL=UL/U1的幅频特性及相频特性曲线； 4）设L=1H ，C=1F ，

,/10

s rad =ω

改变R 之值，观察各特性曲线的变化情况；

5）撰写MATLAB 课程设计说明书

时间安排：

学习MATLAB 语言的概况 第1天 学习MATLAB 语言的基本知识 第2、3天 学习MATLAB 语言的应用环境，调试命令，绘图能力 第4、5天 课程设计 第6-9天 答辩 第10天

指导教师签名： 年 月 日 系主任（或责任教师）签名： 年 月 日

摘要

MATLAB语言具备高效、可视化及推理能力强等特点，是目前工程界流行最广的科学计算语言。特别是在电子通信领域，MATLAB常常被用于进行电路、信号与系统、数字信号处理等多个方面的理论验证与演算求解。将MATLAB软件引入到电路分析中,大大地提高了计算精度和工作效率,为电路分析提供了一个有效的辅助工具，是电子工程人员不可或缺的辅助工具软件。

本次课程设计基于MATLAB强大的图形处理功能、符号运算功能和数值计算，着重对于线性电路中较有代表的RLC串联电路的频率响应进行分析，着重训练MATLAB在电路分析的应用，能够运用相关软件进行数学模型建立、相关参量求解、结果呈现与分析。从而达到对MATLAB软件及其程序编写方式的熟悉。

关键字：MATLAB 线性电路频率响应

Abstract

MATLAB language with high efficiency, visualization and reasoning ability and other characteristics, is the current practice of the most widely popular scientific computing language. Especially in the field of electronic communications, MATLAB is often used for circuits, signals and systems, digital signal processing and other aspects of the theory of authentication and routing solution. MATLAB software is introduced into the circuit analysis, greatly improves the accuracy and efficiency. It is an effective auxiliary circuit analysis tools. MATLAB is an indispensable auxiliary tool for electronic engineers.

This course design based on MATLAB powerful graphics capabilities, and numerical computation symbolic operation, focuses on the frequency response of RLC series circuit which represented the linear circuit analysis. Training in the application of MATLAB in circuit analysis, make us be able to use relevant software to mathematical modeling, solve the relevant parameters, present and analyze the results. After the design, we will be able to achieve the MATLAB software and its programming on the way to the familiar.

Keywords: MATLAB frequency response of linear circuits

目录

摘要 (2)

Abstract (2)

目录 (3)

1. MATLAB简介 (4)

1.1.基本功能 (4)

1.2.相关应用 (4)

1.3.特点与优势 (5)

2.电路分析 (6)

3.程序设计 (8)

3.1.程序设计步骤 (8)

3.2.程序流程图 (10)

3.3.调试结果 (10)

4.结果分析 (12)

4.1.幅频响应曲线特点 (12)

4.2.相频响应曲线特点 (13)

4.3.不同R值对于曲线的影响 (13)

5. 心得体会 (14)

6. 参考文献 (15)

7. 附录：程序源代码与实验结果 (16)

7.1.程序源代码 (16)

7.2.实验结果 (17)

1.MATLAB简介

MATLAB是矩阵实验室（Matrix Laboratory）的简称，是美国Mathworks公司出品的商业数学软件，用于算法开发、数据可视化、数据分析以及数值计算的高级技术计算语言和交互式环境，主要包括MATLAB和Simulink两大部分。

1.1.基本功能

MATLAB是由美国Mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。

MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。

MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++ ，JA V A的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进行下载就可以用。

1.2.相关应用

MATLAB 产品族可以用来进行以下各种工作：

1)数值分析

2)数值和符号计算

3)工程与科学绘图

4)控制系统的设计与仿真

5)数字图像处理技术

6)数字信号处理技术

7)通讯系统设计与仿真

8)财务与金融工程

MATLAB 的应用范围非常广，包括信号和图像处理、通讯、控制系统设计、测试和测量、财务建模和分析以及计算生物学等众多应用领域。附加的工具箱（单独提供的专用MATLAB 函数集）扩展了MATLAB 环境，以解决这些应用领域内特定类型的问题。

1.3.特点与优势

MATLAB具有以下特点：

1)高效的数值计算及符号计算功能,能使用户从繁杂的数学运算分析中解脱出来;

2)具有完备的图形处理功能,实现计算结果和编程的可视化;

3)友好的用户界面及接近数学表达式的自然化语言,使学者易于学习和掌握;

4)功能丰富的应用工具箱(如信号处理工具箱、通信工具箱等) ,为用户提供了大量方便

实用的处理工具.

MATLAB与同类产品相比同时还具有以下几点优势：

1)友好的工作平台和编程环境

2)简单易用的程序语言

3)强大的科学计算机数据处理能力

4)出色的图形处理功能

5)应用广泛的模块集合工具箱

6)实用的程序接口和发布平台

7)应用软件开发（包括用户界面）

2. 电路分析

线性电路是指由时不变线性无源元件、线性受控源和独立电源组成的电路。本次课程设计主要针对RLC 串联电路进行分析。如图所示的RLC 串联电路，

在可变频率的正选电压源S U 激励下，由于感抗、容抗随频率变动，所以，电路中的电压、电流响应亦随频率变动。电路的输入阻抗(j )Z ω可表示为：

由于串联电路中同时存在着电感L 电容C ，两者的频率特性不仅相反（感抗与ω成正比，而容抗与ω成反比），而且直接相减（电抗角差180°）。一定存在一个角频率0ω，使感抗和容抗相互完全抵消，此时电路发生谐振。此时有：

这时电阻R 是唯一控制和调节谐振峰的电路元件，从而控制谐振时的电感和电容的电压及其储能状态。L 、C 串联端口相当于短路，但其两端电压都不为零，而是模值相等且反向，相互完全抵消。这时有以下关系：

图1. RLC 串联电路

1(j )()Z R j L C

ω=+ω-

ω0ω=

0000000000011

()()()()()()()()0

X S S X L C L U j j L I j j U j j U j C R CR

U j U j U j ωω=ω-

ω=ω-ωωωω=ω+ω=00()()s U j I j R

ωω=

根据上式定义了品质因数Q

，即001def

L Q R CR ω=

==

ω 为了便于比较不同参数的RLC 串联电路的频率响应之间在性能上的差异，纵、横坐标都采用相对于谐振点的比值（倍率）作为绘制频率特性的坐标系。由于已设定

0()()S S U j U j ω=ω，所以()R U j ω、()L U j ω、()C U j ω与()S U j ω比的模值就表示任一频率时的输出电压与其谐振点的输入电压之比；而横坐标也用与谐振频率0ω的比值表示，即

0ηω=ω。这样，在1η=处谐振，都在同一个相对尺度下来比较相互频率特性的差异，这一共同的尺度也表示各谐振电路的偏谐程度，以下将用新的坐标系来分析频率响应特性。

从以上所得到的结论，可以推知不同元件两端电压作为输出变量的网络函数如下：

根据课设题目给出的条件，即1L H =、1C F =、01rad ω=，上式可化为：

通过以上三个网络函数，就可以利用MATLAB 做出所要求得的RLC 串联电路的频率响应曲线，即幅频响应曲线与相频响应曲线。

22

()1()1

1

()

()1()

()()(1)

(1)()()1

1

(1)

(1)

R R S C C S L L S U j R H j U j R j L jQ C

U j jQ H j jQ U j U j jQ H j U j jQ ηηη

ηηηηηηη

ηω=

=

=

ω+ω-

+-ω-=

=

+-=

=

+-

2()1

()

1

()(1)

()1

()

R C L R

H j R j H j j R j H j R j ω=

+ω-ωω=

ω-ω-ω

ω=

+ω-ω

3.程序设计

3.1.程序设计步骤

在MATLAB中进行程序编写完成对于RLC串联电路频率响应分析需要以下几个步骤：选取适当的频率范围与作图采样周期，确定电阻R取值，相关函数表达式的编写，图形输出参数设置。

3.1.1.频率范围与作图采样周期确定

通过对相关函数表达式的分析得知，函数的变化情况在整个频率范围内并不相同。其ω=周围变化特别明显，而在ω趋于0或无穷时，函数几乎为线性函数；因此使用在1

LOGSPACE()将ω的范围定在0.1到10之间，其取样间隔为在10的对数空间均分为1000份。

3.1.2.确定电阻R取值

本次实验，需要对不同的R取值下频率响应曲线进行比较。为了方便参数的修改，在进行代码编写时，未确定R取值，而是在调用函数时提示输入R的值。此处运用了MATLAB 的INPUT()命令，这个命令允许在程序执行到此处时，显示相关提示信息，并从键盘取的相应用户输入值。

3.1.3.相关函数表达式编写

MATLAB的表达式录入相当方便，遵循了数学，特别是矩阵运算的书写规则，易于接受与理解。需要注意的是，MATLAB是以矩阵作为运算单元，在进行一些运算，例如乘、除、乘方时，需要按照矩阵书写规则表明为“.*”“./”“.^”。结合前文中得到的表达式，代码编写为：

Hr=R./(R+j.*(w-1./w));

Hc=1./(j.*w*R-(w.^2-1));

Hl=j*w./(R+j.*(w-1./w));

3.1.

4.图形输出参数设置

为了更好的体现幅频变化关系，在幅频响应曲线的纵坐标采用对数表达法，所有图像的和坐标也采用对数表达法，因此使用了SEMILOGX()函数，并且将函数的绝对值处理为dB的表示方法。曲线共分为幅频响应、相频响应两大类；R、C、L三个部分，因此曲线排列为2行3列。考虑到为了在多次选取R值后，结果的对比方便，在输出图形后加入HOLD ON 命令，使得在同一个作图位置重新作图时不清除原有图形。XLABEL、YLABEL可以为曲线的横纵坐标加注释。具体带外如下：

subplot(2,3,1),semilogx(w,20*log10(abs(Hr))),hold on

xlabel('R的频率响应|'),ylabel('dB')

subplot(2,3,2),semilogx(w,20*log10(abs(Hc))),hold on

xlabel('C的频率响应|'),ylabel('dB')

subplot(2,3,3),semilogx(w,20*log10(abs(Hl))),hold on

xlabel('L的频率响应|'),ylabel('dB')

subplot(2,3,4),semilogx(w,angle(Hr)),hold on

xlabel('R的相频响应|')

subplot(2,3,5),semilogx(w,angle(Hc)),hold on

xlabel('C的相频响应')

subplot(2,3,6),semilogx(w,angle(Hl)),hold on

xlabel('L的相频响应')

3.2. 程序流程图

图2.程序设计流程图

3.3. 调试结果

当R=1时，曲线的输出结果为：

10

10

10

R 的幅频响应

d B

10

10

10

C 的幅频响应

d B

10

10

10

L 的幅频响应

d B

10

10

10

R 的相频响应

10

10

10

C 的相频响应

10

10

10

L 的相频响应

图3.R=1时频率响应曲线

当R=0.5时，曲线的输出结果为：

10

10

10

R 的幅频响应

d B

10

10

10

C 的幅频响应

d B

10

10

10

L 的幅频响应

d B

10

10

10

R 的相频响应

10

10

10

C 的相频响应

10

10

10

L 的相频响应

图4.R=0.5时频率响应曲线

当R=0.2时，曲线的输出结果为：

图5.R=0.2时频率响应曲线

当R=0.1时，曲线的输出结果为：

10

10

10

R 的幅频响应

d B

10

10

10

C 的幅频响应

d B

10

10

10

L 的幅频响应

d B

10

10

10

R 的相频响应

10

10

10

C 的相频响应

10

10

10

L 的相频响应

图6.R=0.1时频率响应曲线

10

10

10

R 的幅频响应

d B

10

10

10

C 的幅频响应

d B

10

10

10

L 的幅频响应

d B

10

10

10

R 的相频响应

10

10

10

C 的相频响应

10

10

10

L 的相频响应

4.结果分析

4.1.幅频响应曲线特点

从以上图中看出，R、C、L的幅频响应曲线有着许多的不同。这里要引入一个通频带的概念。工程上把幅频响应降低3dB内的频率范围叫做通频带。它们的通频带分布在频率空间内，彼此互不重叠；它们的最值的大小也不相同；曲线的变化趋势也不一样。但是它们也有一个相同点，就是曲线的最大值都在1

ω=这一点处。曲线的变化趋势的原因可以通过幅频函数对频率的导数式得到答案。在这次课程设计中，我们通过对上面的得到的曲线分析如下。

4.1.1.电阻R幅频响应曲线特点

ω=为中心出现谐振峰，随着ω向两边移动，曲线电阻R幅频响应曲线整体呈现为以1

下降，并且下降的速率放缓。实际应用中可以起到带通滤波的作用，即将输入的低频信号与高频信号滤除，只将谐振点附近的频率的信号输出，可用作频率选择器。

4.1.2.电容C幅频响应曲线特点

电容C幅频响应曲线在低频区域的输出基本为水平线，在谐振点附近形成较大的突起，然后开始迅速的近乎线性下落。实际应用中可以起到低通滤波器的作用，即将输入的高频信号滤除，只输出低频信号。

4.1.3.电感L幅频响应曲线特点

电感L的幅频响应曲线与电容的有许多相似之处，但是它的曲线刚好是电容频响曲线的翻折。其开始迅速的近乎线性上升，在谐振点附近形成较大的突起，然后在高频区域的输出基本为水平线。实际应用中可以起到高通滤波器的作用，即将输入的低频信号滤除，只输出高频信号。

4.2.相频响应曲线特点

相频响应的曲线变化趋势相同，只是相位关系有着较大的差别。以电阻R的为例。在低频区，输出电压比输入电压相位超前90°，在高频区，输出电压笔输入电压相位滞后90°。在谐振点时他们同相。电容的相频响应比电阻的滞后90°，电感的相频响应比电阻的超前90°

4.3.不同R值对于曲线的影响

正如上文中所说，电阻R是唯一控制和调节谐振峰的电路元件，从而控制谐振时的电感和电容的电压及其储能状态。因为在C、L确定时，它决定了品质因数Q。而Q又决定着曲线的变化速率，当Q越大，即R越小时，曲线变化越快，图线越陡峭。当Q>1时，电感电容两端在谐振点将出现比输入电压高Q倍的过电压，即图中的那个突起。在高压电的电路系统中（如电力系统），这种过电压非常高，可能会危及系统的安全，必须采取必要的防范措施。但在低电压的电路系统中，如无线电接收系统中，则要利用谐振时出现的过电压来获得较大的输入信号。

5.心得体会

本次MATLAB课程设计让我更熟悉地掌握该软件的功能，灵活运用MATLAB,加强对MATLAB软件强大的图形处理功能、符号运算功能和数值计算功能，掌握利用MATLAB 进行建模并做出相关函数曲线，从而运用MATLAB进行电路分析。在熟悉掌握MATLAB 程序和操作的同时培养了我的独立思考能力，钻研精神，解决问题能力和动手能力。

虽然在此之前，在《信号与系统》、《数字信号处理》等课程的实验中运用过MATLAB 软件并完成简单图像操作，但知道的只是最基本的应用，对于MATLAB在图像处理上的应用，我没有完全理解。本次课程设计中通过查阅资料，阅读程序并读写程序对MATLAB 电路分析应用有了更深的了解。我同时也认识到了MATLAB功能并不只是图形的绘制及波形的处理，有着很多方面的运用，如绘制函数，处理音频、图像数据，创建用户界面等功能，实为一个功能强大的软件。

每一次课设都会给我一种新的体验与感受，这次课程设计一样，平时都是啃课本，很单调很乏味，但是，这次课设给了我们动手和自主学习的机会，当然，在这过程中不免会遇到困难，当时也会很困惑、沮丧，但问题总要解决的，所以自己还是会去想办法，我想这也是一种收获。还有，经过了课程设计，我也深刻体会到MATLAB功能的强大以及学好MATLAB的重要性，在今后的学习生活中，我会继续学习使用它。

6.参考文献

[1] 邱关源.电路原理（五版）.北京:高等教育出版社,2006

[2] 陈怀琛，吴大正，高西全.MATLAB及在电子信息课程中的应用（第3版）.北京:

电子工业出版社,2006.3

[3] 周开利,邓春晖.MA TLAB基础及应用教程.北京:北京大学出版社，2007

[4] 张威.MATLAB6.0基础与编程与入门。西安:西安电子科技大学出版社，2007

[5] 张智星.MATLAB程序设计与应用。北京:清华大学出版社，2002

7.附录：程序源代码与实验结果

7.1.程序源代码

w=logspace(-1,1,1000); %确定频率向量的取值

R=input('R='); %确定选取的电阻R阻值Hr=R./(R+j.*(w-1./w)); %电阻R频率响应函数

Hc=1./(j.*w*R-(w.^2-1)); %电容C频率响应函数

Hl=j*w./(R+j.*(w-1./w)); %电感L频率响应函数

subplot(2,3,1),semilogx(w,20*log10(abs(Hr))),hold on xlabel('R的频率响应|'),ylabel('dB')

subplot(2,3,2),semilogx(w,20*log10(abs(Hc))),hold on xlabel('C的频率响应|'),ylabel('dB')

subplot(2,3,3),semilogx(w,20*log10(abs(Hl))),hold on xlabel('L的频率响应|'),ylabel('dB')

subplot(2,3,4),semilogx(w,angle(Hr)),hold on

xlabel('R的相频响应|')

subplot(2,3,5),semilogx(w,angle(Hc)),hold on

xlabel('C的相频响应')

subplot(2,3,6),semilogx(w,angle(Hl)),hold on

xlabel('L的相频响应')

7.2. 实验结果

10-110

10

1

-40

-30-20-10

0R 的幅频响应

d

B

10-110

10

1

-40

-20020

40C 的幅频响应

d

B

10

-110

10

1

-40

-20020

40L 的幅频响应

d

B

10-110

10

1

-2-1

012R 的相频响

应

10-110

10

1

-4-3-2-10C 的相频响

应

10

-110

10

1

01234L 的相频响

应

图7.实验结果汇总图

本科生课程设计成绩评定表

指导教师签字：

年月日

第十一章电路的频率响应 习题答案

第十一章电路的频率响应 习题 一、选择题 串联谐振电路的 Q 值越高，则 (D ) (A) 电路的选择性越差，电路的通频带越窄 (B) 电路的选择性越差，电路的通频带越宽 (C) 电路的选择性越好，电路的通频带越宽 (D ) 电路的选择性越好，电路的通频带越窄 串联电路谐振时，L 、C 储存能量的总和为 (D ) (A) W = W L + W C = 0 (B) 22 1 LI W W W C L =+= (C) 2 2 1C C L CU W W W =+= (D ) 2C C L CU W W W =+= 3.R L C 串联电路发生串联谐振时，下列说法不． 正确的是： (D ) A ．端电压一定的情况下，电流为最大值 B ．谐振角频率LC 10＝ ω C ．电阻吸收有功功率最大 D ．阻抗的模值为最大 4. RLC 串联电路在0f 时发生谐振。当电源频率增加到02f 时，电路性质呈 （B ） A. 电阻性 B . 电感性 C. 电容性 D. 视电路元件参数而定 5.下面关于RLC 串联谐振电路品质因数的说法中，不正确的是 （D ） A. 品质因数越高，电路的选择性越好 B. 品质因数高的电路对非谐振频率的电流具有较强的抵制能力 C. 品质因数等于谐振频率与带宽之比 D . 品质因数等于特性感抗电压有效值与特性容抗电压有效值之比 串联谐振电路品质因数Q=100，若U R =10V ，则电源电压Us 、电容两端电压U C 分别为 ( A ) 、1000V B. 1000V 、10V C. 100V 、1000V D. 1000V 、100V 二、判断题

1.图示电路，R << 0L，保持U S 一定，当发生谐振时，电流表的读数最小。 （×） 串联电路发生谐振时，电源输出的有功功率与无功功率均为最大。（×） 3.图示RLC串联电路，S闭合前的谐振频率与品质因数为f0与Q， S闭合后 的谐振频率与品质因数为f 0'与Q '，则 f f' =，Q < Q '。（×） 并联的交流电路中，当改变电路频率出现谐振时，则此时电路端口的阻抗值最小。（×） 4.若RLC串联谐振电路的电感增加至原来的4倍（R、C不变），则谐振角频率应变为原来的2倍。（×） 三填空题 1.图示电路，当发生串联谐振时，其谐振频率f 0= ( C M L L) 2 ( 2 1 2 1 + + π )。 2.电感L= 50mH与电容C= 20F并联，其谐振角频率 = ( 1000rad/s )；其并联谐振时的阻抗Z = ( )。 串联电路如下图所示，则电路的谐振角频率 = ( 500rad/s )，电路的品质因数Q = ( 100 )。

模态分析和频率响应分析的目的

有限元分析类型 一、nastran中的分析种类 （1）静力分析 静力分析是工程结构设计人员使用最为频繁的分析手段，主要用来求解结构在与时间无关或时间作用效果可忽略的静力载荷（如集中载荷、分布载荷、温度载荷、强制位移、惯性载荷等）作用下的响应、得出所需的节点位移、节点力、约束反力、单元内力、单元应力、应变能等。该分析同时还提供结构的重量和重心数据。 （2）屈曲分析 屈曲分析主要用于研究结构在特定载荷下的稳定性以及确定结构失稳的临界载荷，NX Nastran中的屈曲分析包括两类：线性屈曲分析和非线性屈曲分析。 （3）动力学分析 NX Nastran在结构动力学分析中有非常多的技术特点，具有其他有限元分析软件所无法比拟的强大分析功能。结构动力分析不同于静力分析，常用来确定时变载荷对整个结构或部件的影响，同时还要考虑阻尼及惯性效应的作用。 NX Nastran的主要动力学分析功能：如特征模态分析、直接复特征值分析、直接瞬态响应分析、模态瞬态响应分析、响应谱分析、模态复特征值分析、直接频率响应分析、模态频率响应分析、非线性瞬态分析、模态综合、动力灵敏度分析等可简述如下： ?正则模态分析 正则模态分析用于求解结构的固有频率和相应的振动模态，计算广义质量，正则化模态节点位移，约束力和正则化的单元力及应力，并可同时考虑刚体模态。 ?复特征值分析 复特征值分析主要用于求解具有阻尼效应的结构特征值和振型，分析过程与实特征值分析类似。此外

Nastran的复特征值计算还可考虑阻尼、质量及刚度矩阵的非对称性。 ?瞬态响应分析（时间-历程分析） 瞬态响应分析在时域内计算结构在随时间变化的载荷作用下的动力响应，分为直接瞬态响应分析和模态瞬态响应分析。两种方法均可考虑刚体位移作用。 直接瞬态响应分析 该分析给出一个结构随时间变化的载荷的响应。结构可以同时具有粘性阻尼和结构阻尼。该分析在节点自由度上直接形成耦合的微分方程并对这些方程进行数值积分，直接瞬态响应分析求出随时间变化的位移、速度、加速度和约束力以及单元应力。 模态瞬态响应分析 在此分析中，直接瞬态响应问题用上面所述的模态分析进行相同的变换，对问题的规模进行压缩，再对压缩了的方程进行数值积分，从而得出与用直接瞬态响应分析类型相同的输出结果。 ?随机振动分析 该分析考虑结构在某种统计规律分布的载荷作用下的随机响应。例如地震波，海洋波，飞机超过建筑物的气压波动，以及火箭和喷气发动机的噪音激励，通常人们只能得到按概率分布的函数，如功率谱密度（PSD）函数，激励的大小在任何时刻都不能明确给出，在这种载荷作用下结构的响应就需要用随机振动分析来计算结构的响应。NX Nastran中的PSD可输入自身或交叉谱密度，分别表示单个或多个时间历程的交叉作用的频谱特性。计算出响应功率谱密度、自相关函数及响应的RMS值等。计算过程中，NX Nastran不仅可以像其他有限元分析那样利用已知谱，而且还可自行生成用户所需的谱。 ?响应谱分析 响应谱分析（有时称为冲击谱分析）提供了一个有别于瞬态响应的分析功能，在分析中结构的激励用各个小的分量来表示，结构对于这些分量的响应则是这个结构每个模态的最大响应的组合。 ?频率响应分析 频率响应分析主要用于计算结构在周期振荡载荷作用下对每一个计算频率的动响应。计算结果分实部和虚部两部分。实部代表响应的幅度，虚部代表响应的相角。 直接频率响应分析 直接频率响应通过求解整个模型的阻尼耦合方程，得出各频率对于外载荷的响应。该类分析在频域中主要求解两类问题。第一类是求结构在一个稳定的周期性正弦外力谱的作用下的响应。结构可以具有粘性阻尼和结构阻尼，分析得到复位移、速度、加速度、约束力、单元力和单元应力。这些量可以进行正则化以获得传递函数。 第二类是求解结构在一个稳态随机载荷作用下的响应。此载荷由它的互功率谱密度定义。而结构载荷由上面所提到的传递函数来表征。分析得出位移、加速度、约束力或单元应力的自相关系数。该分析也对自功率谱进行积分而获得响应的均方根值。 模态频率响应 模态频率响应分析和随机响应分析在频域中解决的两类问题与直接频率响应分析解决相同的问题。

频响频响分析方法总结

频响频响分析方法总结-标准化文件发布号：（9556-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII

频响分析，或者叫稳态动力学分析在abaqus中包括以下三种方法： 直接稳态动力学分析（direct solution steady state dynamic analysis） 模态稳态动力学分析（mode based steady state dynamic analysis） 子空间稳态动力学分析（subspace projection steady state dynamic analysis） 1）直接稳态动力学 优点：在直接稳态动力学分析中，系统的稳态谐波响应是通过对模型的原始方程直接积分计算出来的。如果分析的对象存在非对称刚度、包含模态阻尼以外的其他阻尼或者必须考虑粘弹性材料特性（频变特性），则不能提取特征模态的情况下，可以应用直接法进行稳态响应的计算和分析。 缺点：进行直接稳态动力学分析不需要提取系统的特征模态，而是在每个频率点对整个模型进行复杂的积分运算。因此，对于具有大阻尼和频变特性的模型，应用直接法比模态分析方法精确，但是耗时较多。 2）模态稳态动力学分析 模态稳态动力学分析方法是基于模态叠加法求解系统的稳态响应。因此，在求解稳态响应之前必须先提取无阻尼系统的特征模态，也就是在说必须在step steady state dynamics,modal前加一步step frequency。另外，必须确定需要保留的特征模态，以确保能够精确描述系统的动力学特性，也就是说如果是进行0-1000hz的分析，step frequency的number of eigenvalues requested选定的阶数的模态频率必须大于1000hz，简单的作法是这里选all……，下面的maximum……填入1000。 模态稳态动力学分析的特点：相较于直接法和子空间法分析速度快，耗时最少，计算精度低于直接法和子空间法，不适合于分析具有大阻尼特性的模型，不适合于分析具有粘弹性材料（频变特性）的模型。 3）子空间稳态动力学分析 子空间稳态动力学分析的基本思想是：首先提取无阻尼、对称系统的特征模态，并选取适当的特征向量组成特征模态子空间，然后将稳态动力学方程组投影到特征模态子空间上，通过直接法求解子空间的稳态动力学方程。 我的感觉是子空间法是直接法和模态法的折中，它的特点是模型可以定义任意形式的阻尼，可以处理具有非对称刚度矩阵的模型，可以处理具有频变特性的模型，计算时间和精度也是在直接法和模态法的中间。

matlab多元线性回归模型

云南大学数学与统计学实验教学中心 实验报告 一、实验目的 1.熟悉MATLAB的运行环境. 2.学会初步建立数学模型的方法 3.运用回归分析方法来解决问题 二、实验内容 实验一：某公司出口换回成本分析 对经营同一类产品出口业务的公司进行抽样调查,被调查的13家公司,其出口换汇成本与商品流转费用率资料如下表。试分析两个变量之间的关系,并估计某家公司商品流转费用率是6.5%的出口换汇成本. 实验二：某建筑材料公司的销售量因素分析 下表数据是某建筑材料公司去年20个地区的销售量（Y，千方），推销开支、实际帐目数、同类商品

竞争数和地区销售潜力分别是影响建筑材料销售量的因素。1）试建立回归模型，且分析哪些是主要的影响因素。2）建立最优回归模型。 提示：建立一个多元线性回归模型。

三、实验环境 Windows 操作系统; MATLAB 7.0. 四、实验过程 实验一：运用回归分析在MATLAB 里实现 输入：x=[4.20 5.30 7.10 3.70 6.20 3.50 4.80 5.50 4.10 5.00 4.00 3.40 6.90]'; X=[ones(13,1) x]; Y=[1.40 1.20 1.00 1.90 1.30 2.40 1.40 1.60 2.00 1.00 1.60 1.80 1.40]'; plot(x,Y,'*'); [b,bint,r,rint,stats]=regress(Y,X,0.05); 输出： b = 2.6597 -0.2288 bint = 1.8873 3.4322 -0.3820 -0.0757 stats = 0.4958 10.8168 0.0072 0.0903 即==1,0?6597.2?ββ，-0.2288,0?β的置信区间为[1.8873 3.4322],1,?β的置信区间为[-0.3820 -0.0757]； 2r =0.4958, F=10.8168, p=0.0072 因P<0.05, 可知回归模型 y=2.6597-0.2288x 成立. 1 1.5 2 2.5 散点图 估计某家公司商品流转费用率是6.5%的出口换汇成本。将x=6.5代入回归模型中，得到 >> x=6.5; >> y=2.6597-0.2288*x y = 1.1725

matlab建立多元线性回归模型并进行显著性检验及预测问题

matlab建立多元线性回归模型并进行显着性检验及预测问题 例子; x=[143 145 146 147 149 150 153 154 155 156 157 158 159 160 162 164]'; X=[ones(16,1) x]; 增加一个常数项Y=[88 85 88 91 92 93 93 95 96 98 97 96 98 99 100 102]'; [b,bint,r,rint,stats]=regress(Y,X) 得结果：b = bint = stats = 即对应于b的置信区间分别为[，]、[,]; r2=, F=, p= p<, 可知回归模型y=+ 成立. 这个是一元的，如果是多元就增加X的行数！ function [beta_hat,Y_hat,stats]=regress(X,Y,alpha) % 多元线性回归(Y=Xβ+ε)MATLAB代码 %？ % 参数说明 % X：自变量矩阵，列为自变量，行为观测值 % Y：应变量矩阵，同X % alpha：置信度，[0 1]之间的任意数据 % beta_hat：回归系数 % Y_beata：回归目标值，使用Y-Y_hat来观测回归效果 % stats：结构体，具有如下字段 % =[fV,fH]，F检验相关参数，检验线性回归方程是否显着 % fV：F分布值，越大越好，线性回归方程越显着 % fH：0或1，0不显着；1显着(好) % =[tH,tV,tW]，T检验相关参数和区间估计，检验回归系数β是否与Y有显着线性关系 % tV：T分布值，beta_hat(i)绝对值越大,表示Xi对Y显着的线性作用% tH：0或1，0不显着；1显着 % tW：区间估计拒绝域，如果beta(i)在对应拒绝区间内，那么否认Xi对Y显着的线性作用 % =[T,U,Q,R]，回归中使用的重要参数 % T：总离差平方和，且满足T=Q+U % U：回归离差平方和 % Q：残差平方和 % R∈[0 1]：复相关系数，表征回归离差占总离差的百分比，越大越好% 举例说明 % 比如要拟合y=a+b*log(x1)+c*exp(x2)+d*x1*x2，注意一定要将原来方程线化% x1=rand(10,1)*10; % x2=rand(10,1)*10; % Y=5+8*log(x1)+*exp(x2)+*x1.*x2+rand(10,1); % 以上随即生成一组测试数据 % X=[ones(10,1) log(x1) exp(x2) x1.*x2]; % 将原来的方表达式化成Y=Xβ，注意最前面的1不要丢了

Matlab多元线性回归

Matlab多元线性回归 [ b , bint , r , rint , stats ]=regress ( y , x ) ， 其中b 是回归方程中的参数估计值，bint 是b 的置信区间，r 和rint 分别表示残差及残差对应的置信区间。StatS 数组包含三个数字，分别是相关系数，F 统计量及对应的概率p 值。 拟合结果： Y=b(1)x(1)+b(2)x(2)+b(3)x(3)+…+b(n)x(n) b(1)是系数,x(1)为全1的一个列向量。 注意：不是插值。 x=[1097 1284 1502 1394 1303 1555 1917 2051 2111 2286 2311 2003 2435 2625 2948 3155 3372];%因变量时间序列数据 y=[698 872 988 807 738 1025 1316 1539 1561 1765 1762 1960 1902 2013 2446 2736 2825];%自变量时间序列数据 X=[ones(size(x')),x']; [b,bint,r,rint,stats]=regress(y',X,0.05);%调用一元回归分析函数 rcoplot(r,rint)%画出在置信度区间下误差分布. 举例：

x = 1 2 4 9 1 4 3 7 1 5 9 0 1 9 1 8 >> y=[10 3 90 48]'; >> [ b , bint , r , rint , stats ]=regress ( y , x ) 得到的结果 b = -186.8333 16.0238 21.8571 8.5952 bint = NaN NaN NaN NaN

(实验三)连续时间LTI系统的频域分析汇总

实验三 连续时间LTI 系统的频域分析 一、实验目的 1、掌握系统频率响应特性的概念及其物理意义； 2、掌握系统频率响应特性的计算方法和特性曲线的绘制方法，理解具有不同频率响应特性的滤波器对信号的滤波作用； 3、学习和掌握幅度特性、相位特性以及群延时的物理意义； 4、掌握用MA TLAB 语言进行系统频响特性分析的方法。 基本要求：掌握LTI 连续和离散时间系统的频域数学模型和频域数学模型的MATLAB 描述方法，深刻理解LTI 系统的频率响应特性的物理意义，理解滤波和滤波器的概念，掌握利用MATLAB 计算和绘制LTI 系统频率响应特性曲线中的编程。 二、实验原理及方法 1 连续时间LTI 系统的频率响应 所谓频率特性，也称为频率响应特性，简称频率响应（Frequency response ），是指系统在正弦信号激励下的稳态响应随频率变化的情况，包括响应的幅度随频率的变化情况和响应的相位随频率的变化情况两个方面。 上图中x(t)、y(t)分别为系统的时域激励信号和响应信号，h(t)是系统的单位冲激响应，它们三者之间的关系为：)(*)()(t h t x t y =，由傅里叶变换的时域卷积定理可得到： )()()(ωωωj H j X j Y = 3.1 或者： ) () ()(ωωωj X j Y j H = 3.2 )(ωj H 为系统的频域数学模型，它实际上就是系统的单位冲激响应h(t)的傅里叶变换。即 ? ∞ ∞ --= dt e t h j H t j ωω)()( 3.3 由于H(j ω)实际上是系统单位冲激响应h(t)的傅里叶变换，如果h(t)是收敛的，或者说 是绝对可积（Absolutly integrabel ）的话，那么H(j ω)一定存在，而且H(j ω)通常是复数，

线性控制系统的频率响应分析

一．实验目的 1．了解和掌握对数幅频曲线和相频曲线（波德图）、幅相曲线（奈奎斯特图）的构造及绘制方法。 2．二阶开环系统中的相位裕度和幅值穿越频率的计算。 二．实验内容及要求 1．一阶惯性环节的频率特性曲线测试。 2．二阶开环系统的频率特性测试，研究表征系统稳定程度的相位裕度和 幅值穿越频率对系统的影响。 三、实验主要仪器设备和材料 1．labACT自控/计控原理实验机一台 2．数字存储示波器一台 四、实验方法、步骤及结果测试 1．一阶惯性环节的频率特性曲线 惯性环节的频率特性测试模拟电路见图4-1。 图4-1 惯性环节的频率特性测试模拟电路 实验步骤：注：‘S ST'不能用“短路套”短接！ （1）将数/模转换器（B2）输出OUT2作为被测系统的输入。 （2）按图4-1安置短路套及测孔联线。 （3）运行、观察、记录： ①运行LABACT程序，选择自动控制菜单下的线性控制系统的频率响应分析-实验项目，选择一阶系统，再选择开始实验，点击开始，实验机将自动产生0.5Hz~64Hz多个频率信号，测试被测系统的频率特性，等待将近十分钟，测试结束。 ②测试结束后，可点击界面下方的“频率特性”选择框中的任意一项进行切换，将显示被测系统的对数幅频、相频特性曲线（伯德图）和幅相曲线（奈 奎斯特图），同时在界面上方将显示点取的频率点的L、、Im、Re等相关数

据。如点击停止，将停止示波器运行，不能再测量数据。 ③分别改变惯性环节开环增益与时间常数，观察被测系统的开环对数幅频曲线、相频曲线及幅相曲线，在幅频曲线或相频曲线上点取相同的频率点，测量、记录数据于实验数据表中。 实验数据表1：改变惯性环节开环增益，（T=0.05，C=1u，R2=50K） 实验数据表2： 改变惯性环节时间常数， K=1（R1=50K、R2=50K） 2．二阶开环系统的频率特性曲线 二阶系统模拟电路图的构成如图4-2所示。

matlab建立多元线性回归模型并进行显著性检验及预测问题

matlab建立多元线性回归模型并进行显著性检 验及预测问题 例子; x=[143 145 146 147 149 150 153 154 155 156 157 158 159 160 162 164]'; X=[ones(16,1) x]; 增加一个常数项 Y=[88 85 88 91 92 93 93 95 96 98 97 96 98 99 100 102]'; [b,bint,r,rint,stats]=regress(Y,X) 得结果：b = bint = stats = 即对应于b的置信区间分别为[，]、[,]; r2=, F=, p= p<, 可知回 归模型 y=+ 成立. 这个是一元的，如果是多元就增加X的行数！ function [beta_hat,Y_hat,stats]=regress(X,Y,alpha) % 多元线性回归(Y=Xβ+ε)MATLAB代码 % % 参数说明 % X：自变量矩阵，列为自变量，行为观测值 % Y：应变量矩阵，同X % alpha：置信度，[0 1]之间的任意数据 % beta_hat：回归系数 % Y_beata：回归目标值，使用Y-Y_hat来观测回归效果 % stats：结构体，具有如下字段 % =[fV,fH]，F检验相关参数，检验线性回归方程是否显著 % fV：F分布值，越大越好，线性回归方程 越显著 % fH：0或1，0不显著；1显著(好) % =[tH,tV,tW]，T检验相关参数和区间估计，检验回归系数β是 否与Y有显著线性关系 % tV：T分布值，beta_hat(i)绝对值越大, 表示Xi对Y显著的线性作用 % tH：0或1，0不显著；1显著 % tW：区间估计拒绝域，如果beta(i)在对 应拒绝区间内，那么否认Xi对Y显著的线性作用 % =[T,U,Q,R]，回归中使用的重要参数 % T：总离差平方和，且满足T=Q+U % U：回归离差平方和 % Q：残差平方和 % R∈[0 1]：复相关系数，表征回归离差占总 离差的百分比，越大越好 % 举例说明 % 比如要拟合 y=a+b*log(x1)+c*exp(x2)+d*x1*x2，注意一定要将原来方程 线化 % x1=rand(10,1)*10;

(完整版)第十一章电路的频率响应

第十一章 电路的频率响应 11-1 网络函数 11-2 RLC 串联电路的谐振 11-3 RLC 串联电路的频率响应 11-4 RLC 并联谐振电路 11-5 波特图 11-6 滤波器简介 重点 1. 网络函数 2. 串、并联谐振的概念 11-1 网络函数 当电路中激励源的频率变化时，电路中的感抗、容抗将跟随频率变化，从而导致电路的工作状态亦跟随频率变化。因此，分析研究电路和系统的频率特性就显得格外重要。 频率特性 电路和系统的工作状态跟随频率而变化的现象，称为电路和系统的频率特性，又称频率响应。 1. 网络函数H (j ω)的定义 在线性正弦稳态网络中，当只有一个独立激励源作用时，网络中某一处的响应（电压或电流）与网络输入之比，称为该响应的网络函数。 def (j )(j )(j ) R H E ωωω=

2. 网络函数H(j ω)的物理意义 ⑴ 驱动点函数 激励是电流源，响应是电压 策动点阻抗 激励是电压源，响应是电流 策动点导纳 ⑵ 转移函数(传递函数) 激励是电压源 转移导纳 转移电压比 (j ) I ω(j U 1(U 1(j )I ω(j )(j )(j ) U H I ωωω= (j )(j )(j ) I H U ωωω= 21(j )(j )(j )I H U ωωω= 21(j ) (j )(j ) U H U ωωω=

激励是电流源 转移阻抗 转移电流比 注意 ①H(j ω)与网络的结构、参数值有关，与输入、输出变量的类型以及端口对的相互位置有关，与输入、输出幅值无关。因此网络函数是网络性质的一种体现。 ②H(j ω) 是一个复数，它的频率特性分为两个部分： 幅频特性 ：模与频率的关系 ()H j ωω - 相频特性：幅角与频率的关系 ()j ?ωω - ③网络函数可以用相量法中任一分析求解方法获得。 例1-1 求图示电路的网络函数 2 S I U ? ? 和 L S U U ? ? 解：列网孔方程解电流 _ 2 I 1 I 21(j ) (j )(j ) U H I ωωω= 21(j ) (j )(j ) I H I ωωω= 12s 12(2j )22(4j )0 I I U I I ωω?+-=??-++=??s 2224(j )j6U I ωω = ++

多元回归分析matlab剖析

回归分析MATLAB 工具箱 一、多元线性回归 多元线性回归：p p x x y βββ+++=...110 1、确定回归系数的点估计值： 命令为：b=regress(Y , X ) ①b 表示???? ?? ????????=p b βββ?...??10 ②Y 表示????????????=n Y Y Y Y (2) 1 ③X 表示??? ??? ????? ???=np n n p p x x x x x x x x x X ...1......... .........1 (12) 1 22221 11211 2、求回归系数的点估计和区间估计、并检验回归模型： 命令为：[b, bint,r,rint,stats]=regress(Y ,X,alpha) ①bint 表示回归系数的区间估计. ②r 表示残差. ③rint 表示置信区间. ④stats 表示用于检验回归模型的统计量,有三个数值：相关系数r 2、F 值、与F 对应的概率p. 说明：相关系数2 r 越接近1,说明回归方程越显著；)1,(1-->-k n k F F α时拒绝0H ,F 越大,说明回归方程越显著；与F 对应的概率p α<时拒绝H 0,回归模型成立. ⑤alpha 表示显著性水平(缺省时为0.05) 3、画出残差及其置信区间. 命令为：rcoplot(r,rint) 例1.如下程序. 解：(1)输入数据. x=[143 145 146 147 149 150 153 154 155 156 157 158 159 160 162 164]'; X=[ones(16,1) x]; Y=[88 85 88 91 92 93 93 95 96 98 97 96 98 99 100 102]'; (2)回归分析及检验. [b,bint,r,rint,stats]=regress(Y ,X) b,bint,stats 得结果：b = bint =

基于MATLAB的线性电路频率响应特性分析 上传

课程设计任务书 学生姓名： 专业班级： 指导教师： 工作单位： 题 目: 基于MATLAB 的线性电路频率响应特性分析 初始条件： MATLAB 软件，微机 主要任务： 利用MATLAB 强大的图形处理功能、符号运算功能和数值计算功能，实现线性电路频率响应特性的仿真波形。 1）绘出RLC 串联电路中AuR=UR/U1的幅频特性及相频特性曲线； 2）绘出RLC 串联电路中AuC=UC/U1的幅频特性及相频特性曲线； 3）绘出RLC 串联电路中AuL=UL/U1的幅频特性及相频特性曲线； 4）设L=1H ，C=1F ， ,/10 s rad =ω 改变R 之值，观察各特性曲线的变化情况； 5）撰写MATLAB 课程设计说明书 时间安排： 学习MATLAB 语言的概况 第1天 学习MATLAB 语言的基本知识 第2、3天 学习MATLAB 语言的应用环境，调试命令，绘图能力 第4、5天 课程设计 第6-9天 答辩 第10天 指导教师签名： 年 月 日 系主任（或责任教师）签名： 年 月 日

摘要 MATLAB语言具备高效、可视化及推理能力强等特点，是目前工程界流行最广的科学计算语言。特别是在电子通信领域，MATLAB常常被用于进行电路、信号与系统、数字信号处理等多个方面的理论验证与演算求解。将MATLAB软件引入到电路分析中,大大地提高了计算精度和工作效率,为电路分析提供了一个有效的辅助工具，是电子工程人员不可或缺的辅助工具软件。 本次课程设计基于MATLAB强大的图形处理功能、符号运算功能和数值计算，着重对于线性电路中较有代表的RLC串联电路的频率响应进行分析，着重训练MATLAB在电路分析的应用，能够运用相关软件进行数学模型建立、相关参量求解、结果呈现与分析。从而达到对MATLAB软件及其程序编写方式的熟悉。 关键字：MATLAB 线性电路频率响应 Abstract MATLAB language with high efficiency, visualization and reasoning ability and other characteristics, is the current practice of the most widely popular scientific computing language. Especially in the field of electronic communications, MATLAB is often used for circuits, signals and systems, digital signal processing and other aspects of the theory of authentication and routing solution. MATLAB software is introduced into the circuit analysis, greatly improves the accuracy and efficiency. It is an effective auxiliary circuit analysis tools. MATLAB is an indispensable auxiliary tool for electronic engineers. This course design based on MATLAB powerful graphics capabilities, and numerical computation symbolic operation, focuses on the frequency response of RLC series circuit which represented the linear circuit analysis. Training in the application of MATLAB in circuit analysis, make us be able to use relevant software to mathematical modeling, solve the relevant parameters, present and analyze the results. After the design, we will be able to achieve the MATLAB software and its programming on the way to the familiar. Keywords: MATLAB frequency response of linear circuits

Matlab多变量回归分析教程

本次教程的主要内容包含： 一、多元线性回归 2# 多元线性回归：regress 二、多项式回归 3# 一元多项式：polyfit或者polytool 多元二项式：rstool或者rsmdemo 三、非线性回归 4# 非线性回归：nlinfit 四、逐步回归 5# 逐步回归：stepwise 一、多元线性回归 多元线性回归： 1、b=regress(Y, X ) 确定回归系数的点估计值

2、[b, bint,r,rint,stats]=regress(Y,X,alpha)求回归系数的点估计和区间估计、并检验回归模型 ①bint表示回归系数的区间估计. ②r表示残差 ③rint表示置信区间 ④stats表示用于检验回归模型的统计量,有三个数值：相关系数r2、F值、与F对应的概率p 说明：相关系数r2越接近1，说明回归方程越显著；时拒绝H0，F越大，说明回归方程越显著；与F对应的概率p<α时拒绝H0 ⑤alpha表示显著性水平(缺省时为0.05) 3、rcoplot(r,rint)画出残差及其置信区间 具体参见下面的实例演示 4、实例演示，函数使用说明 (1)输入数据 1.>>x=[143 145 146 147 149 150 153 154 155 156 157 158 159 160 162 164]'; 2.>>X=[ones(16,1) x]; 3.>>Y=[88 85 88 91 92 93 93 95 96 98 97 96 98 99 100 102]'; 复制代码 (2)回归分析及检验 1. >> [b,bint,r,rint,stats]=regress(Y,X) 2. 3. b = 4. 5. -1 6.0730 6.0.7194 7. 8. 9.bint =

matlab多元非线性回归教程

matlab 回归（多元拟合）教程 前言 1、学三条命令 polyfit(x,y,n)---拟合成一元幂函数（一元多次） regress(y,x)----可以多元， nlinfit(x,y,’fun ’,beta0) (可用于任何类型的函数，任意多元函数，应用范围最主，最万能的) 2、同一个问题，这三条命令都可以使用，但结果肯定是不同的，因为拟合的近似结果，没有唯一的标准的答案。相当于咨询多个专家。 3、回归的操作步骤： 根据图形（实际点），选配一条恰当的函数形式（类型）---需要数学理论与基础和经验。（并写出该函数表达式的一般形式，含待定系数）------选用某条回归命令求出所有的待定系数。所以可以说，回归就是求待定系数的过程（需确定函数的形式） 一、回归命令 一元多次拟合polyfit(x,y,n);一元回归polyfit;多元回归regress---nlinfit(非线性) 二、多元回归分析 对于多元线性回归模型(其实可以是非线性，它通用性极高)： e x x y p p ++ ++ = βββ 1 10 设变量12,,,p x x x y 的n 组观测值为12(,, ,)1,2, ,i i ip i x x x y i n = 记 ??????? ??=np n n p p x x x x x x x x x x 2 1 22221 1121111 1，?? ?? ? ?? ??=n y y y y 21，则?????? ? ??=p ββββ 10 的估计值为排列方式与线性代数中的线性方程组相同（），拟合成多元函数---regress 使用格式：左边用b=[b, bint, r, rint, stats]右边用=regress(y, x)或regress(y, x, alpha) ---命令中是先y 后x, ---须构造好矩阵x(x 中的每列与目标函数的一项对应) ---并且x 要在最前面额外添加全1列/对应于常数项

[频响] 频响分析方法总结

频响分析，或者叫稳态动力学分析在abaqus中包括以下三种方法： 直接稳态动力学分析（direct solution steady state dynamic analysis） 模态稳态动力学分析（mode based steady state dynamic analysis） 子空间稳态动力学分析（subspace projection steady state dynamic analysis） 1）直接稳态动力学 优点：在直接稳态动力学分析中，系统的稳态谐波响应是通过对模型的原始方程直接积分计算出来的。如果分析的对象存在非对称刚度、包含模态阻尼以外的其他阻尼或者必须考虑粘弹性材料特性（频变特性），则不能提取特征模态的情况下，可以应用直接法进行稳态响应的计算和分析。 缺点：进行直接稳态动力学分析不需要提取系统的特征模态，而是在每个频率点对整个模型进行复杂的积分运算。因此，对于具有大阻尼和频变特性的模型，应用直接法比模态分析方法精确，但是耗时较多。 2）模态稳态动力学分析 模态稳态动力学分析方法是基于模态叠加法求解系统的稳态响应。因此，在求解稳态响应之前必须先提取无阻尼系统的特征模态，也就是在说必须在step steady state dynamics,modal 前加一步step frequency。另外，必须确定需要保留的特征模态，以确保能够精确描述系统的动力学特性，也就是说如果是进行0-1000hz的分析，step frequency的number of eigenvalues requested选定的阶数的模态频率必须大于1000hz，简单的作法是这里选all……，下面的maximum……填入1000。 模态稳态动力学分析的特点：相较于直接法和子空间法分析速度快，耗时最少，计算精度低于直接法和子空间法，不适合于分析具有大阻尼特性的模型，不适合于分析具有粘弹性材料（频变特性）的模型。 3）子空间稳态动力学分析 子空间稳态动力学分析的基本思想是：首先提取无阻尼、对称系统的特征模态，并选取适当的特征向量组成特征模态子空间，然后将稳态动力学方程组投影到特征模态子空间上，通过直接法求解子空间的稳态动力学方程。 我的感觉是子空间法是直接法和模态法的折中，它的特点是模型可以定义任意形式的阻尼，可以处理具有非对称刚度矩阵的模型，可以处理具有频变特性的模型，计算时间和精度也是在直接法和模态法的中间。 直接法在定义边界条件时通过选项*boundary的amplitude参数来引用频变幅值，但这里默认的好像是位移，如果我有的是加速度或者速度数据，想用直接法进行分析应该如何设定呢，希望知道的大神能相告。 模态法和子空间法不能使用*boundary选项定义边界条件的运动，而只能通过选项*base motion来定义边界条件的运动。

matlab多元非线性回归及显着性分析(实例)

matlab多元非线性回归及显著性分析 给各位高手：小弟有一些数据需要回归分析（非线性）及显著性检验（回归模型，次要项，误差及失拟项纯误差，F值和P值），求大侠帮助，给出程序，不胜感激。 模型： DA TA=... %DA TA前三列是影响因子，第四列为响应值 [2 130 75 48.61; 2 110 75 56.43; 2 130 45 61.32; 2 110 45 65.28; 1 110 45 55.80; 1 130 75 45.65; 1 110 75 50.91; 1 130 45 67.94; 1.5 120 60 74.15; 1.5 120 60 71.28; 1.5 120 60 77.95; 1.5 120 60 74.16; 1.5 120 60 75.20; 1.5 120 85 35.65; 1.5 140 60 48.66; 1.5 120 30 74.10; 1.5 100 60 6 2.30; 0.5 120 60 66.00; 2.5 120 60 75.10]; 回归分析过程： （1）MATLAB编程步骤1：首先为非线性回归函数编程，程序存盘为user_function.m function y=user_function(beta,x) b0 = beta(1); b1 = beta(2); b2 = beta(3); b3 = beta(4); x0 = x(:,1); x1 = x(:,2); x2 = x(:,3); x3 = x(:,4); y=b0*x0+b1*x1.^2+b2*x2.^2+b3*x3.^2; （2）MATLAB编程步骤2：编写非线性回归主程序，程序运行时调用函数user_function x=[1 2 130 75 48.61;

频响分析

radioss频响分析 材料 属性T=*** 1.loadcollector spc DOF 123456 钻柱 井壁 2. A DOF1=2.54mm DAREA DAREA 载荷激励 SPCD 强制位移、速度激励、加速度激励 如果是SPCD，则激励处还需添加相应自由度的SPC约束 3. B card image=TABLED1 x(1)=0,y(1)=1,x(2)=1000,y(2)=1 如果是激励曲线，则从utility——>table creat中导入 4. OMEGA card image=FREQi 勾选FREQ1,F1=20,DF=20,NDF=49 5. RLOAD2card image=RLOAD2 EXCITED——>A TB——>B TP——>φ DELAY——>τ DPHASE——>θ TYPE=LOAD 如果有好几个载荷，则用DLOAD组合 6.loadstep type=freq.resp(direct) SPC——>SPC DLOAD——>RLOAD2 FREQ——>OMEGA 7.定义set type=SET_GRID 6.control cards Displacements format=HG, DISP_FORM=PHASE, DISP_OPT=SID PARAM coupmass:yes G=0.06 OUTPUT keyword=HGFREQ FREQ=ALL 6.loadstep type=freq.resp(direct) 6.loadstep type=freq.resp(direct) 直接频响 模态频响还需设置EIGRL卡片 汽车白车身 输入：白车身与底盘相连的点 输出：方向盘，底板、座椅…

(研究生-数理统计)多元线性回归及显著性检验Matlab程序(完美版)

多元线性回归及显著性检验Matlab程序（完美版） 一、说明： 1、本程序是研究生教材《数理统计》（杨虎、刘琼、钟波编著）例4.4.1（P133）的Matlab 编程解答程序。教材上的例题只做了回归方程显著性分析和一次回归系数显著性分析（剔除x1后没有再检验x2和x3）。 2、本程序在以上的基础之上,还分别检验了x2和x3,并且计算精度更高。 3、本程序可根据用户的需要,在输入不同的显著性水平α之下得到相应的解答。 4、本程序移植性强,对于其他数据,只需要改变excel中的数据即可。 5、本程序输出的可读性强,整洁美观。 二、数据入下（将数据存入excel表格,文件名为jc_p133_example.xls。注意数据是按 ）：

三、完整程序如下： %----------------------------by ggihhimm---------------------------- %《数理统计》杨虎、刘琼、钟波编著例4.4.1 多元线性回归及显著性检验完整解答 % 输入需要的显著水平α（默认α=0.02）,计算出不同结果（见运行结果） % 该程序也适合其他维数的数据分析（只需改变excel表格中的数据即可） %----------------------------by ggihhimm---------------------------- clear;clc; data=xlsread('jc_p133_example.xls','sheet1'); xi=data(:,1:end-1); [n,k]=size(data); k=k-1; index_of_xi_array=ones(1,k); X=[ones(n,1) xi]; Y=data(:,end); fprintf('第1次计算结果：\r') beta_mao=((X'*X)\X'*Y)'; fmt_str0=''; for i0=1:k+1 fmt_str0=[fmt_str0 'β' num2str(i0-1) ' = %0.4f\r']; end fprintf(fmt_str0,beta_mao) fprintf('\r')