微电网鲁棒性研究（Matlab代码实现）

原创

荔枝科研社 2022-08-16 00:53:52 ©著作权

文章标签 算法 matlab matlab代码参考文献权重 文章分类 代码人生

©著作权归作者所有：来自51CTO博客作者荔枝科研社的原创作品，请联系作者获取转载授权，否则将追究法律责任

💥💥💥💞💞💞欢迎来到本博客❤️❤️❤️💥💥💥

📝目前更新：🌟🌟🌟电力系统相关知识，期刊论文，算法，机器学习和人工智能学习。
🚀支持：🎁🎁🎁如果觉得博主的文章还不错或者您用得到的话，可以关注一下博主，如果三连收藏支持就更好啦！这就是给予我最大的支持！

🎉🎉欢迎您的到来🎉🎉

⛅⛅⛅ 📃个人主页：科研室🌈🌈🌈

📚📚📚📋所有代码目录：电气工程科研社👨‍💻👨‍💻👨‍💻

【现在公众号名字改为：荔枝科研社】
👨‍🎓博主课外兴趣：中西方哲学，送予读者：
👨‍💻做科研，涉及到一个深在的思想系统，需要科研者逻辑缜密，踏实认真，但是不能只是努力，很多时候借力比努力更重要，然后还要有仰望星空的创新点和启发点。当哲学课上老师问你什么是科学，什么是电的时候，不要觉得这些问题搞笑，哲学就是追究终极问题，寻找那些不言自明只有小孩子会问的但是你却回答不出来的问题。在我这个专栏记录我有空时的一些哲学思考和科研笔记：科研和哲思。建议读者按目录次序逐一浏览，免得骤然跌入幽暗的迷宫找不到来时的路，它不足为你揭示全部问题的答案，但若能让人胸中升起一朵朵疑云，也未尝不会酿成晚霞斑斓的别一番景致，万一它居然给你带来了一场精神世界的苦雨，那就借机洗刷一下原来存放在那儿的“真理”上的尘埃吧。
或许，雨过云收，神驰的天地更清朗.......🔎🔎🔎
📋📋📋本文目录如下：⛳️⛳️⛳️

目录
1 概述
2 微电网
3 运行结果
4 Matlab代码实现
5 参考文献
6 写在最后

1 概述

在电工装备的制造、装配、使用的过程中,不可避免地存在诸如制造偏差、材料属性差异、工作环境变化等不确定因素。不确定因素往往使得电工装备的实际参数偏离原有的设计值,从而使产品性能下降,甚至发生事故。为了缓解不确定性因素对电工装备造成的负面影响,同时平衡电工装备的可靠性、鲁棒性和其他性能指标之间的关系,因此需要在电工产品的优化设计阶段考虑工程中的不确定性。在考虑鲁棒性和可靠性的前提下,本文为了解决电气工程领域检验电磁场分析方法的标准问题、电工装备中的永磁魔环以及永磁同步电机中的优化问题,主要完成以下工作。首先,在不确定性因素影响下,为了保证电工装备的可靠性与鲁棒性,本文研究了一种科学的最优权重辅助的兼顾可靠性和鲁棒性的优化设计方法（ω-RBRDO）。这种方法兼顾了可靠性和鲁棒性两种指标,通过使约束条件在一定概率下满足,来保证所寻求设计方案的可靠性,同时为了寻求使鲁棒性指标及目标性能这两个指标最优的比例系数,引入了最优权重的思想。

2 微电网

这一部分讲来很多，感兴趣可以搜索：

微电网鲁棒性研究（Matlab代码实现）_matlab_02

3 运行结果

微电网鲁棒性研究（Matlab代码实现）_参考文献_03

微电网鲁棒性研究（Matlab代码实现）_matlab代码_04

微电网鲁棒性研究（Matlab代码实现）_算法_05

微电网鲁棒性研究（Matlab代码实现）_参考文献_06

微电网鲁棒性研究（Matlab代码实现）_权重_07

4 Matlab代码实现

本文仅展现部分代码，全部代码见：🍞正在为您运送作品详情

%% 灵敏度函数的奇异值
figure(1)
S = (eye(2,2) - P0*K_s)^-1;
sigma(minreal(ss(S)))
title('S 的奇异值')
grid on


%% 互补灵敏度函数的奇异值
figure(2)
T = P0*K_s*(eye(2,2) - P0*K_s)^-1;
sigma(minreal(ss(T)))
title('T 的奇异值')
grid on
%saveas(gcf,'T_sigma.png')

%% M11的奇异值
figure(3)
M = lft(Q_ss,K_s);
M11 = M(1:2,1:2);
sigma(minreal(ss(M11)))
grid on
title('M11的奇异值')
%saveas(gcf,'M11_sigma.png')

%% M的奇异值
figure(4)
sigma(minreal(lft(Q,K_s)))
title('M的奇异值')
grid on
%saveas(gcf,'M_sigma.png')

%% 鲁棒稳定性
rbst_s = Wa*K_s*(eye(2,2) - P0*K_s)^-1;
disp('鲁棒稳定性')
norm(rbst_s,inf)

%% 性能
nmnl_p = Ws*((eye(2)+P0*(eye(2,2) - P0*K_s))^-1);
disp('性能')
norm(nmnl_p,inf)

%% 鲁棒性能
rbst_p = [rbst_s;nmnl_p];
disp('鲁棒性能')
norm(rbst_p,inf)

%% 输出
A_k = K_ss.A;
B_k = K_ss.B;
C_k = K_ss.C;
D_k = K_ss.D;
save('shuju','A_k','B_k','C_k','D_k')
return