基于matlab的自控课程设计

基于matlab的自控课程设计摘要串联滞后-超前校正兼有滞后校正和超前校正的优点，即已校正系统的响应速度较快，超调量较小，抑制高频噪声的性能也较好。当校正系统不稳定，且要求校正后系统的响应速度，相角裕度和稳态精度较高时，以采用串联滞后-超前校正为宜。其基本原理是利用滞后-超前网络的超前部分来增大系统的相角裕度，同时利用滞后部分来改善系统的稳态性能。此次课程设计就是利用MATLAB对一单位反馈系统进行滞后-超前校正。通过运用MATLAB的相关功能，绘制系统校正前后的伯德图、根轨迹和阶跃响应曲线，并计算校正后系统的时域性能指标。关键字：超前-滞后校正MATLAB伯德图时域性能指标目录摘要2设计要求及方法5设计要求5设计方法5滞后-超前校正设计目的和原理5滞后-超前设计目的5滞后=超前设计原理6滞后-超前校正的设计过程7校正前参数73.1.1用MATLAB绘制校正前系统的伯德图73.1.2用MATLAB求校正前系统的幅值裕量和相位裕量83.1.3用MATLAB绘制校正前系统的根轨迹9 3.1.4对校正前系统进行仿真分析 10滞后-超前校正设计参数计算113.2.1利用MATLAB计算出滞后校正器的传递函数113.2.2利用MATLAB计算出超前校正器的传递函数11滞后-超前校正的验证13用MATLAB求校正后系统的幅值裕量和相位裕量13用MATLAB求校正后系统的伯德图14用MATLAB求校正后系统的根轨迹15用MATLB对校正后的系统进行仿真分析16传递函数特征根的计算17校正前系统的传递函数的特征根17校后前系统的传递函数的特征根18系统动态分析18校正前系统的动态性能分析18校正后系统的动态性能分析23系统的根轨迹分析26校正前系统的根轨迹分析26校正后系统的根轨迹分析28系统的幅相特性30校正前系统的幅相特性30校正后系统的幅相特性31系统的对数幅频特性及对数相频特性32校正前系统的对数幅频特性及对数相频特性32校正后系统的对数幅频特性及对数相频特性33心得体会35参考文献351设计要求及方法1.1设计要求K已知单位负反馈系统的开环传递函数G(S) 0 ，试用频率法设计S(S1)(0.25S1)串联校正装置，要求校正后系统的静态速度误差系数K5s1，系统的相角裕度v450，校正后的剪切频率2rads。C1.2设计方法应用频率法确定滞后超前校正参数的步骤：根据稳态性能指标，绘制未校正系统的伯德图；选择校正后的截止频率；c确定校正参数；确定滞后部分的参数T；2确定超前部分的参数T；1将滞后部分和超前部分的传递函数组合在一起，即得滞后-超前校正的传递函数；绘制校正后的伯德图，检验性能指标。滞后-超前校正设计目的和原理2.1滞后-超前校正设计目的所谓校正就是在系统不可变部分的基础上，加入适当的校正元部件，使系统满足给定的性能指标。校正方案主要有串联校正、并联校正、反馈校正和前馈校正。确定校正装置的结构和参数的方法主要有两类：分析法和综合法。分析法是针对被校正系统的性能和给定的性能指标，首先选择合适的校正环节的结构，然滞后-超前校正的频域设计实际是超前校正和滞后校正频域法设计的综合，基本方法是利用滞后校正将系统校正后的穿越频率调整到超前部分的最大相角处的频率。具体方法是先合理地选择截止频率，先设计滞后校正部分，再根c据已经选定的设计超前部分。滞后-超前校正的设计过程3.1校正前系统的参数根据初始条件，调整开环传递函数：Gss1K00.25s1s根据系统的静态速度误差系数K5S1，取K=5S1。则满足初始条件的最小Kv V值时的开环传递函数为Gss150.25s1s3.1.1用MATLAB绘制校正前系统的伯德图程序：num=[5];den=[0.25,1.25,1,0];bode(num,den)grid得到的伯德图如图2所示。图2校正前系统的伯德图3.1.2用MATLAB求校正前系统的幅值裕量和相位裕量用命令margin(G)可以绘制出G的伯德图，并标出幅值裕量、相位裕量和对应的频率。用函数[kg,r,wg,wc]=margin(G)可以求出G的幅值裕量、相位裕量和幅值穿越频率。程序：num=[5];den=[0.25,1.25,1,0];G=tf(num,den);margin(G)[kg,r,wg,wc]=margin(G)得到的幅值裕量和相位裕量如图3所示。图3校正前系统的幅值裕量和相位裕量运行结果：kg=1.0000r=7.3342*10^-6wg=2.0000wc=2.0000即幅值裕量dBh02.65.0lg20相位裕量=(7.3342*10^-6)o穿越频率ωx=2rad/s截止频率ωcrad/s=23.1.3用MATLAB绘制校正前系统的根轨迹MATLAB中专门提供了绘制根轨迹的有关函数。[p,z]=pzmap(num,den)的功能是绘制连续系统的零、极点图。[r,k]=rlocus(num,den)的功能是绘制0k部分的根轨迹。程序：num=[5];den=[0.25,1.25,1,0];rlocus(num,den)得到校正前系统的根轨迹如图得到校正前系统的根轨迹如图4所示。图4校正前系统的根轨迹3.1.4对校正前系统进行仿真分析Simulink是可以用于连续、离散以及混合的线性、非线性控制系统建模、仿真和分析的软件包，并为用户提供了用方框图进行建模的图形接口，很适合于控制系统的仿真。仿真后得到的结果如图5和图6所示。图5校正前系统的仿真图图6校正前系统仿真的阶跃响应曲线3.2滞后-超前校正设计参数计算3.2.1利用MATLAB计算出滞后校正器的传递函数滞后校正器的传递函数为:TsTssGc11)(1根据题目要求，取校正后系统的剪切频率10,/3sradwc。并编写求滞后校正器传递函数的MATLAB程序MATLAB程序如下：wc=3;k0=5;nl=1;dl=conv(conv([1,0],[1,1]),[0.25,1]);beta=10;T=1/(0.1*wc);betat=beta*T;Gc1=tf([T,1],[betat,1])程序结果为：即对于校正后系统的滞后校正补偿器传递函数为：133.331333.3)(sssGc3.2.2利用MATLAB语言计算出超前校正器的传递函数要计算出校正后系统的传递函数，就编写求校正器的传递函数的MATLAB程序，其中调用了求超前校正器传递函数的函数leadc()，leadc.m保存在matlab7.0\work\文件夹下。若已知系统的开环传递函数与要求校正后的相角稳定裕度或剪切频率，求系统串联滞后校正器传递函数时，就可以调用此函数。leadc.m编制如下：function[Gc]=leadc(key,sope,vars)%MATLABFUNCTIONPROGRAMleadc.m%ifkey==1gama=vars(1);gama1=gama+5;[mag,phase,w]=bode(sope);[mu,pu]=bode(sope,w);gam=gama1*pi/180;alpha=(1-sin(gam))/(1+sin(gam));adb=20*log10(mu);am=10*log10(alpha);wc=spline(adb,w,am);T=1/(wc*sqrt(alpha));alphat=alpha*T;Gc=tf([T1],[alphat1]);elseifkey==2wc=vars(1);num=sope.num{1};den=sope.den{1};na=polyval(num,j*wc);da=polyval(den,j*wc);g=na/da;g1=abs(g);h=20*log10(g1);a=10^(h/10);wm=wc;T=1/(wm*(a)^(1/2));alphat=a*T;Gc=tf([T1],[alphat1]);elseifkey==3gama=vars(1);wc=vars(2);gama1=gama+5;num=sope.num{1};den=sope.den{1};ngv=polyval(num,j*wc);dgv=polyval(den,j*wc);g=ngv/dgv;thetag=angle(g);thetag_d=thetag*180/pi;mg=abs(g);gama_rad=gama1*pi/180;z=(1+mg*cos(gama_rad-thetag))/(-wc*mg*sin(gama_rad-thetag));p=(cos(gama_rad-thetag)+mg)/(wc*sin(gama_rad-thetag));nc=[z,1];dc=[p,1];Gc=tf(nc,dc);End其中key=1时，为var=gama，是根据要求校正后的相角稳定裕度计算滞后校正器；当key=2时，为var=wc，则是根据要求校正后的剪切频率计算校正器。MATLAB程序如下：wc=3;n1=conv([0,5],[3.333,1]);dl=conv(conv([1,0],[1,1]),conv([0.25,1],[33.33,1]));sope=tf(n1,dl);Gc]=leadc(2,sope,[wc])[程序结果为：即对于校正后系统的超前校正补偿器传递函数为：101412.01867.7)(sssGc综上，校正后系统的开环床底函数为：ssssssssGsGc2345259.344.42924.81177.0559.558.129)()(3.3滞后-超前校正后的验证由于校正过程中，多处采用的是近似计算，可能会造成滞后-超前校正后得到的系统的传递函数不满足题目要求的性能指标。所以需要对滞后-超前校正后的系统进行验证。下面用MATLAB求已校正系统的相角裕量和幅值裕量。3.3.1用MATLAB求校正后系统的幅值裕量和相位裕量程序：num=[129.8,55.59,5]den=[0.1177,8.924,42.4,34.59,1,0]G=tf(num,den);margin(G)[kg,r,wg,wc]=margin(G)得到的校正后系统的幅值裕量和相位裕量如图7所示。图7校正后系统的幅值裕量和相位裕量运行结果：kg=22.2950r=61.8620wg=18.1141wc=2.9766即校正后系统的幅值裕量dBh96.262950.22lg20相位裕量=61.862o穿越频率ωx=18.1141rad/s截止频率ωc=2.9766rad/s满足要求。（注：如果验证结果不满足要求，则需要重新选择校正后的截止频率，重复上述过程，直到满足要求为止。）3.3.2用MATLAB绘制校正后系统的伯德图程序：num=[129.8,55.59,5]den=[0.1177,8.924,42.4,34.59,1,0]bode(num,den)gridgrid得到的伯德图如图8所示。图8校正后系统的伯德图3.3.3用MATLAB绘制校正后系统的根轨迹程序：num=[129.8,55.59,5]den=[0.1177,8.924,42.4,34.59,1,0]rlocus(num,den)得到的校正后系统的根轨迹如图9所示。图9校正后系统的根轨迹3.3.4用MATLAB对校正前后的系统进行仿真分析用Simulink对校正后的系统仿真。仿真后得到的结果如图10和图11所示。图10校正后系统的仿真图图11校正后系统仿真的阶跃响应曲线4传递函数特征根的计算4.1校正前系统的传递函数的特征根校正前的开环传递函数为：)125.0)(1(5)(ssssG故该系统的闭环特征方程为：0525.125.023sss求特征根的MATLAB程序为：p=[0.25,1.25,1,5];roots(p)结果为：由于校正前系统单位负反馈的特征方程有右半平面的根，故校正前的闭环系统不稳定。4.24.2校正后系统的传递函数的特征根校正后的开环传递函数为：ssssssssGsGc2345259.344.42924.81177.0559.558.129)()(故该系统的闭环特征方程为：0559.5639.1644.42924.81177.02345sssss求特征根的MATLAB程序为p=[0.1177,8.924,42.4,164.39,56.59,5];roots(p)结果为：由于校正后系统单位负反馈的特征方程没有右半平面的根，故校正后的闭环系统稳定。5系统动态性能的分析5.1校正前系统的动态性能分析校正前的开环传递函数为：)125.0)(1(5)(ssssG⑴单位脉冲响应MATLAB程序为：num=[5];den=[0.25,1.25,1,0];s=tf(num,den);closep=feedback(s,1);figure(1)impulse(closep)title('校正前单位脉冲响应')图5-1校正前系统的单位脉冲响应⑵单位阶跃响应MATLAB程序为：num=[5];den=[0.25,1.25,1,0];s=tf(num,den);closep=feedback(s,1);figure(1)step(closep)title('校正前单位阶跃响应')图5-2校正前系统的单位阶跃响应由阶跃响应求动态性能参数要计算出阶跃响应动态性能参数，就编写求解阶跃响应动态性能参数的MATLAB程序，其中调用了函数perf()，perf.m保存在matlab7.0\work\文件夹下，其中key=1时，表示选择5%误差带，当key=2时表示选择2%误差带。y，t是对应系统阶跃响应的函数值与其对应的时间。函数返回的是阶跃响应超调量sigma(即σ)、峰值时间tp、调节时间ts。perf.m编制如下：function[sigma,tp,ts]=perf(key,y,t)%MATLABFUNCTIONPROGRAMperf.m%%Countsgmaandtp[mp,tf]=max(y);cs=length(t);yss=y(cs);sigma=(mp-yss)/ysstp=t(tf)%Counttsi=cs+1;n=0;whilen==0,i=i-1;ifkey==1,ifi==1,n=1;elseify(i)>1.05*yss,n=1;end;elseifkey==2,ifi==1,n=1;elseify(i)>1.02*yss,n=1;end;endend;t1=t(i);cs=length(t);j=cs+1;n=0;whilen==0,j=j-1;ifkey==1,ifj==1,n=1;elseify(j)<0.95*yss,n=1;end;elseifkey==2,ifj==1,n=1;elseify(j)<0.98*yss,n=1;end;end;end;t2=t(j);ift2<tpift1>t2;ts=t1endelseift2>tp,ift2<t1,ts=t2elsets=t1endendMATLABMATLAB程序为：globalyts1=tf([0,5],[0.25,1.25,1,0]);sys=feedback(s1,1);figure(1)step(sys)y,t]=step(sys);[perf(1,y,t)程序结果为：所以σ%=1.9267%，tp=8.3897s，ts=8.7124s。在Simulink窗口里菜单方式下的单位斜坡响应的动态结构图如下：图5-3校正前系统的单位斜坡响应的动态结构仿真图校正前单位斜坡响应曲线如下所示：图5-4校正前系统的单位斜坡响应5.2校正后系统的动态性能分析校正后的开环传递函数为：ssssssssGsGc2345259.344.42924.81177..0559.558.129)()(⑴单位脉冲响应MATLAB程序为：num=[129.8,55.59,5]den=[0.1177,8.924,42.4,34.59,1,0]s=tf(num,den);closep=feedback(s,1);figure(1)impulse(closep)title('校正后单位脉冲响应')图5-5校正后系统的单位脉冲响应⑵单位阶跃响应MATLAB程序为：num=[129.8,55.59,5]den=[0.1177,8.924,42.4,34.59,1,0]s=tf(num,den);closep=feedback(s,1);figure(1)step(closep)title('校正后单位阶跃响应')')')图5-6校正后系统的单位阶跃响应由阶跃响应求动态性能参数MATLAB程序为：globalytnum=[120.1,70.02,10]den=[0.5063,12.05,49.84,39.3,1,0]s=tf(num,den);sys=feedback(s,1);figure(1)step(sys)[y,t]=step(sys);pref(1,y,t)程序结果为：所以σ%=0.0568%，tp=1.0394s，ts=1.0394s。由图可知，系统稳态误差011)(1hess显然，显然，校正后的系统动态性能比校正前的系统动态性能好得多，改善了系统的平稳性、快速性和稳态精度。⑶单位斜坡响应在Simulink窗口里菜单方式下的单位斜坡响应的动态结构图如下：图5-7校正后系统的单位斜坡响应的动态结构仿真图校正前单位斜坡响应曲线如下所示：图5-8校正后系统的单位斜坡响应单位脉冲、阶跃、斜坡响应曲线的相互对应关系是：单位脉冲响应的积分是单位阶跃响应曲线，单位阶跃响应的积分是单位斜坡响应。6系统的根轨迹分析6.1校正前系统的根轨迹分析校正前的开环传递函数为：)125.0)(1(5)(ssssGMATLABMATLAB程序如下：num=[5];den=[0.25,1.25,1,0];rlocus(num,den)图6-1校正前系统的根轨迹确定分离点坐标：图6-2校正前根轨迹的分离点坐标分离点坐标d=-0.465确定与虚轴交点的坐标：确定与虚轴交点的坐标：图6-3校正前根轨迹与虚轴的交点坐标与虚轴的交点坐标：i21、i22。K*=5.02校正前系统稳定时增益K*的变化范围是K*<5.02，而题目中的K*=10，校正前的闭环系统不稳定。6.2校正后系统的根轨迹分析校正后的开环传递函数为：ssssssssGsGc2345259.344.42924.81177.0559.558.129)()(MATLAB程序如下：num=[129.8,55.59,5]den=[0.1177,8.924,42.4,34.59,1,0]rlocus(num,den)图6-4校正后系统的根轨迹确定分离点坐标：图6-5校正后根轨迹的分离点坐标分离点坐标d=-2.4确定与虚轴交点的坐标：确定与虚轴交点的坐标：图6-6校正后根轨迹与虚轴的交点坐标与虚轴的交点坐标：i9.171、i9.172。增益K*=100校正前系统稳定时增益K*的变化范围是K*<90.3，而题目中的K*=5，校正后的闭环系统稳定。7系统的幅相特性7.1校正前系统的幅相特性校正前的开环传递函数为：)125.0)(1(5)(ssssGMATLAB程序如下：num=[5];den=[0.25,1.25,1,0];s1=tf(num,den);nyquist(s1)图7-1校正前系统的Nyquist曲线由于开环传递函数中含有一个积分环节，所以从0到0顺时针补画一圈，再由上图可知，Nyquist曲线顺时针围绕点(-1，j0)两圈，所以,2R，而0P，所以校正前闭环系统不稳定。7.2校正后系统的幅相特性校正后的开环传递函数为：ssssssssGsGc2345259.344.42924.81177.0559.558.129)()(MATLAB程序如下：num=[120.1,70.2,10];den=[0.5063,12.05,49.84,39.3,1,0];s=tf(num,den);nyquist(s)图7-2校正后系统的Nyquist曲线由于开环传递函数中含有一个积分环节，所以从0到0顺时针补画一圈，再由上图可知，Nyquist曲线顺时针围绕点(-1，j0)0圈，所以,0R，而0P，所以000RPZ，所以校正后闭环系统稳定。8系统的对数幅频特性及对数相频特性8.1校正前系统的对数幅频特性及对数相频特性校正前的开环传递函数为：)125.0)(1(5)(ssssGMATLAB程序如下：MATLAB程序如下：num=[5];den=[0.25,1.25,1,0];G=tf(num,den);margin(G)[[kg,r,wg,wc]=margin(G)得到的幅值裕量和相位裕量如图3所示。图3校正前系统的幅值裕量和相位裕量运行结果：kg=1.0000r=7.3342*10^-6wg=2.0000wc=2.0000即幅值裕量dBh02.65.0lg20相位裕量=(7.3342*10^-6)o穿越频率ωx=2rad/s截止频率ωc=2rad/s由于开环传递函数中含有一个积分环节，应从对数相频特性曲线上0处开始，用虚线向上补画90角。所以，在开环对数幅频dB0)(L的频率范围内，对应的开环对数相频特性曲线)(对线有一次负穿越，即1N，则2)(2NNR，2)2(0RPZ，所以校正前闭环系统不稳定。8.2校正后系统的对数幅频特性及对数相频特性校正后的开环传递函数为：sssssssssGsGc2345259.344.42924.81177.0559.558.129)()(MATLAB程序如下：MATLAB程序如下：num=[129.8,55.59,5]den=[0.1177,8.924,42.4,34.59,1,0]G=tf(num,den);margin(G)[kg,r,wg,wc]=margin(G)得到的校正后系统的幅值裕量和相位裕量如图7所示。图7校正后系统的幅值裕量和相位裕量运行结果：kg=22.2950r=61.8620wg=18.1141wc=2.9766即校正后系统的幅值裕量dBh96.262950.22lg20相位裕量=61.862o穿越频率ωx=18.1141rad/s截止频率ω=2.9766rad/sc由于开环传递函数中含有一个积分环节，应从对数相频特性曲线上0处开始，用虚线向上补画90角。所以，在开环对数幅频L()0dB的频率范围，对应的开环对数相频特性曲线()对线无穿越，即N0、N0，则R2(NN)0，ZPR000，。 并且校正后的系统的相位裕量61.862o≥45o、静态速度误差系数K=5≥2、剪切频率w2rad/s均符合题目设计的要求，故系统校正到此全部完。 c9心得体会此次课程设计的内容对一个单位反馈系统进行滞后-超前校正。回顾此次实践的整个过程，虽然只有短短的几天，但是真的在这个自己独立学习的过程中学到了好多东西。课程设计开始阶段比较顺利，但是做到计算校正后系统的时域性能指标这里时，遇到了不小麻烦，不会用MATLAB编程得校正参数。后来，在同学的帮助下，找到了计算校正参数的那段程序，于是就求出了校正后的传递函数。这次的课程设计，让我更更深一步的了解MATLAB这个十分有用的软件，与理论相结合。与此同时，通过此次课程设计，让我加深了对系统进行滞后-超前设计过程的理解，还掌握了用MATLAB编程计算系统时域性能指标和系统幅值裕量、相位裕量的方法。总而言之，这次的课程设计的确让我受益匪浅。参考文献[1]陈杰主.MATLAB宝典.电子工业出版.2007.[2]李国勇.智能控制与其MATLAB实现.电子工业出版社.2005.[3]胡寿松.自动控制原理.科学出版社.2007.[4]王万良.自动控制原理.高等教育出版社.2008.[5]黄坚.自动控制原理及其应用.高等教育出版社.2004.[6]程鹏.自动控制原理[M].北京：高等教育出版社，2009[7]徐薇莉.自动控制理论与设计[M]上海：上海交通大学出版社,2001[8]欧阳黎明.MATLAB控制系统设计[M].北京：国防工业出版社,2001[9]吴天明，谢小竹等.MATLAB电力系统设计与分析.国防工业出版社，2004毕业设计（论文）原创性声明和使用授权说明原创性声明本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。作者签名：日期：指导教师签名：日期：使用授权说明本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。作者签名：日期：

学位论文原创性声明本人郑重声明：所呈交的论文是本人在导师的指导下独立进行研究所取得的研究成果。除了文中特别加以标注引用的内容外，本论文不包含任何其他个人或集体已经发表或撰写的成果作品。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。本人完全意识到本声明的法律后果由本人承担。作者签名：日期：年月日学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定，同意学校保留并向国家有关部门或机构送交论文的复印件和电子版，允许论文被查阅和借阅。本人授权大学可以将本学位论文的全部或部分内容编入有关数据库进行检索，可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文。涉密论文按学校规定处理。作者签名：日期：年月日导师签名：日期：年月日指导教师评阅书建议成绩：建议成绩：□优□良□中□及格□不及格（在所选等级前的□内画“√”）指导教师：（签名）单位：（盖章）年月日

评阅教师评阅书评阅教师评价：一、论文（设计）质量1、论文（设计）的整体结构是否符合撰写规范？□优□良□中□及格□不及格2、是否完成指定的论文（设计）任务（包括装订及附件）？□优□良□中□及格□不及格二、论文（设计）水平1、论文（设计）的理论意义或对解决实际问题的指导意义□优□良□中□及格□不及格2、论文的观念是否有新意？设计是否有创意？□优□良□中□及格□不及格3、论文（设计说明书）所体现的整体水平□优□良□中□及格□不及格建议成绩：□优□良□中□及格□不及格（在所选等级前的□内画“√”）评阅教师：（签名）单位：（盖章）年月日教研室（或答辩小组）及教学系意见 --教学系意见：教学系意见：系主任：（签名）年月日

学位论文原创性声明本人郑重声明：所呈交的学位论文，是本人在导师的指导下进行的研究工作所取得的成果。尽我所知，除文中已经特别注明引用的内容和致谢的地方外，本论文不包含任何其他个人或集体已经发表或撰写过的研究成果。对本文的研究做出重要贡献的个人和集体，均已在文中以明确方式注明并表示感谢。本人完全意识到本声明的法律结果由本人承担。学位论文作者（本人签名）：年月日学位论文出版授权书本人及导师完全同意《中国博士学位论文全文数据库出版章程》、《中国优秀硕士学位论文全文数据库出版章程》(以下简称“章程”)，愿意将本人的学位论文提交“中国学术期刊（光盘版）电子杂志社”在《中国博士学位论文全文数据库》、《中国优秀硕士学位论文全文数据库》中全文发表和以电子、网络形式公开出版，并同意编入CNKI《中国知识资源总库》，在《中国博硕士学位论文评价数据库》中使用和在互联网上传播，同意按“章程”规定享受相关权益。论文密级：□公开□保密（___年__月至__年__月）(保密的学位论文在解密后应遵守此协议)作者签名：_______导师签名：___