南昌大学自动控制原理实验报告实验一

合集下载

自动控制原理实验报告(一、二阶系统的电子模拟及时域响应的动态测试等三个实验)

自动控制原理实验报告(一、二阶系统的电子模拟及时域响应的动态测试等三个实验)

自动控制原理实验报告作者姓名学科专业机械工程及自动化班级学号X X年10月27日实验一一、二阶系统的电子模拟及时域响应的动态测试一、实验目的1、了解一、二阶系统阶跃响应及其性能指标与系统参数之间的关系。

2、学习在电子模拟机上建立典型环节系统模型的方法。

3、学习阶跃响应的测试方法。

二、实验内容1、建立一阶系统的电子模型,观测并记录在不同时间常数T时的跃响应曲线,并测定其过渡过程时间TS。

2、建立二阶系统的电子模型,观测并记录在不同阻尼比ζ时的跃响应曲线,并测定其超调量σ%及过渡过程时间TS。

三、实验原理1、一阶系统阶跃响应性能指标的测试系统的传递函数为:()s()1C s KR s Ts φ=+()=模拟运算电路如下图:其中21R K R =,2T R C =;在实验中,始终保持21,R R =即1K =,通过调节2R 和C 的不同取值,使得T 的值分别为0.25,0.5,1。

记录实验数据,测量过度过程的性能指标,其中按照经验公式取3s t T=2、二阶系统阶跃响应性能指标的测试系统传递函数为:令ωn=1弧度/秒,则系统结构如下图:二阶系统的模拟电路图如下:在实验过程中,取22321,1R C R C ==,则442312R R C R ζ==,即4212R C ζ=;在实验当中取123121,1R R R M C C F μ===Ω==,通过调整4R 取不同的值,使得ζ分别为0.25,0.5,1;记录所测得的实验数据以及其性能指标,其中经验公式为3.5%100%,s net σζω=⨯=.四、试验设备:1、HHMN-1型电子模拟机一台。

2、PC机一台。

3、数字万用表一块。

4、导线若干。

五、实验步骤:1、熟悉电子模拟机的使用,将各运算放大器接成比例器,通电调零。

2、断开电源,按照实验说明书上的条件和要求,计算电阻和电容的取值,按照模拟线路图搭接线路,不用的运算放大器接成比例器。

3、将D/A输出端与系统输入端Ui连接,将A/D1与系统输出端UO连接(此处连接必须谨慎,不可接错)。

实验报告-自动控制原理

实验报告-自动控制原理
________________________________________________________________________________
________________________________________________________________________________
〖分析பைடு நூலகம்:______________________________________________________________________
_______________________________________________________________________________
说明:特征参数为比例增益K和微分时间常数T。
1)R2=R1=100KΩ, C2=0.01µF,C1=1µF;特征参数实际值:K=______,T=________。
波形如下所示:
2)R2=R1=100KΩ, C2=0.01µF,C1=0.1µF;特征参数实际值:K= 1,T=0.01。
波形如下所示:
四、实验心得体会
实验报告
班级
姓名
学号
所属课程
《自动控制原理》
课时
2
实践环节
实验3控制系统的稳定性分析
地点
实字4#318
所需设备
电脑、工具箱
一、实验目的
1.观察系统的不稳定现象。
2.研究系统开环增益和时间常数对稳定性的影响
3.学习用MATLAB仿真软件对实验内容中的电路进行仿真。
2、实验步骤
_______________________________________________________________________________

自动控制原理实验报告

自动控制原理实验报告

《自动控制原理》实验报告学院:专业:班级:姓名:学号:实验报告实验一实验名称:典型环节响应实验姓名:实验组别:实验日期:年月日成绩:一、实验目的1、学习设计构成典型环节的模拟电路,掌握典型环节的特性以及电路参数对特性的影响。

2、学习典型环节响应的测量方法,对比实验结果与理论分析。

二、实验设备1、计算机一台,实验软件一套。

2、实验箱一套。

3、面包板、导线、电阻、电容、运算放大器等器件若干。

三、实验原理设计构成下列典型环节的模拟电路,测量其阶跃响应1、比例环节:利用运算放大器、电阻、电容等元件设计一模拟电路,输入输出之间满足如下运动方程:C(t)=KR(t) 传递函数G(s)=K2、惯性环节利用运算放大器、电阻、电容等元件设计一模拟电路,输入输出之间满足如下运动方程:传递函数G(s)=K/(Ts+1)3、积分环节利用运算放大器、电阻、电容等元件设计一模拟电路,输入输出之间满足如下运动方程:传递函数G(s)=K/s4、微分环节利用运算放大器、电阻、电容等元件设计一模拟电路,输入输出之间满足如下运动方程:传递函数G(s)=Ks四、实验结果及分析K=R1/R2=1/2比例环节模拟电路图比例环节阶跃响应惯性环节模拟电路图K=R2/R1=2T=R2*C=200K*1uf=0.2惯性环节阶跃响应积分环节模拟电路图T=RC=100K*1uf=0.1积分环节阶跃响应微分环节模拟电路微分环节阶跃响应微分环节阶跃响应仿真实验报告实验 二 实验名称: 二阶系统及其阶跃响应实验 姓名: 实验组别: 实验日期: 年 月 日 成绩:一、实验目的1、学习设计构成二阶系统的模拟电路,掌握二阶系统的特性以及电路参数对特性的影响。

2、学习二阶系统阶跃响应的测量方法,对比实验结果与理论分析。

二、实验设备1、 计算机一台,实验软件一套。

2、 实验箱一套。

3、 面包板、导线、电阻、电容、运算放大器等器件若干。

三、实验原理设计构成二阶振荡环节的模拟电路,测量其阶跃响应利用运算放大器、电阻、电容等元件设计一模拟电路,输入输出之间满足如下运动方程:)()()(2121222t R K K K K dt t dC dt t C d =++τ传递函数212221)(K K s s K K s ++=τψ四、实验结果及分析二阶系统模拟电路图二阶系统模拟电路结构图R1=200K R2=200K R=100K C取1uf和0.1uf两个值W=1/T=1/RC=1/RCK/2=R2/2*R1二阶系统阶跃响应仿真结果实验报告实验三实验名称:系统频率特性的测量姓名:实验组别:实验日期:年月日成绩:一、实验目的1、加深了解系统及元件频率特性的物理意义。

自动控制原理实验报告

自动控制原理实验报告

自动控制原理实验报告实验一典型环节的电路模拟与软件仿真研究一.实验目的1.通过实验熟悉并掌握实验装置和上位机软件的使用方法。

2.通过实验熟悉各种典型环节的传递函数及其特性,掌握电路模拟和软件仿真研究方法二.实验内容1.设计各种典型环节的模拟电路。

2.完成各种典型环节模拟电路的阶跃特性测试,并研究参数变化对典型环节阶跃特性的影响。

3.在上位机界面上,填入各个环节的实际(非理想)传递函数参数,完成典型环节阶跃特性的软件仿真研究,并与电路模拟研究的结果作比较。

三.实验步骤1.熟悉实验箱,利用实验箱上的模拟电路单元,参考本实验附录设计并连接各种典型环节(包括比例、积分、比例积分、比例微分、比例积分微分以及惯性环节)的模拟电路。

注意实验接线前必须先将实验箱上电,以对运放仔细调零。

然后断电,再接线。

接线时要注意不同环节、不同测试信号对运放锁零的要求。

在输入阶跃信号时,除比例环节运放可不锁零(G可接-15V)也可锁零外,其余环节都需要考虑运放锁零。

2.利用实验设备完成各典型环节模拟电路的阶跃特性测试,并研究参数变化对典型环节阶跃特性的影响。

无上位机时,利用实验箱上的信号源单元U2所输出的周期阶跃信号作为环节输入,即连接箱上U2的“阶跃”与环节的输入端(例如对比例环节即图1.1.2的Ui),同时连接U2的“锁零(G)”与运放的锁零G。

然后用示波器观测该环节的输入与输出(例如对比例环节即测试图1.1.2的Ui和Uo)。

注意调节U2的周期阶跃信号的“频率”电位器RP5与“幅值”电位器RP2,以保证观测到完整的阶跃响应过程。

有上位机时,必须在熟悉上位机界面操作的基础上,充分利用上位机提供的虚拟示波器与信号发生器功能。

为了利用上位机提供的虚拟示波器与信号发生器功能,接线方式将不同于上述无上位机情况。

仍以比例环节为例,此时将Ui连到实验箱 U3单元的O1(D/A通道的输出端),将Uo连到实验箱 U3单元的I1(A/D 通道的输入端),将运放的锁零G连到实验箱 U3单元的G1(与O1同步),并连好U3单元至上位机的USB2.0通信线。

自动控制原理实验报告

自动控制原理实验报告

自动控制原理实验报告实验一、典型环节的时域响应一.实验目的1. 熟悉并掌握TD-ACC+( TD-ACS设备的使用方法及各典型环节模拟控制电路的构成方法。

2. 熟悉各种典型环节的理想阶跃曲线和实际阶跃响应曲线。

对比差异、分析原因。

3. 了解参数变化对典型环节动态特性的影响。

二.实验设备PC机一台,TD-ACC+( TD-ACS实验系统一套。

三.实验内容1. 比例环节2. 积分环节3. 比例积分环节4. 惯性环节5. 比例微分环节6. 比例积分微分环节四、实验感想在本次实验后,我了解了典型环节的时域响应方面的知识,并且通过实践,实现了时域响应相关的操作,感受到了实验成功的喜悦。

实验二、线性系统的矫正一、目的要求1.掌握系统校正的方法,重点了解串联校正。

2.根据期望的时域性能指标推导出二阶系统的串联校正环节的传递函数、仪器设备PC机一台,TD-ACC+或TD-ACS)教学实验系统一套三、原理简述所谓校正就是指在使系统特性发生变接方式可分为馈回路之内采用的测点之后和放1.原系统的结构框图及性能指标对应的模拟电路图2.期望校正后系统的性能指标3.串联校正环节的理论推导四、实验现象分析校正前:校正后:校正前:校正后:六、实验心得次实验让我进一步熟悉了TD-ACC实验系统的使用,进一步学习了虚拟仪器,更加深入地学习了自动控制原理,更加牢固地掌握了相关理论知识,激发了我理论学习的兴趣。

实验三、线性系统的频率响应分析、实验目的1 .掌握波特图的绘制方法及由波特图来确定系统开环传函2 .掌握实验方法测量系统的波特图。

、实验设备PC机一台,TD-ACC系列教学实验系统一套三、实验原理及内容(一)实验原理1 .频率特性当输入正弦信号时,线性系统的稳态响应具有随频率(3由0变至%)而变化的特性。

频率响应法的基本思想是:尽管控制系统的输入信号不是正弦函数,而是其它形式的周期函数或非周期函数,但是,实际上的周期信号,都能满足狄利克莱条件,可以用富氏级数展开为各种谐波分量;而非周期信号也可以使用富氏积分表示为连续的频谱函数。

自动控制原理实验报告 (1)

自动控制原理实验报告 (1)

实验1 控制系统典型环节的模拟实验(一)实验目的:1.掌握控制系统中各典型环节的电路模拟及其参数的测定方法。

2.测量典型环节的阶跃响应曲线,了解参数变化对环节输出性能的影响。

实验原理:控制系统模拟实验采用复合网络法来模拟各种典型环节,即利用运算放大器不同的输入网络和反馈网络模拟各种典型环节,然后按照给定系统的结构图将这些模拟环节连接起来,便得到了相应的模拟系统。

再将输入信号加到模拟系统的输入端,并利用计算机等测量仪器,测量系统的输出,便可得到系统的动态响应曲线及性能指标。

实验内容及步骤实验内容:观测比例、惯性和积分环节的阶跃响应曲线。

实验步骤:分别按比例,惯性和积分实验电路原理图连线,完成相关参数设置,运行。

①按各典型环节的模拟电路图将线接好(先接比例)。

(PID先不接)②将模拟电路输入端(U i)与阶跃信号的输出端Y相连接;模拟电路的输出端(Uo)接至示波器。

③按下按钮(或松开按钮)SP时,用示波器观测输出端的实际响应曲线Uo(t),且将结果记下。

改变比例参数,重新观测结果。

④同理得积分和惯性环节的实际响应曲线,它们的理想曲线和实际响应曲线。

实验数据实验二控制系统典型环节的模拟实验(二)实验目的1.掌握控制系统中各典型环节的电路模拟及其参数的测定方法。

2.测量典型环节的阶跃响应曲线,了解参数变化对环节输出性能的影响。

实验仪器1.自动控制系统实验箱一台2.计算机一台实验原理控制系统模拟实验采用复合网络法来模拟各种典型环节,即利用运算放大器不同的输入网络和反馈网络模拟各种典型环节,然后按照给定系统的结构图将这些模拟环节连接起来,便得到了相应的模拟系统。

再将输入信号加到模拟系统的输入端,并利用计算机等测量仪器,测量系统的输出,便可得到系统的动态响应曲线及性能指标。

实验内容及步骤内容:观测PI,PD和PID环节的阶跃响应曲线。

步骤:分别按PI,PD和PID实验电路原理图连线,完成相关参数设置,运行①按各典型环节的模拟电路图将线接好。

自动控制原理实验报告

自动控制原理实验报告

实验一典型环节的模拟研究及阶跃响应分析1、比例环节可知比例环节的传递函数为一个常数:当Kp 分别为,1,2时,输入幅值为的正向阶跃信号,理论上依次输出幅值为,,的反向阶跃信号。

实验中,输出信号依次为幅值为,,的反向阶跃信号,相对误差分别为1.8%,2.2%,0.2%.在误差允许范围内可认为实际输出满足理论值。

2、 积分环节积分环节传递函数为:〔1〕T=0.1(0.033)时,C=1μf(0.33μf),利用MATLAB ,模拟阶跃信号输入下的输出信号如图:与实验测得波形比较可知,实际与理论值较为吻合,理论上时的波形斜率近似为时的三倍,实际上为,在误差允许范围内可认为满足理论条件。

3、 惯性环节惯性环节传递函数为:K = R f /R 1,T = R f C,(1) 保持K = R f /R 1= 1不变,观测秒,秒〔既R 1 = 100K,C = 1μf ,μf 〕时的输出波形。

利用matlab 仿真得到理论波形如下:时t s 〔5%〕理论值为300ms,实际测得t s =400ms 相对误差为:〔400-300〕/300=33.3%,读数误差较大。

K 理论值为1,实验值,相对误差为〔〕/2.28=7%与理论值较为接近。

时t s 〔5%〕理论值为30ms,实际测得t s =40ms 相对误差为:〔40-30〕/30=33.3% 由于ts 较小,所以读数时误差较大。

K 理论值为1,实验值,相对误差为〔〕/2.28=7%与理论值较为接近(2) 保持T = R f s 不变,分别观测K = 1,2时的输出波形。

K=1时波形即为〔1〕中时波形K=2时,利用matlab 仿真得到如下结果:t s 〔5%〕理论值为300ms,实际测得t s =400ms相对误差为:〔400-300〕/300=33.3% 读数误差较大K 理论值为2,实验值, 相对误差为〔〕/2=5.7%if i o R RU U -=1TS K)s (R )s (C +-=与理论值较为接近。

自动控制原理实验报告一

自动控制原理实验报告一

【范例 2-3】
【解】
将系统 G(s)=
转化为部分分式展开式 ������������ +������������������ +������
������������ +������������+������
MATLAB 程序为: >> num=[1,5,6];den=[1,2,1,0] >> [r,p,k]=residue(num,den) r = -5 -2 6 p = -1 -1 0 k = []
������������(������+������) ������+������.������ ������+������ (������+������)
【范例 2-2】
已知系统传递函数 G(s)= 的零极点模型 【解】MATLAB 程序为: >> k=10 >> z=[-5] >> p=[-0.5,-2,-3] >> sys=zpk(z,p,k) Zero/pole/gain:
试建立控制系统
3
控制系统的数学模型
10 (s+5) ------------------(s+0.5) (s+2) (s+3)
【自我实践 2-2】
建立控制系统的零极点模型 1) 2) G(s)= G(s)=
������ ������+������−������ (������+������+������) ������������ ������+������ ������+������ (������������ +������) ������ ������ ������+������ (������������ +������������ +������)
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

南昌大学实验报告
学生姓名:王瑾然学号:6101113031 专业班级:电气131班
实验类型:■ 验证□ 综合□ 设计□ 创新实验日期:实验成绩:
一、实验项目名称
实验3.1.1 典型环节的模拟研究
二、实验要求
1.了解和掌握各典型环节模拟电路的构成方法、传递函数表达式及输出时域函
数表达式。

2.观察和分析各典型环节的阶跃响应曲线,了解各项电路参数对典型环节动态
特性的影响。

三、主要仪器设备及耗材
1.计算机一台(Wind ows XP操作系统)
2.AEDK-labACT自动控制理论教学实验系统一套
bACT6_08软件一套
四、实验内容和步骤
1.观察比例环节的阶跃响应曲线
(1)打开虚拟示波器的界面,点击开始,按下信号发生器的阶跃信号按钮(0→+4V阶跃),用示波器观测A6输出端(Uo)的实际响应曲线Uo(t)。

(2)改变比例系数,重新观测结果,填入实验报告。

2.观察惯性环节的阶跃响应曲线
(1)打开虚拟示波器的界面,点击开始,用示波器观测A6输出端,按下信号发生器的阶跃信号按钮时,等待完整波形出来后,移动虚拟示波器横游标到4V×0.632处,得到与惯性的曲线的交点,再移动虚拟示波器两根纵游标,从阶跃开始到曲线的交点,量得惯性环节模拟电路时间常数T。

A6输出端的实际响应曲线Uo(t)。

(2)改变时间常数及比例系数,重新观测结果,填入实验报告。

3.观察积分环节的阶跃响应曲线
(1)打开虚拟示波器的界面,点击开始,用示波器观测A6输出端(Uo),调节调宽电位器使宽度从0.3秒开始调到积分输出在虚拟示波器顶端为止。

(2)等待完整波形出来后,移动虚拟示波器横游标到0V处,再移动另一根横游标到ΔV=1V处,得到与积分的曲线的交点,再移动虚拟示波器两根纵游标,从阶跃开始到曲线的交点,量得积分环节模拟电路时间常数Ti。

A6输出端(Uo)的实际响应曲线Uo(t)。

(3)改变时间常数,重新观测结果,填入实验报告。

4、观察比例积分环节的阶跃响应曲线
(1)打开虚拟示波器的单迹界面,点击开始,用示波器观测A6输出端(Uo)。

(2)待完整波形出来后,移动虚拟示波器横游标到1V(与输入相等)处,再移
动另一根横游标到ΔV=Kp×输入电压处,得到与积分曲线的两个交点。

(3)再分别移动示波器两根纵游标到积分曲线的两个交点,量得积分环节模拟电路时间常数Ti 。

(4)改变时间常数及比例系数,重新观测结果,填入实验报告。

5.观察比例微分环节的阶跃响应曲线
(1)打开虚拟示波器的界面,点击开始,用示波器观测系统的A6输出端(Uo ),响应曲线见图3-1-10。

等待完整波形出来后,把最高端电压(4.77V )减去稳态输出电压(0.5V ),然后乘以0.632,得到ΔV=2.7V 。

(2)移动虚拟示波器两根横游标,从最高端开始到ΔV=2.7V 处为止,得到与微分的指数曲线的交点,再移动虚拟示波器两根纵游标,从阶跃开始到曲线的交点,量得Δt=0.048S 。

(3)已知KD=10,则图3-1-9的比例微分环节模拟电路微分时间常数:
0.48S t K T D D =⨯=∆
R2/R1比例=5:2
6.观察PID (比例积分微分)环节的响应曲线
(1)打开虚拟示波器的单迹界面,点击开始,用示波器观测A6输出端(Uo )。

(2)等待完整波形出来后,移动虚拟示波器两根横游标使之ΔV=Kp ×输入电压,得到与积分的曲线的两个交点。

(3)再分别移动示波器两根纵游标到积分的曲线的两个交点,量得积分环节模拟电路时间常数Ti 。

(4)改变时间常数及比例系数,重新观测结果,填入实验报告。

五、实验数据及结果分析
1、比例环节
10
()()()O i U S R G S K U S R === (1)取1100R K =Ω 理论值:10
0.5 2o i R K U KU V R ==== 实际值:如图,1100R K =Ω时,我们发现测出的 2.07o U V =,与理论值2V 差距不大,在误差允许范围内理论值与测量值是相符合的。

(2)取1200R K =Ω 理论值:10
1 4o i R K U KU V R ==== 实际值:如图,1200R K =Ω时,我们发现测出的 3.95o U V =,与理论值4V 差距不大,在误差允许范围内理论值与测量值是相符合的。

2、惯性环节
110()(),,()1O i U S R K G S K R C U S S R ττ====+
取1200R K =Ω,1C F μ= 理论值:10
1 0.2R K s R τ=== 实际值:如图,1200R K =Ω时,我们发现测出的0.200s τ=,与理论值0.2s 相同,在误差允许范围内理论值与测量值是相符合的。

3、积分环节
0()1() ()O i U S G S R C U S S
ττ=== (1)取0100R K =Ω,1C F μ= 理论值:0.1s τ= 实际值:如图,100o R K =Ω时,我们发现测出的0.100s τ=,与理论值0.1s 相同,在误差允许范围内理论值与测量值是相符合的。

(2)取200o R K =Ω,1C F μ= 理论值:0.2s τ= 实际值:如图,200o R K =Ω时,我们发现测出的0.200s τ=,与理论值0.2s 相同,在误差允许范围内理论值与测量值是相符合的。

4、比例积分环节
110(S)1(S)(1),(S)O i U R G K K R C U S R ττ==+==
(1)取1200o R R K ==Ω,2C F μ=
理论值: 1 0.4K s τ==
实际值:如图,200o R K =Ω时,我们发现测出的0.410s τ=,与理论值0.4s 大致相同,在误差允许范围内理论值与测量值是相符合的。

(2)取1200o R R K ==Ω,0.5C F μ= 理论值: 1 0.1K s τ==
实际值:如图,200o R K =Ω时,我们发现测出的0.130s τ=,与理论值0.1s 大致相同,在误差允许范围内理论值与测量值是相符合的。

5、比例微分环节
121233120() D R R R R T T R C R C K K R R R ττ
+=+===+ 取0.5i U V =,0360R K =Ω,1180R K =Ω,2180R K =Ω,310R K =Ω,5C F μ= 理论值:1K =,=0.5T ,100.0050.05s τ=⨯=,10D K = 实际值:如图,0.5i U V =,0360R K =Ω,1180R K =Ω,2180R K =Ω,310R K =Ω时,我们发现测出的0.05s τ=,与理论值0.05s 相同,在误差允许范围内理论值与测量值是相符合的。

6、比例积分微分环节
121232121120
() () d i P R R R R T R C T R R C K R R R +=+=+=+ 取010R K =Ω,110R K =Ω,210R K =Ω,3200R =Ω,11C F μ=,21C F μ=
理论值:121()=0.02i T R R C s =+ 实际值: 我们可以从图中发现,实测值的0.10i T s =,误差较大。

六、实验误差分析
综合实验数据我们可以发现,在比例环节、惯性环节、积分环节和比例积分环节中,实验数据实测值与理论值误差极小,都处于5%以下。

但是最后一个比例积分微分环节(PID )中,误差很大。

分析误差较大的原因,大致有以下几个:
(1)实验由于微分的时间太短,较难捕捉到,把虚拟示波器的波形扩展到最大后但仍无法显示微分信号,因此,我们选择用一般的示波器观察,由于仪器之间存在干扰,因此存在较大的误差。

(2)实验所使用的电阻和电容实际值与标定值差距较大。

(3)仪器线路损坏。

七、实验小结
本次试验是第一次自动控制原理实验,因此需要我们前期较多的准备。

王老师首先给我们介绍了试验箱的使用方法和软件的操作指南,叫我们熟悉了实验的基本操作方法。

本次试验是典型环节的模拟研究,要求我们观察和分析各典型环节的阶跃响应曲线,了解各项电路参数对典型环节动态特性的影响。

在试验中,我们分别对六个典型环节进行试验,通过改变每个环节的电阻、电容等参数,以观察不同参数下的实验数据,并且与实际只作对比,得出最后的结论。

由于实验是验证性试验,因此本次试验比较容易,误差也较小。

有助于我们对理论知识进行对比和学习。

八、参考资料
1、《自动控制理论》,夏德钤主编,机械工业出版社
2、《自动控制理论》,胡寿松,航空工业出版社。

相关文档
最新文档