MATLAB基本语句及语法.docVIP

四川大学电气信息学院 实 验 报 告 书 实验名称： Matlab基本语句及语法 　 学生姓名 　 学 号 年级 指导教师 一、实验题目： 实测电压的有效值及各次谐波的提取 二、实验准备： （1）通常以电压整周期的方均根值来衡量电压的大小。电压方均根值的离散计算公式为： 式中，N为一个周期内的采样点数； uk为第k点的电压瞬时值，V。 （2）中国国家标准GB/T 17626.7-2008和IEC 61000-4-7:2002推荐了谐波的测量方法。推荐采用加矩形窗的DFT算法，窗口宽度应选为10个（50Hz系统）或12个（60Hz系统）周期。某一次谐波有效值以及它临近频谱分量有效值的方和根，如下式（50Hz系统）： 式中，Gg,n表示谐波群的有效值；Ck+i表示DFT输出频谱谱线分量的有效值；谐波数次n=k/N；k是傅里叶频谱分量的个数；N是时间窗截断信号周期数（50Hz系统中N=10），DFT输出频谱中各谱线的间隔为5Hz。 实验内容： 将.mat文件中或excel表格中的实测电压读入matlab。 读取ZQ77.mat文件，编制.m文件，绘出电压瞬时值和有效值随时间变化的曲线。绘制两张图形：第一张为实测电压随时间变化的图形，时间（ms），纵坐标为：电压（V），三相电压分别用红、蓝、绿色曲线（实线）表示。第二张为电压有效值图形，横坐标为：周波（cycle），纵坐标为：电压有效值（V）。三相电压分别用红（实线）、蓝（虚线）、绿色（点线）曲线表示。 实验程序： clc,clear newData1 = load(-mat, ZQ77); vars = fieldnames(newData1); for i = 1:length(vars) assignin(base, vars{i}, newData1.(vars{i})); end plot(t,Va,r); hold on plot(t,Vb,b); hold on plot(t,Vc,g); title(实测电压随时间变化的波形); xlabel(时间(ms)); ylabel(电压(V)); figure for 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); end cycle=1:22; plot(cycle,Varms,r);hold on plot(cycle,Vbrms,b--);hold on plot(cycle,Vcrms,g.);hold on title(电压有效值随周波变化的波形); xlabel(周波（cycle）); ylabel(电压有效值(V)); 实验结果： 读取“003.xls”中的数据，编制.m文件，根据IEC标准，用FFT计算实测的三相电流的各频谱分量有效值，提取实测电流的各次谐波。绘制一张各次谐波图（用subplot语句将该图片分为3个子图，三相电流分别绘在三张子图中），含基频、3次、5次、7次、9次、11次、13次、15次、17次和19次谐波电流柱状图。横坐标为频次，纵坐标为：电流（A）。 实验程序： clc,clear [~, ~, raw, dates] = xlsread(C:\Users\kwy\Desktop\temps\003.xls,Sheet1,A2:G2561,,@convertSpreadsheetExcelDates); raw = raw(:,[2,3,4,5,6,7]); dates = dates(:,1); data = reshape([raw{:}],size(raw)); 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