【正文】 ubplot(3,1,1)。plot(t,y(:,1))。% 下面指定縱橫軸標(biāo)簽xlabel(39。Time[s]39。)ylabel(39。Id[A]39。)title(39。有阻尼三相突然短路特性39。)hold onsubplot(3,1,2)。plot(t,y(:,3))。% 下面指定縱橫軸標(biāo)簽xlabel(39。Time[s]39。)ylabel(39。Ifd[A]39。)hold onw=314。ia=y(:,1).*sin(w*t)+y(:,2).*cos(w*t)。subplot(3,1,3)。plot(t,ia)% 下面指定縱橫軸標(biāo)簽xlabel(39。Time[s]39。)ylabel(39。Ia[A]39。)三相同步發(fā)電機(jī)三相突然短路及轉(zhuǎn)子繞組M文件：% 同步發(fā)電機(jī)有阻尼繞組轉(zhuǎn)子繞組短路,定子繞組突然加上對稱電壓% 將該M函數(shù)定義為zzrzdlfunction dydt=zzrzdl(t,y)% 下面輸入電機(jī)基本數(shù)據(jù)：r=。Rfd=。Rkd=。Rkq=。Ufd=24。w=377。Ll=。Lmd=。Lmq=。Llfd=。Llkd=。Llkq=。Ld=Lmd+Ll。Lq=Lmq+Ll。Mafd0=Lmd。Makd0=Lmd。Makq0=Lmq。Lfd=Llfd+Lmd。Lkd=Llkd+Lmd。Lkq=Llkq+Lmq。Mfkd=Lmd。% 下面輸入電感系數(shù)矩陣：L=[ Ld, 0, Lmd, Lmd, 0。 0, Lq, 0, 0, Lmq。 Lmd, 0, Lfd, Lmd, 0。 Lmd, 0, Lmd, Lkd, 0。 0, Lmq, 0, 0, Lkq] 。 G=[ 0, Lq, 0, 0, 3/2*Makq0。 Ld, 0, Mafd0, Makd0, 0。 0, 0, 0, 0, 0。 0, 0, 0, 0, 0。 0, 0, 0, 0, 0]。% 下面輸入電阻矩陣：R=[ r, 0, 0, 0, 0。 0, r, 0, 0, 0。 0, 0, Rfd, 0, 0。 0, 0, 0, Rkd, 0。 0, 0, 0, 0, Rkq]。% 下面輸入電壓向量：Udq0=[100,0,0,0,0]39。% 下面列寫微分方程：dydt=L\(Udq0w*G*yR*y)。% 編寫求解空載建立電壓微分方程的函數(shù)(zzrzdlfz)有阻尼轉(zhuǎn)子繞組短路定子繞組突然加對稱電壓function zzrzdlfz %設(shè)定時間和初始值tspan = [0 10]。y0 = [0。0。0。0。0]。% 方程求解[t,y]=ode113(39。zzrzdl39。,tspan,y0)。%繪制定子電流和勵磁電流變化規(guī)律subplot(3,1,1)。plot(t,y(:,1))xlabel(39。Time[s]39。)ylabel(39。Id[A]39。)title(39。有阻尼轉(zhuǎn)子繞組短路定子繞組突然加對稱電壓特性39。)hold onsubplot(3,1,2)。plot(t,y(:,3))。xlabel(39。Time[s]39。)ylabel(39。Ifd[A]39。)hold onsubplot(3,1,3)。plot(t,y(:,4))xlabel(39。Time[s]39。)ylabel(39。Ikd[A]39。)無阻尼篇：三相同步發(fā)電機(jī)空載特性M文件：% 編寫同步發(fā)電機(jī)無阻尼繞組空載建立電壓過程微分方程的M函數(shù)% 將該M函數(shù)定義為wzlkzfunction dydt=wzlkz(t,y)% 下面輸入電機(jī)基本數(shù)據(jù)：r=。Rfd=5。9013E01。Rkd=。Rkq=。Ufd=24。w=377。Ll=。Lmd=。Lmq=。Llfd=。Llkd=。Llkq=。Ld=Lmd+Ll。Lq=Lmq+Ll。Mafd0=Lmd。Makd0=Lmd。Makq0=Lmq。Lfd=Llfd+Lmd。Lkd=Llkd+Lmd。Lkq=Llkq+Lmq。Mfkd=Lmd。% 下面輸入電感系數(shù)矩陣：L=[ Ld, 0, Lmd。 0, Lq, 0。 Lmd, 0, Lfd ]。 G=[ 0, Lq, 0 。 Ld, 0, Mafd0。 0, 0, 0 ]。 % 下面輸入電阻矩陣：R=[r,0,0。 0,r,0。0,0,Rfd]。% 下面輸入電壓向量：Uabc=[0,0,Ufd]39。% 下面列寫微分方程：dydt=L\(Uabcw*G*yR*y)。% 編寫求解空載建立電壓微分方程的函數(shù)(zlkzfz)無阻尼空載特性曲線function wzlkzfz %設(shè)定時間和初始值tspan = [0 10]。y0 = [0。0。0]。% 方程求解[t,y]=ode113(39。wzlkz39。,tspan,y0)。%空載建立電壓過程中勵磁電流的變化規(guī)律。ifd=y(:,3)。plot(t,ifd)% 下面指定縱橫軸標(biāo)簽xlabel(39。Time[s]39。) ylabel(39。Ifd[A]39。) title(39。無阻尼空載特性39。)三相同步發(fā)電機(jī)三相突然短路特性M文件：% 編寫同步發(fā)電機(jī)無阻尼繞組三相突然短路微分方程的M函數(shù)% 將該M函數(shù)定義為wsanxtrdlfunction dydt=wsanxtrdl(t,y)% 下面輸入電機(jī)基本數(shù)據(jù)：r=。Rfd=。Rkd=。Rkq=。Ufd=24。w=377。Ll=。Lmd=。Lmq=。Llfd=。Llkd=。Llkq=。Ld=Lmd+Ll。Lq=Lmq+Ll。Mafd0=Lmd。Makd0=Lmd。Makq0=Lmq。Lfd=Llfd+Lmd。Lkd=Llkd+Lmd。Lkq=Llkq+Lmq。Mfkd=Lmd。% 下面輸入電感系數(shù)矩陣：L=[ Ld, 0, Lmd, 。 0, Lq, 0, 。 Lmd, 0, Lfd ]。 G=[ 0, Lq, 0, 。 Ld, 0, Mafd0, 。 0, 0, 0, ]。 % 下面輸入電阻矩陣：R=[ r, 0, 0, 。 0, r, 0, 。 0, 0, Rfd, ]。% 下面輸入電壓向量：Udq0=[100,10,24]39。% 下面列寫微分方程：dydt=L\(Udq0w*G*yR*y)。% 編寫求解空載建立電壓微分方程的函數(shù)(wsanxtrdlfz)無阻尼三相突然短路function wsanxtrdlfz %設(shè)定時間和初始值tspan = [0 10]。y0 = [0。0。10]。% 方程求解[t,y]=ode23tb(39。wsanxtrdl39。,tspan,y0)。%繪制定子電流和勵磁電流變化規(guī)律subplot(3,1,1)。plot(t,y(:,3))。xlabel(39。Time[s]39。)ylabel(39。Id[A]39。)title(39。無阻尼三相突然短路特性39。)hold onw=377。ia=y(:,1).*sin(w*t)+y(:,2).*cos(w*t)。subplot(3,1,2)。plot(t,ia) xlabel(39。Time[s]39。) ylabel(39。Ia[A]39。)hold onsubplot(3,1,3)。plot(t,y(:,3))。 xlabel(39。Time[s]39。) ylabel(39。Ifd[A]39。)三相同步發(fā)電機(jī)三相突然短路及轉(zhuǎn)子繞組M文件：% 同步發(fā)電機(jī)無阻尼繞組轉(zhuǎn)子繞組短路,定子繞組突然加上對稱電壓% 將該M函數(shù)定義為wzzrzdlfunction dydt=wzzrzdl(t,y)% 下面輸入電機(jī)基本數(shù)據(jù)：r=。Rfd=。Rkd=。Rkq=。Ufd=24。w=377。Ll=。Lmd=。Lmq=。Llfd=。Llkd=。Llkq=。Ld=Lmd+Ll。Lq=Lmq+Ll。Mafd0=Lmd。Makd0=Lmd。Makq0=Lmq。Lfd=Llfd+Lmd。Lkd=Llkd+Lmd。Lkq=Llkq+Lmq。Mfkd=Lmd。% 下面輸入電感系數(shù)矩陣：L=[ Ld, 0, Lmd, 。 0, Lq, 0, 。 Lmd, 0, Lfd ]。 G=[ 0, Lq, 0, 。 Ld, 0, Mafd0, 。 0, 0, 0, ]。 % 下面輸入電阻矩陣：R=[ r, 0, 0, 。 0, r, 0, 。 0, 0, Rfd, ]。% 下面輸入電壓向量：Udq0=[100,0,0]39。% 下面列寫微分方程：dydt=L\(Udq0w*G*yR*y)。% 編寫求解空載建立電壓微分方程的函數(shù)(wzzrzdlfz)無阻尼轉(zhuǎn)子繞組短路定子繞組突然加對稱電壓function wzzrzdlfz%設(shè)定時間和初始值tspan = [0 10]。y0 = [0。0。0]。% 方程求解[t,y]=ode113(39。wzzrzdl39。,tspan,y0)。%繪制定子電流和勵磁電流變化規(guī)律subplot(2,1,1)plot(t,y(:,3))xlabel(39。Time[s]39。)ylabel(39。Ifd[A]39。)title(39。無阻尼轉(zhuǎn)子繞組短路定子繞組突然加對稱電壓特性39。)hold onsubplot(2,1,2)w=377。ia=y(:,1).*sin(w*t)+y(:,2).*cos(w*t)。plot(t,ia) xlabel(39。Time[s]39。) ylabel(39。Ia[A]39。)