微机继电保护matlab算法仿真(有源程序).doc
- 文档编号:30861205
- 上传时间:2024-05-02
- 格式:DOC
- 页数:7
- 大小:145KB
微机继电保护matlab算法仿真(有源程序).doc
《微机继电保护matlab算法仿真(有源程序).doc》由会员分享,可在线阅读,更多相关《微机继电保护matlab算法仿真(有源程序).doc(7页珍藏版)》请在冰豆网上搜索。
微
机
继
电
保
护
算
法
仿
真
电控学院
一.两点乘积算法仿真
(1)matlab中编写的程序
N=12;
t=(0:
0.02/N:
0.02)';
m=size(t);
y=sin(2*pi*50*t);
y1=[zeros(N/4,1);y(1:
m-N/4)];
ym=sqrt(y.^2+y1.^2);
subplot(3,2,1)
plot(t,y,'r.',t,y1,'xb');
legend('y(k)','y(k-T/4)');
title('两点乘积算法N=12');
subplot(3,2,2)
holdon;
plot(t,ym,'-r');
axis([0,0.02,0,1.2]);
xlabel('t/s');
ylabel('ym');
title('两点乘积算法算的有效值N=12');
text(0.01,0.6,'N=12');
N=16;
t=(0:
0.02/N:
0.02)';
m=size(t);
y=sin(2*pi*50*t);
y1=[zeros(N/4,1);y(1:
m-N/4)];
ym=sqrt(y.^2+y1.^2);
subplot(3,2,3)
plot(t,y,'r.',t,y1,'xb');
legend('y(k)','y(k-T/4)');
title('两点乘积算法N=16');
subplot(3,2,4)
holdon;
plot(t,ym,'-r');
axis([0,0.02,0,1.2]);
xlabel('t/s');
ylabel('ym');
title('两点乘积算法算的有效值N=16');
text(0.01,0.6,'N=16');
N=24;
t=(0:
0.02/N:
0.02)';
m=size(t);
y=sin(2*pi*50*t);
y1=[zeros(N/4,1);y(1:
m-N/4)];
ym=sqrt(y.^2+y1.^2);
subplot(3,2,5)
plot(t,y,'r.',t,y1,'xb');
legend('y(k)','y(k-T/4)');
title('两点乘积算法N=24');
subplot(3,2,6)
holdon;
plot(t,ym,'-r');
title('两点乘积算法算的有效值N=24');
axis([0,0.02,0,1.2]);
xlabel('t/s');
ylabel('ym');
text(0.01,0.6,'N=24');
(2)仿真出的波形
输出ym=sqrt(y.^2+y1.^2)
m=size(t);得到离散的两点
y1=[zeros(N/4,1);y(1:
m-N/4)]
N=12,即采样频率f=600HZ,原始信号离散化
输入信号y=sin(2*pi*50*t)
(3)流程图
二.傅里叶算法仿真
(1)matlab中编写的程序
T=0.02;
t1=0.02;
N=12;
Ts=T/N;
t=0:
Ts:
4*T;
y=(exp(-t/t1)-cos(2*pi*50*t));
subplot(3,2,1);
plot(t,y,'.r');
xlabel('t/ms');
ylabel('y(t)');
title('输入信号N=12');
a=1;
i=1:
N;
bs=sin(2*pi*i/N);
bc=cos(2*pi*i/N);
ys=filter(bs,a,y);
yc=filter(bc,a,y);
ym=2*abs(complex(ys,yc))/N;
subplot(3,2,2);
holdon
plot(t,ym)
xlabel('t/ms');
ylabel('ym');
title('傅里叶算法计算的有效值N=12');
holdon
T=0.02;
t1=0.02;
N=16;
Ts=T/N;
t=0:
Ts:
4*T;
y=(exp(-t/t1)-cos(2*pi*50*t));
subplot(3,2,3);
plot(t,y,'.r');xlabel('t/ms');
ylabel('y(t)');
title('输入信号N=16');
a=1;
i=1:
N;
bs=sin(2*pi*i/N);
bc=cos(2*pi*i/N);
ys=filter(bs,a,y);
yc=filter(bc,a,y);
ym=2*abs(complex(ys,yc))/N;
subplot(3,2,4);
plot(t,ym)
xlabel('t/ms');
ylabel('ym');
title('傅里叶算法计算的有效值N=16');
T=0.02;
t1=0.02;
holdon
N=24;
Ts=T/N;
t=0:
Ts:
4*T;
y=(exp(-t/t1)-cos(2*pi*50*t));
subplot(3,2,5);
plot(t,y,'.r');
xlabel('t/ms');
ylabel('y(t)');
title('输入信号N=24');
a=1;
i=1:
N;
bs=sin(2*pi*i/N);
bc=cos(2*pi*i/N);
ys=filter(bs,a,y);
yc=filter(bc,a,y);
ym=2*abs(complex(ys,yc))/N;
subplot(3,2,6);
plot(t,ym)
xlabel('t/ms');
ylabel('ym');
title('傅里叶算法计算的有效值N=24');
(2)仿真出的波形
(3)流ys=filter(bs,a,y);
yc=filter(bc,a,y);
bs=sin(2*pi*i/N);
bc=cos(2*pi*i/N);
输入信号y=(exp(-t/t1)-cos(2*pi*50*t))
N即每周波的采样点数,原始信号离散化
输出ym=2*abs(complex(ys,yc))/N
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 微机 保护 matlab 算法 仿真 有源 程序