MATLAB基本语句及语法

1、精选优质文档-倾情为你奉上四川大学电气信息学院实验报告书实验名称： Matlab基本语句及语法 学生姓名 学 号 年级 指导教师 一、实验题目：实测电压的有效值及各次谐波的提取二、实验准备：（1）通常以电压整周期的方均根值来衡量电压的大小。电压方均根值的离散计算公式为：式中，N为一个周期内的采样点数； uk为第k点的电压瞬时值，V。（2）中国国家标准GB/T 17626.7-2008和IEC 61000-4-7:2002推荐了谐波的测量方法。推荐采用加矩形窗的DFT算法，窗口宽度应选为10个（50Hz系统）或12个（60Hz系统）周期。某一次谐波有效值以及它临近频谱分量有效值的方和根，如下式（

2、50Hz系统）：式中，Gg,n表示谐波群的有效值；Ck+i表示DFT输出频谱谱线分量的有效值；谐波数次n=k/N；k是傅里叶频谱分量的个数；N是时间窗截断信号周期数（50Hz系统中N=10），DFT输出频谱中各谱线的间隔为5Hz。3、 实验内容：1. 将.mat文件中或excel表格中的实测电压读入matlab。2. 读取ZQ77.mat文件，编制.m文件，绘出电压瞬时值和有效值随时间变化的曲线。绘制两张图形：第一张为实测电压随时间变化的图形，时间（ms），纵坐标为：电压（V），三相电压分别用红、蓝、绿色曲线（实线）表示。第二张为电压有效值图形，横坐标为：周波（cycle），纵坐标为：电压有效

3、值（V）。三相电压分别用红（实线）、蓝（虚线）、绿色（点线）曲线表示。（1） 实验程序：clc,clearnewData1 = load('-mat', 'ZQ77');vars = fieldnames(newData1);for i = 1:length(vars) assignin('base', varsi, newData1.(varsi);endplot(t,Va,'r');hold onplot(t,Vb,'b');hold onplot(t,Vc,'g');title('实测

4、电压随时间变化的波形');xlabel('时间(ms)');ylabel('电压(V)');figurefor i=0:21 Vasum=0;Vbsum=0;Vcsum=0; for j=i*96+1:1:(i+1)*96 Vasum=Vasum+Va(j)2;Vbsum=Vbsum+Vb(j)2;Vcsum=Vcsum+Vc(j)2; end Varms(i+1)=sqrt(Vasum/96);Vbrms(i+1)=sqrt(Vbsum/96);Vcrms(i+1)=sqrt(Vcsum/96);endcycle=1:22;plot(cycle,Var

5、ms,'r');hold onplot(cycle,Vbrms,'b-');hold onplot(cycle,Vcrms,'g.');hold ontitle('电压有效值随周波变化的波形');xlabel('周波（cycle）');ylabel('电压有效值(V)');（2） 实验结果：3. 读取“003.xls”中的数据，编制.m文件，根据IEC标准，用FFT计算实测的三相电流的各频谱分量有效值，提取实测电流的各次谐波。绘制一张各次谐波图（用subplot语句将该图片分为3个子图，三相电流分别

6、绘在三张子图中），含基频、3次、5次、7次、9次、11次、13次、15次、17次和19次谐波电流柱状图。横坐标为频次，纵坐标为：电流（A）。（1） 实验程序：clc,clear, , raw, dates = xlsread('C:UserskwyDesktoptemps003.xls','Sheet1','A2:G2561','',convertSpreadsheetExcelDates);raw = raw(:,2,3,4,5,6,7);dates = dates(:,1);data = reshape(raw:,size(r

7、aw);XData = datetime(dates:,1.', 'ConvertFrom', 'Excel');VaVolt = data(:,1);VbVolt = data(:,2);VcVolt = data(:,3);IaAmp = data(:,4);IbAmp = data(:,5);IcAmp = data(:,6);clearvars data raw dates;N=length(IaAmp);n=0:N-1;for n=1:2:19 k=n*10+1;sum=0; for j=-4:1:4 sum=sum+(Ybbb(k+j)2;

8、end Q(n)=sqrt(Ybbb(k-5)2/2+sum+(Ybbb(k+5)2/2);endsubplot(3,1,2)bar(Q);xlim(0 20);xlabel('频次');ylabel('电流(A)');title('B相电流的各次谐波');fs=12800;Ya=fft(IaAmp,N);Yb=fft(IbAmp,N);Yc=fft(IcAmp,N);f=n*fs/N;Yaa=Ya*2/N;Ybb=Yb*2/N;Ycc=Yc*2/N;subplot(2,1,1);Yaaa=abs(Yaa)/sqrt(2);Ybbb=abs(Yb

9、b)/sqrt(2);Yccc=abs(Ycc)/sqrt(2);for n=1:2:19. k=n*10+1;sum=0; for j=-4:1:4 sum=sum+(Yaaa(k+j)2; end G(n)=sqrt(Yaaa(k-5)2/2+sum+(Yaaa(k+5)2/2);endsubplot(3,1,1)bar(G);xlim(0 20);xlabel('频次');ylabel('电流(A)');title('A相电流的各次谐波');for n=1:2:19 k=n*10+1;sum=0; for j=-4:1:4 sum=sum+(Yccc(k+j)2; end W(n)=sqrt(Yccc(k-5)2/2+sum+(Yccc(k+5)2/2);endsubplot(3,1,3)bar(W);xlim(0 20);xlabel('频次');ylabel('电流(A)');title('C相电流的各次谐波');（2） 实验结果：4、 实验心得：