2021.9.12科研日志

一、工作内容

1.论文

今天躺平了一会，准备在仿真上做一做文章，所以推进了一点Matpower的文献。这篇名为《MATPOWER Steady-State Operations Planning and Analysis Tools for Power Systems Research and Education》基本上是所有使用Matpower做仿真的文献中必定会引用的文献，个人感觉可能躲不过，必须要了解一下，于是开了个头。

2.仿真

目前看的不少文献都是针对多通道的系统开发的算法，所以打算先把多通道的滤波器写出来，然后把最基本的卡方分布的检测写一下，在加点简单的攻击，以求把最简单的攻击检测用代码实现了，也便于在后面仿真的时候直接或者间接的套用这个框架。
9.13写完了最简单的不相关的两通道的传感器融合在这里。

二、心得

matlab中的plot的用法
1、画单个带点的黑色实线：

plot(t,x(t),'-k.');

其中后面单引号内的三个参数有以下可选项：

标记符	线型	标记符	点型	标记符	颜色
-	实线	+	加号	r	红
–	双划线	o	空心圆	g	绿
:	虚线	*	星号	b	蓝
:.	点划线	.	实心圆	c	蓝绿
		x	叉号	m	紫红
		s	正方形	y	黄
		d	菱形	k	黑
		^	上三角	w	白
		v	下三角
		>	右三角
		<	左三角
		p	五角星
		h	六边形

再举个例子，画个带上三角的蓝色虚线：

plot(t,x(t),':b^');

实际上，如果是一维的数据，可以不使用如t来做索引，比如上面的使用以下形式效果一样：

plot(x,'-k.');

2、在一张图片里绘制多副可以使用subplot命令：

figure;
subplot(2,2,1)
plot(w(1,:));xlabel('采样时间');ylabel('噪声');
title('第1状态值过程噪声');
subplot(2,2,2)
plot(w(2,:));xlabel('采样时间');ylabel('噪声');
title('第2状态值过程噪声');
subplot(2,2,3)
plot(v_1);xlabel('采样时间');ylabel('噪声');
title('第1状态值第1传感器量测噪声');
subplot(2,2,4)
plot(v_2);xlabel('采样时间');ylabel('噪声');
title('第2状态值第2传感器量测噪声');

3、标题使用title()命令

figure;
t=2:N;
plot(t,x(1,t),'-k.',t,z_1(t),'-b.',t,xkf_1(1,t),'-r.',t,xkf(1,t),'-g.');
legend('第1状态值','第1量测值','局部滤波器1估计值','全局滤波器估计值');
xlabel('采样时间');ylabel('位置');
title('局部滤波器1跟踪状态');

4、图片叠加 hold on,box on;用法

plot(t,x(2,t),'-k.',t,z_2(t),'-b.',t,xkf_2(2,t),'-r.',t,xkf(2,t),'-g.');

可以写成

hold on,box on;
plot(t,x(2,t),'-k.');
plot(t,z_2(t),'-b.');
plot(t,xkf_2(2,t),'-r.');
plot(t,xkf(2,t),'-g.');

三、未解决的问题

不相关的传感器的融合方法非常简单，但是融合后的提升非常有限，而且实际上攻击检测背景中的多通道系统也是存在大量冗余传感器的，言外之意就是传感器存在大量的耦合（相关），因此还是需要将相关传感器的联邦滤波写出来，才能实现攻击检测的功能。（直接剔除受攻击的传感器。）