sinwave.m

电力系统中负序信号的Matlab仿真,包含simlink例程一个,负序仿真程序三个.

GAIN=[0 6 12 18 20 26 32 38];
REFOUT=1.25;            %AD REF Voltage :1.25V
AD73360MASTERCLK=8000;  %AD输入时钟kHz
AD73360DR=3;            %抽样频率分频位0-3 (1):DMCLK/1024=8kHz
AD73360SCD=3;           %串口分频位0-3 (0):DMCLK/8=1MHz
AD73360DMCLK=5;         %主时钟分频系数0-5 (0):DMCLK=8MHz
VOLTIGS=1;       %电压增益位0-8 (0):0dB
CURRIGS=1;       %电流增益位0-8 (0):0dB
%
compcoeff=20000;        %Compression coefficent;为适合HEX模式，浮点数转换成16位整型数
%fftdebug=2;             %调试模式1-2,
spfreq=6400;   %Sample Frequency ;抽样频率
fundaWave=50;           %Fundnamental Wave frqency ;基波频率
%hncf=[1 2 3 4 5 6 7]; %Harmnoic coefficent ;谐波系数表
%hnmg=0.1*[10 1 4 1 7 1 9]%Harmnoic magintude;谐波幅值表
%debug start JCL
hncf=[1 2 3 4 5 6 7]; %Harmnoic coefficent ;谐波系数表
hnmg=0.1*[10 1 5 1 3 2 4]%Harmnoic magintude;谐波幅值表
%debug end
time=20;                %Sample time ;采样时间 ms
t=0:1/spfreq:time/1000; 	%
i=diag(hnmg')*sin(2*pi*(fundaWave.*hncf')*t+pi*(0/2));
m=0:64;
i=compcoeff*sum(i);  	%
y=int16(i);
z=fft(i,128);
h=angle(z);
subplot(3,1,1);
plot(t,y);
subplot(3,1,2);
stem(abs(z));
subplot(3,1,3);
stem(h);
fid = fopen('sinWave.txt', 'wt');
fprintf(fid ,'%s\n','1651 2 0 0 0');%File Header *.dat
fprintf(fid ,'%d\n',y);
fed=fclose(fid);
copyfile('sinWave.txt','sinWave.dat');%The MATLAB can't driect out *.dat
afd=sprintf('%s','sinWave.dat out complete');