簡(jiǎn)介Power System ToolboxPST是一套運(yùn)行于MATLAB環(huán)境的電力系統(tǒng)專(zhuān)業(yè)工具箱面向科研人員、工程師與學(xué)生用于解決復(fù)雜電力網(wǎng)絡(luò)的建模、穩(wěn)態(tài)/動(dòng)態(tài)/暫態(tài)仿真及控制策略驗(yàn)證問(wèn)題。壓縮包共有211個(gè)文件其中209個(gè)m腳本/函數(shù)和2個(gè)pdf文檔整體僅1.1MBm文件涵蓋設(shè)備模型、仿真主程序、示例數(shù)據(jù)與繪圖工具pdf文檔提供使用與部署說(shuō)明。核心功能支持發(fā)電機(jī)、變壓器、線路、負(fù)荷、斷路器、繼電保護(hù)裝置等模型詳細(xì)程度可從簡(jiǎn)化靜態(tài)到復(fù)雜動(dòng)態(tài)靈活選擇能夠模擬正常運(yùn)行時(shí)的電壓、功率分布及故障后的系統(tǒng)穩(wěn)定性。工具箱還提供勵(lì)磁控制器、調(diào)速器、無(wú)功補(bǔ)償器等控制器的設(shè)計(jì)與測(cè)試環(huán)境配合波形圖、功率流圖、頻率響應(yīng)曲線等可視化圖表并支持用戶(hù)自定義模型和Simulink集成便于新能源并網(wǎng)、保護(hù)配置優(yōu)化等應(yīng)用擴(kuò)展。目前已有761人學(xué)習(xí)下載借助內(nèi)置腳本和示例數(shù)據(jù)電力系統(tǒng)專(zhuān)業(yè)師生或從業(yè)者可快速建立仿真環(huán)境開(kāi)展教學(xué)實(shí)驗(yàn)與課題研究。1. Power System Toolbox 到底是干什么的三個(gè)問(wèn)題幫你確認(rèn)它值不值得裝Power System Toolbox 這類(lèi) MATLAB 仿真工具集最容易被誤解成“又一個(gè)畫(huà)潮流圖的插件”。實(shí)際上它是一套面向電力系統(tǒng)動(dòng)態(tài)分析的命令行框架能算潮流、能搭動(dòng)態(tài)模型、能跑暫態(tài)穩(wěn)定還能把你自己的控制策略寫(xiě)成一段回調(diào)函數(shù)塞進(jìn)仿真循環(huán)里。所以先別急著找圖形界面——如果你關(guān)心的是“我這個(gè)控制器掛到雙機(jī)系統(tǒng)上故障切除后轉(zhuǎn)子還能不能穩(wěn)住”那它就比一堆只做靜態(tài)潮流計(jì)算的工具箱對(duì)路得多。它適合三類(lèi)人做機(jī)電暫態(tài)算法驗(yàn)證的研究生做保護(hù)或控制策略研發(fā)的工程師以及要在論文里復(fù)現(xiàn)“某個(gè)典型系統(tǒng)”的仿真曲線、又不想從零手寫(xiě)微分方程求解器的開(kāi)發(fā)者。不適合只想要一個(gè)“點(diǎn)一下出一張單線圖”的可視化工具的人。反直覺(jué)的一點(diǎn)是這套工具真正值錢(qián)的地方不是它自帶的幾個(gè)測(cè)試系統(tǒng)而是那套“矩陣定義系統(tǒng)、函數(shù)驅(qū)動(dòng)計(jì)算、回調(diào)注入控制”的松耦合結(jié)構(gòu)。你的核心資產(chǎn)是算例數(shù)據(jù)和模型參數(shù)不是界面。2. 從數(shù)據(jù)格式入手母線、支路、發(fā)電機(jī)矩陣怎么排才不出錯(cuò)2.1 bus 矩陣母線編號(hào)、電壓和功率注入怎么排這類(lèi)工具箱最通用的數(shù)據(jù)組織方式是把整個(gè)電力系統(tǒng)拆成若干組固定列數(shù)的矩陣。第一組是母線矩陣每一行代表一條母線每一列對(duì)應(yīng)一個(gè)物理量。列的順序一般沿襲一套歷史約定母線編號(hào)、電壓幅值初值、相角初值、有功負(fù)荷、無(wú)功負(fù)荷、并聯(lián)導(dǎo)納、所屬區(qū)域編號(hào)、節(jié)點(diǎn)類(lèi)型以及電壓約束上下限。% 構(gòu)造一個(gè)三節(jié)點(diǎn)測(cè)試系統(tǒng)的 bus 矩陣 % 列含義編號(hào), V0(pu), theta0(rad), P_load(pu), Q_load(pu), Gsh(pu), Bsh(pu), 區(qū)域, Vmax, Vmin bus [ 1 1.00 0.00 1.20 0.40 0 0 1 1.10 0.90; % 節(jié)點(diǎn)1PQ負(fù)荷1.20.4j 2 1.00 0.00 0.50 0.20 0 0 1 1.10 0.90; % 節(jié)點(diǎn)2小負(fù)荷 3 1.00 0.00 0.30 0.10 0 0 1 1.10 0.90; % 節(jié)點(diǎn)3接發(fā)電機(jī) ];這里每個(gè)數(shù)字都是標(biāo)幺值pu基準(zhǔn)容量一般取全系統(tǒng)統(tǒng)一的 MVA 值比如 100MVA。第 1、2 列的電壓和相角是潮流迭代的初值不是約束值真正的電壓約束在最后兩列。新手最容易踩的坑是把“初值”當(dāng)成“給定值”導(dǎo)致發(fā)電機(jī)母線電壓永遠(yuǎn)迭代不動(dòng)。2.2 line 矩陣R/X/B/變比對(duì)應(yīng)哪一列第二組是支路矩陣每條線路或每臺(tái)變壓器一行。列順序有很強(qiáng)的一致性支路編號(hào)、首端母線號(hào)、末端母線號(hào)、串聯(lián)電阻 R、串聯(lián)電抗 X、對(duì)地電納 B、變壓器變比以及相移。單位同樣是標(biāo)幺值R/X 是基于系統(tǒng)基準(zhǔn)值的串聯(lián)阻抗B 是線路充電電納的一半還是全量不同工具箱習(xí)慣不一但最常見(jiàn)的是全量電納除以 2 再填入。% 三節(jié)點(diǎn)系統(tǒng)的三條支路 % 列含義編號(hào), from, to, R(pu), X(pu), B/2(pu), 變比k, 相移(度) line [ 1 1 2 0.020 0.060 0.030 0 0.0; % 線路1-2 2 2 3 0.030 0.090 0.040 0 0.0; % 線路2-3 3 3 1 0.040 0.120 0.050 1 0.0; % 變壓器3-1變比1.0 ];注意變壓器和線路混在同一張表里時(shí)變比列的填法有講究如果是標(biāo)幺變比 1.0直接填 1如果是非標(biāo)幺變比比如 10.5/121kV 的實(shí)際抽頭需要折算到基準(zhǔn)電壓之后再填。X 為負(fù)不一定是你填錯(cuò)了可能是串補(bǔ)電容但如果 R 和 X 填反潮流會(huì)直接給你一個(gè)“換相失敗”式的離譜相角。2.3 發(fā)電機(jī)與動(dòng)態(tài)模型矩陣潮流模型和動(dòng)態(tài)模型的銜接發(fā)電機(jī)矩陣是這兩類(lèi)模型之間的橋梁。潮流計(jì)算只需要發(fā)電機(jī)母線號(hào)、有功出力、無(wú)功出力、有功上下限、無(wú)功上下限這幾個(gè)字段而暫態(tài)穩(wěn)定仿真還需要給出發(fā)電機(jī)的動(dòng)態(tài)參數(shù)比如慣性常數(shù)、暫態(tài)電抗、阻尼系數(shù)以及勵(lì)磁系統(tǒng)和調(diào)速器對(duì)應(yīng)的模型編號(hào)。常見(jiàn)做法是分兩張表一張管潮流注入一張管動(dòng)態(tài)參數(shù)中間用母線號(hào)關(guān)聯(lián)。% 發(fā)電機(jī)潮流數(shù)據(jù)母線號(hào), Pg(pu), Qg(pu), Pmax, Pmin, Qmax, Qmin gen [ 3 0.80 0.20 1.20 0.00 0.80 -0.40; % 節(jié)點(diǎn)3的平衡機(jī) ]; % 發(fā)電機(jī)動(dòng)態(tài)參數(shù)母線號(hào), H(s), Xd(pu), Xq(pu), D, 模型階數(shù) mac [ 3 6.50 0.30 0.35 2.0 2; % 2階經(jīng)典模型 ];mac 表的模型階數(shù)直接決定后面積分器要不要解勵(lì)磁繞組暫態(tài)。2 階經(jīng)典模型只保留轉(zhuǎn)子運(yùn)動(dòng)方程適合看第一擺穩(wěn)定性想研究阻尼特性至少要 4 階以上模型把勵(lì)磁回路動(dòng)態(tài)算進(jìn)去。如果之后發(fā)現(xiàn)暫態(tài)波形衰減得過(guò)分平滑先回來(lái)檢查 D 是不是被填成了 0。3. 最小潮流計(jì)算跑通三組命令把結(jié)果拿到手3.1 準(zhǔn)備算例文件把三組矩陣寫(xiě)進(jìn)工作區(qū)這類(lèi)工具箱通常不要求你把數(shù)據(jù)存成特定格式的文件只要在 MATLAB 工作區(qū)里有bus、line、gen三組矩陣就能直接開(kāi)工。最穩(wěn)妥的落地方式是把系統(tǒng)數(shù)據(jù)寫(xiě)進(jìn)一個(gè)獨(dú)立的腳本文件每次運(yùn)行先執(zhí)行這個(gè)腳本再調(diào)求解器。這樣算例與求解分離換系統(tǒng)時(shí)只替換數(shù)據(jù)文件不碰主循環(huán)。% case3.m —— 三節(jié)點(diǎn)測(cè)試系統(tǒng)數(shù)據(jù)腳本 bus [ ... ]; % 粘貼 2.1 節(jié)的三行母線數(shù)據(jù) line [ ... ]; % 粘貼 2.2 節(jié)的三行支路數(shù)據(jù) gen [ ... ]; % 粘貼 2.3 節(jié)的發(fā)電機(jī)數(shù)據(jù)腳本文件沒(méi)有返回值運(yùn)行之后矩陣直接進(jìn)工作區(qū)后續(xù)函數(shù)按變量名取數(shù)據(jù)。變量名必須嚴(yán)格一致否則求解器會(huì)報(bào)“找不到匹配的 bus 變量”。養(yǎng)成一個(gè)習(xí)慣任何從外部導(dǎo)入的數(shù)據(jù)無(wú)論是文本還是 Excel都先轉(zhuǎn)成這三組矩陣再做后續(xù)計(jì)算不要試圖繞過(guò)這套數(shù)據(jù)結(jié)構(gòu)。3.2 跑起來(lái)牛頓法還是高斯賽德?tīng)栐趺催x以這套框架的常見(jiàn)命令風(fēng)格為例潮流求解的主函數(shù)只需要傳入三組矩陣外加兩個(gè)關(guān)鍵選項(xiàng)迭代算法和收斂容差。牛頓法收斂快、適合病態(tài)系統(tǒng)但對(duì)初值敏感高斯賽德?tīng)柗▽?shí)現(xiàn)簡(jiǎn)單、初值要求低但大系統(tǒng)里收斂慢。實(shí)際工程里我一般直接用牛頓法容差設(shè)在 1e-6最大迭代次數(shù)給 30 次。% 運(yùn)行潮流計(jì)算 [V, converged, iter] psflow(bus, line, gen, ... algorithm, newton, ... % 可選 newton 或 gauss tol, 1e-6, ... % 收斂判據(jù)電壓修正量范數(shù) maxiter, 30); % 迭代上限超過(guò)即判不收斂 % 檢查是否收斂 if ~converged error(潮流不收斂迭代 %d 次后仍未滿足精度, iter); end收斂判據(jù)這里的含義是相鄰兩次迭代的電壓修正向量范數(shù)小于 1e-6。設(shè)得越嚴(yán)計(jì)算越準(zhǔn)但對(duì)病態(tài)電網(wǎng)越容易觸發(fā)上限退出。如果發(fā)現(xiàn)不收斂與其把容差調(diào)松到 1e-4不如先檢查負(fù)荷和發(fā)電機(jī)出力是否嚴(yán)重不平衡——多數(shù)不收斂是初值離解太遠(yuǎn)不是算法不行。3.3 讀結(jié)果電壓、相角、線路潮流分別在哪潮流函數(shù)一般返回一個(gè)復(fù)數(shù)電壓向量V實(shí)部虛部分別對(duì)應(yīng)直角坐標(biāo)下的電壓分量或者用abs(V)和angle(V)轉(zhuǎn)成幅值與相角。要看線路有功和無(wú)功流向需要再用一組函數(shù)基于V和line矩陣回算線路功率。這里最容易踩坑的是V里元素的順序是按母線編號(hào)排的還是按 bus 矩陣的行號(hào)排的。兩者在多數(shù)算例里一致但如果你手動(dòng)刪過(guò) bus 某一行索引就會(huì)錯(cuò)位。% 提取幅值和相角 Vmag abs(V); % 各母線電壓幅值單位 pu Vang angle(V) * 180/pi; % 各母線相角單位 度 % 線路功率回算返回首末端有功/無(wú)功 [Pf, Qf, Pt, Qt] lineflow(bus, line, V); % 打印節(jié)點(diǎn)1的電壓和節(jié)點(diǎn)1-2線路有功 fprintf(V1 %.4f pu, theta1 %.2f deg\n, Vmag(1), Vang(1)); fprintf(P12 %.4f pu\n, Pf(1));線路功率回算是基于節(jié)點(diǎn)電壓和線路阻抗支路的復(fù)功率公式方向定義為首端流向末端。做校核時(shí)有個(gè)快速粗查辦法所有線路首端有功之和加上所有發(fā)電機(jī)出力必須等于所有負(fù)荷有功加全系統(tǒng)網(wǎng)絡(luò)損耗誤差超過(guò) 0.5% 就說(shuō)明某個(gè)矩陣填錯(cuò)了。4. 把潮流結(jié)果接進(jìn)暫態(tài)穩(wěn)定仿真三個(gè)必調(diào)參數(shù)4.1 仿真步長(zhǎng)與總時(shí)長(zhǎng)先定這兩個(gè)數(shù)暫態(tài)穩(wěn)定仿真的總時(shí)長(zhǎng)和步長(zhǎng)不是越大越好也不是越小越好?？倳r(shí)長(zhǎng)要覆蓋“第一擺”到“第二擺”的完整過(guò)程經(jīng)典機(jī)電暫態(tài)取 5 到 10 秒就夠了研究弱阻尼或低頻振蕩要放到 20 秒。步長(zhǎng)則取決于積分算法和你關(guān)心的最高頻率分量梯形隱式積分法在機(jī)電暫態(tài)尺度下一般取 0.01 秒極限試驗(yàn)時(shí)可以到 0.005 秒再小就是純浪費(fèi)時(shí)間。% 暫態(tài)穩(wěn)定求解常見(jiàn)命令風(fēng)格示例 [t, x, Vt] psstab(bus, line, gen, mac, ... tend, 10.0, ... % 仿真總時(shí)長(zhǎng)單位秒 dt, 0.01, ... % 固定步長(zhǎng)單位秒 event, event_seq); % 故障事件序列見(jiàn) 4.2仿真總時(shí)長(zhǎng)和事件時(shí)間有一個(gè)隱蔽關(guān)系如果事件發(fā)生時(shí)刻不在步長(zhǎng)網(wǎng)格上求解器要么就近取整、要么插值兩種處理都會(huì)讓波形產(chǎn)生虛假抖動(dòng)。我自己會(huì)先把所有事件時(shí)刻對(duì)齊到步長(zhǎng)的整數(shù)倍比如 1.00、1.10、1.30而不要寫(xiě) 1.05、1.12 這種步長(zhǎng)整除不了的數(shù)。4.2 故障序列接地、斷線、重合閘的事件表事件序列是暫態(tài)仿真的核心輸入通常用一張表來(lái)定義每一行是一個(gè)事件列包含觸發(fā)時(shí)間、動(dòng)作類(lèi)型、作用母線。以常見(jiàn)實(shí)現(xiàn)為例三相短路接地、故障清除、線路切除、重合閘都是通過(guò)修改導(dǎo)納矩陣來(lái)模擬的動(dòng)作的物理含義是“把某條母線對(duì)地導(dǎo)納改變到多少”。% 事件序列t1.0s 母線3三相短路1.1s清除1.3s重合閘 event_seq [ 1.00 3 3 -1 0.0 1; % t1.0母線3經(jīng)0.0歐接地三相短路 1.10 0 0 -1 0.0 0; % t1.1清除故障恢復(fù)正常拓?fù)?1.30 3 3 -1 0.0 1; % t1.3斷路器重合 ];事件表的時(shí)間列必須嚴(yán)格單調(diào)遞增同一個(gè)時(shí)刻只能有一個(gè)事件。想模擬單相接地需要把三相短路的“接地”改成“某相接地”并重新設(shè)置不對(duì)稱(chēng)系數(shù)多數(shù)機(jī)電暫態(tài)框架默認(rèn)只做正序網(wǎng)絡(luò)單相和三相短路的差異主要體現(xiàn)在故障期間的正序附加阻抗上不要過(guò)份追求不對(duì)稱(chēng)細(xì)節(jié)。4.3 動(dòng)態(tài)模型參數(shù)哪些是必調(diào)哪些可以默認(rèn)發(fā)電機(jī)動(dòng)態(tài)模型參數(shù)表中慣性常數(shù) H、暫態(tài)電抗 Xd、負(fù)序/次暫態(tài)參數(shù)、阻尼系數(shù) D 是最先要確認(rèn)的四個(gè)。經(jīng)典 2 階模型只需要 H、Xd 和 D4 階模型要補(bǔ)勵(lì)磁繞組時(shí)間常數(shù) Td06 階模型還要再加交軸阻尼繞組參數(shù)。勵(lì)磁系統(tǒng)的模型參數(shù)更敏感——調(diào)壓器放大倍數(shù)設(shè)大了暫態(tài)過(guò)程會(huì)高頻振鈴設(shè)小了電壓恢復(fù)慢。參數(shù)典型范圍調(diào)大后的現(xiàn)象調(diào)小后的現(xiàn)象慣性常數(shù) H (s)3.0 ~ 8.0轉(zhuǎn)子擺動(dòng)變慢第一擺幅度變小擺動(dòng)變快更容易失穩(wěn)阻尼系數(shù) D1.0 ~ 5.0振蕩衰減加快穩(wěn)態(tài)偏差變大振蕩衰減慢持續(xù)搖擺暫態(tài)電抗 Xd (pu)0.2 ~ 0.4極限切除時(shí)間變短更難穩(wěn)定穩(wěn)定裕度增大勵(lì)磁放大倍數(shù) K_A20 ~ 200電壓響應(yīng)快可能高頻振蕩電壓恢復(fù)慢動(dòng)態(tài)壓降大這些參數(shù)在真實(shí)發(fā)電機(jī)銘牌和試驗(yàn)報(bào)告里都有據(jù)可查做研究時(shí)建議從典型值起步不要一上來(lái)就調(diào)成極端值。一個(gè)很容易被忽視的細(xì)節(jié)是調(diào)速器模型參數(shù)中調(diào)差系數(shù)對(duì)潮流結(jié)果有影響但對(duì)暫態(tài)第一擺影響不大如果你只關(guān)心故障后 2 秒內(nèi)的行為調(diào)速器用默認(rèn)參數(shù)即可把精力花在勵(lì)磁和發(fā)電機(jī)本體的參數(shù)上。5. 避坑指南暫態(tài)穩(wěn)定仿真最容易翻車(chē)的 4 個(gè)數(shù)據(jù)坑5.1 潮流不收斂但程序不報(bào)錯(cuò)最后全輸出 1.0 pu現(xiàn)象潮流函數(shù)返回converged true但所有母線電壓都是 1.0相角全是 0和物理常識(shí)完全不符。原因典型的“初值等于解”的假收斂。當(dāng)負(fù)荷和發(fā)電機(jī)出力都被填成 0或者系統(tǒng)基準(zhǔn) MVA 和發(fā)電機(jī)數(shù)據(jù)基準(zhǔn) MVA 不一致時(shí)比如發(fā)電機(jī)數(shù)據(jù)按 200MVA 基準(zhǔn)填寫(xiě)系統(tǒng)基準(zhǔn)卻是 100MVA迭代方程在原點(diǎn)附近自然收斂但這個(gè)解沒(méi)有承載任何功率。解決先用功率平衡粗查。把全系統(tǒng) P_load 加起來(lái)減去所有 Pg差值應(yīng)該在網(wǎng)損量級(jí)如果差值幾乎為 0大概率是基準(zhǔn)值錯(cuò)了。把發(fā)電機(jī)有功整體乘 0.5 再跑一次如果電壓和相角有明顯變化說(shuō)明數(shù)據(jù)通路是通的只是基準(zhǔn)需要統(tǒng)一。5.2 暫態(tài)波形在故障清除時(shí)刻出現(xiàn)鋸齒形突變現(xiàn)象搖擺曲線在 1.1 秒附近出現(xiàn)一個(gè)很窄的脈沖或凹坑看起來(lái)像數(shù)值振蕩但步長(zhǎng)已經(jīng)很小了。原因事件發(fā)生時(shí)間和步長(zhǎng)網(wǎng)格不重合。比如步長(zhǎng) 0.01 秒事件寫(xiě)的是 1.105 秒積分器會(huì)就近取整到 1.10 或 1.11再?gòu)?qiáng)行觸發(fā)代數(shù)變量跳變形成的波形就帶細(xì)鋸齒。這在重合閘時(shí)刻尤其明顯因?yàn)閿嗦菲鲃?dòng)作本身就會(huì)造成一個(gè)不連續(xù)的拓?fù)淝袚Q。解決把每個(gè)事件時(shí)刻都對(duì)齊到步長(zhǎng)的整數(shù)倍。我通常先定步長(zhǎng)再反推事件時(shí)間步長(zhǎng) 0.01 秒事件就寫(xiě) 1.00、1.10、1.30需要更精細(xì)的故障時(shí)刻就把步長(zhǎng)改成 0.005 秒讓 1.105 變得合法。另一個(gè)辦法是檢查求解器有沒(méi)有事件插值開(kāi)關(guān)開(kāi)了反而會(huì)引入新的數(shù)值阻尼慎用。5.3 線路電抗為負(fù)功率方向全亂現(xiàn)象潮流算出來(lái)某條線路的有功方向是從負(fù)荷側(cè)流向電源側(cè)且該線路的 X 列是負(fù)數(shù)。原因兩種可能。一種是 X 列本來(lái)填的就是串補(bǔ)電容的容抗為負(fù)正常另一種更常見(jiàn)——R 和 X 填反了或者 bus/line 矩陣?yán)锏牧许樞蚝臀臋n對(duì)不上導(dǎo)致程序把某列當(dāng)成了 X。這類(lèi)錯(cuò)誤最陰險(xiǎn)因?yàn)槌绷髡諛邮諗恐皇墙Y(jié)果不物理。解決先把該線路單獨(dú)拿出來(lái)算一遍手算值已知兩端電壓幅值和相角差手算線路有功傳輸公式 P V1*V2/X * sin(theta1 - theta2)。如果手算是正方向而程序結(jié)果是負(fù)方向就是矩陣列錯(cuò)位。我會(huì)把 bus 和 line 矩陣的表頭注釋永遠(yuǎn)留在腳本第一段不刪除、不壓縮切換數(shù)據(jù)源時(shí)先拿一個(gè)已知算例驗(yàn)證列順序。5.4 動(dòng)態(tài)仿真初始功率和潮流結(jié)果對(duì)不上現(xiàn)象潮流計(jì)算顯示發(fā)電機(jī)出力是 0.8 pu但暫態(tài)仿真的第一幀結(jié)果里同一條母線的注入功率變成了 0.75 pu二者對(duì)不上。原因動(dòng)態(tài)模型初始化沒(méi)有正確吃進(jìn)潮流結(jié)果。發(fā)電機(jī)的內(nèi)電勢(shì)、勵(lì)磁電壓初值、調(diào)速器閥門(mén)開(kāi)度都需要從潮流結(jié)果反推常見(jiàn)實(shí)現(xiàn)里這一步叫“初值計(jì)算”或“初始化”。如果動(dòng)態(tài)參數(shù)表里的 Xd 與潮流計(jì)算用的 Xd 不一致反推出來(lái)的內(nèi)電勢(shì)就會(huì)帶誤差進(jìn)而影響第一幀功率。解決檢查初始化函數(shù)是否基于同一份潮流結(jié)果。我的習(xí)慣是先跑潮流把 V 作為初始化函數(shù)的輸入不要初始化和潮流各算各的。另一個(gè)自查辦法把暫態(tài)仿真的前 0.1 秒結(jié)果取平均正常情況下它應(yīng)該非常接近潮流結(jié)果如果差得明顯優(yōu)先懷疑動(dòng)態(tài)電抗參數(shù)不一致而不是懷疑求解器。調(diào)試時(shí)把這個(gè)差值打印出來(lái)比盯著整條曲線找問(wèn)題快得多。6. 進(jìn)階用法做出自己的算例并用結(jié)果反向驗(yàn)證6.1 把自己的數(shù)據(jù)寫(xiě)成算例文件如果你拿到的是某電網(wǎng)或某設(shè)備廠商的導(dǎo)出數(shù)據(jù)常見(jiàn)格式是帶表頭的 CSV 或 Excel。不要手動(dòng)復(fù)制進(jìn)矩陣——用腳本自動(dòng)拼。先讀表頭確認(rèn)列名再按工具箱要求的列順序重新組織輸出成caseN.m腳本。關(guān)鍵步驟是給每列做一次單位折算發(fā)電機(jī)容量基準(zhǔn)、電壓基準(zhǔn)、阻抗基準(zhǔn)三者必須一致。我會(huì)習(xí)慣在數(shù)據(jù)文件頭部用注釋寫(xiě)明“基準(zhǔn)容量 100MVA基準(zhǔn)電壓按各區(qū)平均電壓取”防止三個(gè)月后回來(lái)看不懂。% 從CSV自動(dòng)生成bus矩陣的示例邏輯 csv readtable(network.csv); bus zeros(height(csv), 10); bus(:,1) csv.bus_number; bus(:,3) (csv.load_kW / 100000) * exp(0); % kW - 標(biāo)幺P bus(:,4) (csv.load_kvar / 100000); % kvar - 標(biāo)幺Q % 列2電壓初值默認(rèn)1.0列9/10為電壓限值默認(rèn)1.1/0.96.2 三個(gè)自查手段功率平衡、電壓越限、曲線單調(diào)性做完仿真先別畫(huà)圖用三個(gè)腳本內(nèi)聯(lián)檢查代替肉眼全網(wǎng)功率平衡檢查、電壓越限掃描、搖擺曲線局部單調(diào)性檢查。功率平衡直接決定潮流數(shù)據(jù)可信度電壓越限掃描幫你發(fā)現(xiàn)會(huì)不會(huì)有接近失穩(wěn)的邊界曲線單調(diào)性則是數(shù)值振蕩的信號(hào)——真實(shí)機(jī)電暫態(tài)不會(huì)出現(xiàn)持續(xù)等幅鋸齒出現(xiàn)就說(shuō)明步長(zhǎng)或初值有問(wèn)題。% 電壓越限掃描 Vfinal abs(Vt(end,:)); if any(Vfinal 1.1) || any(Vfinal 0.9) warning(存在電壓越限節(jié)點(diǎn)max%.3f pu, min%.3f pu, ... max(Vfinal), min(Vfinal)); end % 第一擺單調(diào)性轉(zhuǎn)子角前0.5s內(nèi)不應(yīng)出現(xiàn)頻繁反向 theta1 angle(Vt(:,1)); % 取1號(hào)機(jī)轉(zhuǎn)子角隨時(shí)間序列 dtheta diff(theta1(1:round(0.5/t))); if sum(dtheta 0) length(dtheta) * 0.8 warning(轉(zhuǎn)子角前0.5s內(nèi)有劇烈反向疑似數(shù)值振蕩); end我最早做這類(lèi)仿真的時(shí)候習(xí)慣拿到結(jié)果直接畫(huà)曲線眼睛看著差不多就完事。后來(lái)在一次故障時(shí)間參數(shù)微調(diào)中吃了虧——曲線整體變好看了但功率平衡檢查差了 3%追查半天才發(fā)現(xiàn)是某條線路的 R 填錯(cuò)了一位小數(shù)。從那以后我把這三個(gè)檢查寫(xiě)進(jìn)每次仿真的收尾腳本里不改數(shù)據(jù)、不重跑結(jié)果不合格就拒絕看曲線。這算不上什么聰明技巧但確實(shí)讓我的仿真結(jié)果可信度上了一個(gè)臺(tái)階。希望這篇筆記能幫你少走這段彎路把精力花在真正要研究的控制策略和穩(wěn)定性問(wèn)題上。本文還有配套的精品資源點(diǎn)擊獲取