电机设计matlab程序
%电机设计程序
clear all
format short e
m1=3;p=2;f=50
%1.额定功率
PN=5.5*10^3 ;
%2.额定电压（单位V，三角形接法）
UN=380;UN0=380;
%3.功电流(单位A)
IKW=PN/(m1*UN0)
%4.效率eta按照技术条件的规定eta=0.875
eta=0.855 ;
%5.功率因数cos(phi) =0.84,按照技术条件的规定cos(phi)=0.84
phi=acos(0.84);
cos(phi);
%6.极对数p=2
p=2;
%7.定转子槽数：每极每相槽数取整数。参考类似规格电机取q1=3,则Z1=2m1pq1，再查表10-8选Z2=32，并采用转子斜槽。
q1=3;
Z1=2*m1*p*q1
Z2=32 ;
%8.定转子每极槽数
Zp1=Z1/(2*p)
Zp2=Z2/(2*p)
%9.确定电机的主要尺寸；一般可参考类似电机的主要尺寸来确定Di1和lef.现按10-2中的
KE1=0.0108*log(PN/1000)-0.013*p+0.931
P1=KE1*PN/(eta*cos(phi))
alphap1=0.68;KNm1=1.10;Kdp1=0.96;A1=25000;
Bdelta1=0.69;n1=1450;
V=(6.1/(alphap1*KNm1*Kdp1))*(1/(A1*Bdelta1 ))*(P1/n1)
D1=0.21;
%铁心的有效长度
Di1=0.136;
lef =V/((Di1)^2)
%气隙的确定
%参考类似产品或由经验公式（10-10a），得
lt=0.115;
delta =0.0004
lef=lt + 2*delta
D2=Di1-2*delta
%转子内径先按转轴直径决定（以后再校验转子轭部磁密）
Di2=0.048 ;
%11.极距tau
tau =pi*Di1/(2*p)
%12.定子齿距t1
t1=(pi*Di1/Z1)
%转子齿距t2
t2=(pi*D2/Z2)
bsk=0.01187;
%15.设计定子绕组
Nphi11=eta*cos(phi)*pi*Di1*A1/(m1*IKW)
%取并联支路a1=1,由式（10-15），可得每槽导体数
a1=1;
Ns1=47
%16.每相串联导体数Nphi1
Nphi1=Ns1*Z1/(m1*a1)
%每相串联匝数N1
N1=Nphi1/2
%17.绕组线规设计
%初选定子电密J11=5.0A/mm^2,由式（10-16），计算导线并绕根数和每根导线面积的乘积。
J11=5.0;
%其中定子电流初步估计值
I11=IKW/(eta*cos(phi))
Nt1Ac11=I11/(a1*J11)
%在附录二中选用截面积相近的铜线，高强度漆包线
Nt1=2;d1=0.00090;d=0.00095;Ac11=0.0000006362;Nt1Ac11=0.0000006362*2;
%18. 设定子槽形
%因定子绕组为原导线散嵌，故采用梨形槽，齿部平行。初步取Bt11=1.4T,按式（10-18）估计定子齿宽
Bt11=1.493641;
KFe=0.95 ;
bt1=t1*Bdelta1/(KFe*Bt11)
%初步取Bj11=1.25T,按式（10-19），估计定子轭部计算高度
Bj11=1.25;
hj11=tau*alphap1*Bdelta1/(2*KFe*Bj11)
%按齿宽和定子轭部计算高度的估算值作出定子槽形如图10-24，槽口尺寸参考类似产品决定，齿宽计算如下：
b01=0.0035;
b11=0.0067;
h01=0.0008;
h11=(b11-b01)/2*tan(pi/6);
h21=0.0145-h11
r21=0.0044;
bi11=pi*(Di1+2*h01+2*h21+2*h11)/Z1-2*r21;
bi12=pi*(Di1+2*h01+2*h11)/Z1-b11;
bt1=0.5*(bi12+bi11)
%齿部基本平行，齿宽
%19.槽满率
%曹面积
hs1=h11+h21;h=0.002 ;
As=((2*r21+b11)/2)*(hs1-h)+pi*(r21)^2/2
%按附录三，槽绝缘采用DMD复合绝缘，,槽楔为h=2mm层压板，槽绝缘占面积Ai
deltai=0.00023;
Ai=deltai*(2*hs1+pi*r21)
%槽的有效面积
Aef=As-Ai
%槽满率
sf=Nt1*Ns1*d^2/Aef
%20.绕组系数
alpha1=(p*360/Z1)
Kp1=1;alpha=alpha1*pi/180,q=q1;
Kd1=sin(q*alpha/2)/(q*sin(alpha/2))
Kdp1=Kd1*Kp1
%每相有效串联导体数Nef
Nef=Nphi1*Kdp1
%21.设计转子槽形与转子绕组
%按式（10-39），预计转子导条电流：
Ki=0.89 ;
I21=Ki*I11*(3*Nphi1*Kdp1)/Z2
%初步取转子导条电密JB1=3.5A/mm^2,于是导条截面积
JB1=3.5;
AB1=I21/JB1
%初步取Bt21=1.3T,估算转子齿宽
Bt21=1.3;
bt21=(t2/10*Bdelta1)/(KFe*Bt21)
%初步取Bj21=1.25T,估计转子轭部计算高度
Bj21=1.25;
hj21=tau*alphap1*Bdelta1/(2*KFe*Bj21)
%为获得比较好的起动性能，采用平行槽，作槽形图如图10-25所示，取槽口尺寸b02=1mm,h02=0.5mm
b02=0.001;h02=0.0005 ;
%齿壁不平行的槽形的宽度（按3-3）计算如下：
b12=0.0055;h12=0.0011;h22=0.026-0.0011;
h12=(b12-b02)/2*tan(pi/6);h22=0.023-h12;b22=0.003;
a=0.026;b=0.010;
bt213=pi*(D2-2*2/3*(h02+h12+h22))/Z2-(b22+(b12-b22)/3)
%导条截面积（转子槽面积）
AB=(b02+b12)/2*h12+(b12+b22)/2*h22
%按式（11-41）估计端环电流
IR1=I21*Z2/(2*p*pi)
%端环所需面积
JR1=0.6*JB1;
AR11=IR1/JB1
AR11=(a+b)/2*(h02+h12+h22)*10^6
AR1=input('输入端环面积AR1=')
%其中端环电密JR1=0.6,JB1=2.1A/mm^2.按照工艺要求有所需面积确定端环内径及厚度如图10-25b,得端环面积AE=40*10^(-6)
%(二)磁路计算
%22
KE1=0.918;epsilonl1=1-0.918;
E1=(1-epsilonl1)*UN
%23
Ks(1)=1.2;KNm(1)=1.095;
precision=1;
while precision >0.01
for n=1:3
KNm(n)=input('输入波形系数KNm(n)=')
alphap(n)=input('输入极弧系数alphap(n)=')
phi(n)=E1/(4*KNm(n)*Kdp1*f*N1)
%24
At1=KFe*lt*bt1*Zp1
At2=KFe*lt*bt213*Zp2
%25
hs1=h21+h01+h11+r21;
hs2=h02+h12+h22;
hj11=(D1-Di1)/2-hs1+r21/3
hj21=(D2-Di2)/2-hs2
Aj1=KFe*lt*hj11
Aj2=KFe*lt*hj21
%26
Adelta=tau*lef
%27
Fs(n)=input('输入Fs(n)=')
%28
Bdelta(n)=Fs(n)*phi(n)/Adelta
%29
Bt1=Fs(n)*phi(n)/At1
%30
Bt2=Fs(n)*phi(n)/At2
%31
Ht1(n)=input('输入Ht1(n)=');Ht2(n)=input('输入Ht2(n)=');
%32
Kdelta1=t1*100*(4.4*delta*100+0.75*b01*100)/(t1*100*(4.4*delta*100+0.75*b01* 100)-(b01*100)^2)
Kdelta2=t2*100*(4.4*delta*100+0.75*b02*100)/(t2*100*(4.4*delta*100+0.75*b02* 100)-(b02*100)^2)
Kdelta=Kdelta1*Kdelta2
deltaef=Kdelta*delta
%33
Lt1=(h11+h21)+1/3*r21
Lt2=(h12+h22)
%34
Lj11=pi*(D1-hj11)/(4*p)
Lj21=pi*(Di2+hj21)/(4*p)
%35
mu0=4*pi*10^(-7);
Fdelta(n)=Kdelta*Bdelta(n)*delta/mu0
%36
Ft1(n)=Ht1(n)*Lt1*100
Ft2(n)=Ht2(n)*Lt2*100
%37
Ks(2*n)=(Fdelta(n)+Ft1(n)+Ft2(n))/Fdelta(n)
Ks(2*n+1)=Ks(2*n)-(Ks(2*n)-Ks(2*n-1))/3
precision=abs((Ks(2*n-1)-Ks(2*n))/Ks(2*n))
end
end
Ks=Ks(6);alphap1=alphap(3);phi=phi(3);Bdelta=Bdelta(3);Fdelta=Fdelta(3);Ft1=Ft1( 3);Ft2=Ft2(3);
KNm=KNm(3);Fs=Fs(3);Ht1=Ht1(3);Ht2=Ht2(3);
%迭代的出来的结论
%38
Bj1=0.5*phi/Aj1
%39
Bj2=0.5*phi/Aj2
%40
Hj1=input('输入Hj1=')
Hj2=input('输入Hj2=')
%41
hj11/tau
Cj1=input('根据Bj1和hj11/tau输入Cj1=')
Fj1=Cj1*Hj1*Lj11*10^2
hj21/tau
Cj2=input('根据Bj2和hj21/tau输入Cj2=')
Fj2=Cj2*Hj2*Lj21*10^2
%42
F0=Fdelta+Ft1+Ft2+Fj1+Fj2
%43
Im=2*p*F0/(0.9*m1*N1*Kdp1)
%44
Ims=Im/IKW
%45
Xms=(4*f* mu0*m1/pi)*((N1*Kdp1)^2/(Ks*p))*lef*(tau/deltaef)
Xmss=Xms*IKW/UN
%（三)参数计算
%46
beta=(2*8+1*7)/(3*9)
tauy=pi*(Di1+2*(h01+h11)+h21+r21)*beta/(2*p)
d1=0.015;
lB=lt+2*d1
K0=1.2;
l0=lB+K0*tauy
%47
lE=2*d1+K0*tauy
%48
Cx=4*pi*f*mu0*(N1*Kdp1)^2*lef*PN/(m1*p*UN^2)
%49
KU1=1.0;KL1=1.0;
lambdaU1=h01/b01+2*h11/(b01+b11)
h212r21=h21/(2*r21),b112r21=b11/(2*r21)
lambdaL1=input('输入lambdaL1=')
lambdas1=KU1*lambdaU1+KL1*lambdaL1
%50
Xs1s=2*m1*p/(Z1*Kdp1^2)*(lt/lef)*lambdas1*Cx
%51
sigmas=0.0129;
Xdelta1s=(m1/(pi^2)*(tau/deltaef)*(sigmas/((Kdp1^2)*Ks)))*Cx %52
XE1s=0.47*(lE-0.64*tauy)*Cx/(lef*Kdp1^2)
%53
Xsigma1s=Xs1s+Xdelta1s+XE1s
%54
%55
lambdaU2=h02/b02
hb2=h22/b22,b1b2=b12/b22
lambdaL=input('输入lambdaL=');
lambdaL2=2*h12/(b02+b12)+lambdaL
lambdas2=lambdaU2+lambdaL2
%56
Xs2s=2*m1*p/Z2*(lt/lef)*lambdas2
%57
Z22p=Z2/(2*p)
sigmaR=input('根据Z22p输入sigmaR=')
Xdelta2s=m1*tau*sigmaR/(pi^2*deltaef*Ks)
%58
DR=0.107;
XE2s=(0.2523*Z2*DR/(2*p*lef*2*p))*(2*m1*p/Z2)
%59
Xsks=0.5*(bsk/t2)^2*Xdelta2s
%60
Xsigma2s=(Xs2s+Xdelta2s+Xsks+XE2s)*Cx
%61
Xsigmas=Xsigma1s+Xsigma2s
%62
rhow=0.0217*10^(-6);
R1=rhow*(2*N1*l0)/(Nt1*Ac11*a1)
%63
R1s=R1*IKW/UN
%64
C=1.05;rhoCu=8.9*10^3;
GCu=C*l0*Ns1*Z1*Ac11*Nt1*rhoCu
rhoFe=7.8*10^3;
GFe=KFe*lt*(D1+delta)^2*rhoFe
%65 AR1代替AR
KB=1.04;rhow1=0.0434*10^(-6);
R21=rhow1*(KB*lt/AB+Z2*DR/(2*pi*p^2*AR1*10^(-6)))*(4*m1*(N1*Kdp1)^2/Z 2)
RB1=rhow1*(KB*lt/AB)*(4*m1*(N1*Kdp1)^2/Z2)
RBs=RB1*IKW/UN
RR1=rhow1*Z2*DR/(2*pi*p^2*AR1*10^(-6))*(4*m1*(N1*Kdp1)^2/Z2)
RRs=RR1*IKW/UN
R2s=R21*IKW/UN
%66
eta(1)=0.855;
I1Ps(1)=1/eta(1);
precisioneta=1;
while precisioneta>0.005
for m=1:2
I1Ps(m)=1/eta(2*m-1)
%67
sigma1=1+Xsigma1s/Xmss
IXs(m)=sigma1*Xsigmas*(I1Ps(m))^2*(1+(sigma1*Xsigmas*I1Ps(m))^2)
%68
I1Qs(m)=Ims+IXs(m)
%69
epsilonL(m)=I1Ps(m)*R1s+I1Qs(m)*Xsigma1s
KE=1-(I1Ps(m)*R1s+I1Qs(m)*Xsigma1s)
%70
epsilon0=Ims*Xsigma1s
1-epsilon0
%71
Bt10(m)=(1-epsilon0)/(1-epsilonL(m))*Bt1
Ht10(m)=input('输入Ht10(m)=')
Bt20(m)=(1-epsilon0)/(1-epsilonL(m))*Bt2
Ht20(m)=input('输入Ht20(m)=')
%73
Bj10(m)=(1-epsilon0)/(1-epsilonL(m))*Bj1
Hj10(m)=input('输入Hj10(m)=')
%74
Bj20(m)=(1-epsilon0)/(1-epsilonL(m))*Bj2
Hj20(m)=input('输入Hj20(m)=')
%75
Bdelta0(m)=(1-epsilon0)/(1-epsilonL(m))*Bdelta
%76
Ft10(m)=Ht10(m)*Lt1*100
%77
Ft20(m)=Ht20(m)*Lt2*100
%78
hj11/tau,Bj10
Cj1(m)=input('输入Cj1(m)=')
Fj10(m)=Cj1(m)*Hj10(m)*Lj11*100
%79
hj21/tau,Bj20
Cj2(m)=input('输入Cj2(m)=')
Fj20(m)=Cj2(m)*Hj20(m)*Lj21*100
%80
Fdelta0(m)=Kdelta*delta*Bdelta0(m)/mu0
%81
F00(m)=Fdelta0(m)+Fj20(m)+Fj10(m)+Ft20(m)+Ft10(m) %82
Im0=2*p*F00(m)/(0.9*m1*N1*Kdp1)
%(四)工作性能计算
%83
I1s(m)=(I1Ps(m)^2+I1Qs(m)^2)^(0.5)
I1(m)=I1s(m)*IKW
%84
J1(m)=(I1(m)/(a1*Nt1*Ac11))/10^6
%85
A1(m)=m1*Nphi1*I1(m)/(pi*Di1)
%86
I2s(m)=(I1Ps(m)^2+IXs(m)^2)^(0.5)
I2(m)=I2s(m)*IKW*m1*Nphi1*Kdp1/Z2
IR(m)=I2(m)*Z2/(2*pi*p)
%87
JB(m)=I2(m)/(AB*10^6)
JR(m)=IR(m)/(AR1)
%88
pCu1s(m)=I1s(m)^2*R1s
pCu1(m)=pCu1s(m)*PN
%89
pA12s(m)=I2s(m)^2*R2s
pA12(m)=pA12s(m)*PN
%90
pss=0.02;
ps=pss*PN
%91
pfw=(3/p)^2*(D1)^4*10^4
pfws=pfw/PN
%92
Bj10
phej(m)=input('输入phej(m)=')
k2=2;
Gj=4*p*Aj1*Lj11*rhoFe
pFej(m)=k2*phej(m)*Gj
Bt10
phet(m)=input('输入phet(m)=')
Gt=2*p*At1*Lt1*rhoFe
k1=2;
pFet(m)=k1*phet(m)*Gt
pFe(m)=pFej(m)+pFet(m)
pFes(m)=pFe(m)/PN
%93
sigmaps(m)=pCu1s(m)+pA12s(m)+pss+pfws+pFes(m) %94
PN1s(m)=1+sigmaps(m)
%95
eta(2*m)=1-sigmaps(m)/PN1s(m)
eta(2*m+1)=eta(2*m)+(eta(2*m)-eta(2*m-1))/5
precisioneta=abs((eta(2*m)-eta(2*m-1))/eta(2*m))
end
end
%96
I1s=I1s(2);I1Ps=I1Ps(2);
phi=acos(I1Ps/I1s);
cos(phi)
%97
pFers=((1-1/2)*pFej+(1-1/2.5)*pFet)*10^(-3)/PN*1000
sN=pA12s/(1+pA12s+pFers+pss+pfws)
%98
nN=60*f/p*(1-sN)
%99
Tms=(1-sN)/(2*(R1s+(R1s^2+Xsigmas^2)^(0.5)))
%(五)启动性能
%100
K(1)=3;
%101
precisionist=1;
K(1)=3;Ist(1)=K(1)*Tms*IKW
while precisionist>0.03
for k=1:3
Fst(k)=Ist(2*k-1)*Ns1/a1*0.707*(KU1+Kd1^2*Kp1*Z1/Z2)*(1-epsilon0)^(0.5) beta0=0.64+2.5*(delta/(t1+t2))^(0.5)
BL=mu0*Fst(k)/(2*delta*beta0)
%102
KZ(k)=input('输入启动漏抗饱和KZ（k）=')
%103
cs1(k)=(t1-b01)*(1-KZ(k))
cs2(k)=(t2-b02)*(1-KZ(k))
%105
deltalambdaU1(k)=(h01+0.58*h11)/b01*(cs1(k)/(cs1(k)+1.5*b01)) lambdas1st(k)=KU1*(lambdaU1-deltalambdaU1(k))+KL1*lambdaL1 %106
Xs1sts(k)=lambdas1st(k)/lambdas1*Xs1s
%107
Xdelta1sts(k)=KZ(k)*Xdelta1s
%108
Xsigma1sts(k)=Xs1sts(k)+Xdelta1sts(k)+XE1s
%109
hB=h22;
xi=1.987*10^(-3)*hB*(f/rhow1)^(0.5)
%110
b1b2=b22/b12;
KF(k)=input('输入电阻增加系数KF(k)=')
Kx(k)=input('输入电抗增加系数Kx(k)=')
%111
deltalambdaU2(k)=h02/b02*(cs2(k)/(cs2(k)+b02))
lambdas2st(k)=(lambdaU2-deltalambdaU2(k))+Kx(k)*lambdaL2
%112
Xs2sts(k)=lambdas2st(k)*Xs2s/lambdas2
%113
Xdelta2sts(k)=KZ(k)*Xdelta2s
%114
Xksts(k)=KZ(k)*Xsks
%115
Xsigma2sts(k)=(Xs2sts(k)+Xdelta2sts(k)+Xksts(k)+XE2s)*Cx
%116
Xsigmasts(k)=Xsigma1sts(k)+Xsigma2sts(k)
%117
R2sts(k)=KF(k)*RBs+RRs
%118
Rsts(k)=R1s+R2sts(k)
Zsts(k)=(Rsts(k)^2+(Xsigmasts(k))^2)^(0.5)
%120
Ist(2*k)=IKW/Zsts(k)
precisionist=abs((Ist(2*k)-Ist(2*k-1))/Ist(2*k))
Ist(2*k+1)=input('输入Ist(2*k+1)=')
end
end
%121
ist=Ist(6)/I1(2)
%122
Tsts=R2sts(3)/Zsts(3)^2*(1-sN)
MATLAB程序设计教程(第二版)课后答案(可编辑修改word版)
MATLAB 第二版课后答案unit3-8 unit3 实验指导 1、 n=input('请输入一个三位数：'); a=fix(n/100); b=fix((n-a*100)/10); c=n-a*100-b*10; d=c*100+b*10+a 2（1） n=input('请输入成绩'); switch n case num2cell(90:100) p='A'; case num2cell(80:89) p='B'; case num2cell(70:79) p='C'; case num2cell(60:69) p='D'; otherwise p='E'; end price=p （2）n=input('请输入成绩'); if n>=90&n<=100 p='A'; elseif n>=80&n<=89 p='B'; elseif n>=70&n<=79 p='C'; elseif n>=60&n<=69 p='D'; else p='E'; end price=p （3）try n; catch price='erroe' end 3 n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6];
b=n(1); for m=2:20 if n(m)>a a=n(m); elseif n(m)=0 disp(A([n],:)); elseif n<0 disp(lasterr); else disp(A([6],:)); disp(lasterr); end 7（1） f=[];
三相异步电动机Matlab仿真
中国石油大学胜利学院综合课程设计总结报告 题目：三相异步电机直接启动特性实验模型 学生姓名：潘伟鹏 系别：机械与电气工程系 专业年级： 2012级电气工程专业专升本2班 指导教师：王铭
2013年 6 月 27日
一、设计任务与要求 普通异步电动机直接起动电流达到额定电流的6--7倍，起动转矩能达到额定转矩的1.25倍以上。过高的温度、过快的加热速度、过大的温度梯度和电磁力，产生了极大的破坏力，缩短了定子线圈和转子铜条的使用寿命。但在电网条件和工艺条件允许的情况下，异步电动机也可以直接启动。本次课程设计通过MATLAB软件建模模拟三相异步电动机直接启动时的各个元器件上的电量变化。 参考： 电力系统matlab仿真类书籍 电机类教材 二、方案设计与论证 三相异步电动机直接起动就是利用开关或接触器将电动机的定子绕组直接接到具有额定电压的电网上。 由《电机学》知三相异步电动机的电磁转矩M与直流电动机的电磁转矩有相似的表达形式。它们都与电机结构（表现为转矩常数）和每级下磁通有关，只不过在三相异步电动机中不再是通过电枢的全部电流，而是点数电流的有功分量。三相异步电机电磁转矩的表达式为： （1-1）式中——转矩常数 ——每级下磁通 ——转子功率因数 式(1-1)表明，转子通入电流后，与气隙磁场相互作用产生电磁力，因此，反映了电机中电流、磁场和作用力之间符合左手定则的物理关系，故称为机械特性的物理表达式。该表达式在分析电磁转矩与磁通、电流之间的关系时非常方便。 从三相异步电动机的转子等值电路可知， （1-2） （1-3）将式(1-2)、(1-3)代入(1-1)得：
电机设计方法
第2章电磁场有限元分析简介 电磁场的边值问题实际上是求解给定边界条件下的麦克斯韦（Maxwell）方程组及由方程组深化出的其他偏微分方程问题。从求解问题的技术手段上来说，它可以分为解析求解和数值求解两大类。对于简单模型，有时可以得到方程的解析解。若模型复杂度增加，则往往很难获得模型的解析解。随着计算工具，特别是高速大容量电子计算机的发展，电磁场数值分析已深入到工业生产各个领域，解决问题的面越来越广，分析的问题也日趋复杂。电磁场数值分析是一门综合性的学科，涉及电磁场理论、数值分析、计算方法、计算机基础知识及高级语言等多个方面，但在计算上存在着共性。有限元法是一种常用的数值方法，并有相应的电磁软件问世，其中ANSOFT公司的Maxwell 3D/2D就是非常优秀的电磁分析软件。 本章将对电磁场的基本理论、电磁场有限元的求解及ANSOFT公司的Maxwell 3D/2D 作简单的介绍。至于完整的电磁理论描述，读者可以参考许多教科书。如果读者已熟悉电磁理论，完全可以略过本章，直接从第2章开始学习如何使用Maxwell电磁软件。 1.1电磁场基本理论 1.1.1麦克斯韦方程组 在19世纪中叶，麦克斯韦在总结前人工作的基础上，提出了适用于所有宏观电磁现象的数学模型，称之为麦克期韦方程组。它是电磁场理论的基础，也是工程电磁场数值分析的出发点。 麦克斯韦方程组包括微分和积分两种形式，在此仅给出它们的微分形式，通过它们可以导出能用有限元处理电磁问题的微分方程。 麦克斯韦方程组为 法拉第电磁感应定律 麦克斯韦-安培定律 高斯电通定律 高斯磁通定律 电荷守恒定律
式中，E为电场强度，V/m；D为电通量密度，C/m；H为磁场强度，A/m；B为磁通量密度，T；J为电流密度，A/㎡；P为电荷密度C/m3。 上面5个方程中包含两个旋度方程式（1.1）、式（1.2）和3个散度方程式（1.3）、式（1.4）和式（1.5）。
MATLAB程序设计教程课后答案
实验指导 1、 n=input('请输入一个三位数：'); a=fix(n/100); b=fix((n-a*100)/10); c=n-a*100-b*10; d=c*100+b*10+a 2（1） n=input('请输入成绩'); switch n case num2cell(90:100) p='A'; case num2cell(80:89) p='B'; case num2cell(70:79) p='C'; case num2cell(60:69) p='D'; otherwise p='E'; end price=p （2）n=input('请输入成绩'); if n>=90&n<=100 p='A'; elseif n>=80&n<=89 p='B'; elseif n>=70&n<=79 p='C'; elseif n>=60&n<=69 p='D'; else p='E'; end price=p （3）try n; catch price='erroe' end 3 n=[1,5,56,4,3,476,45,6,3,76,45,6,4,3,6,4,23,76,908,6]; a=n(1)；
for m=2:20 if n(m)>a a=n(m); elseif n(m)=0 disp(A([n],:)); elseif n<0 disp(lasterr); else disp(A([6],:)); disp(lasterr); end 7（1） f=[]; for n=1:40
直流电动机的MATLAB仿真..
第一章课程设计内容及要求 1. 直流电动机的机械特性仿真； 2. 直流电动机的直接起动仿真； 3. 直流电动机电枢串联电阻启动仿真； 4. 直流电动机能耗制动仿真； 5．直流电动机反接制动仿真； 6. 直流电动机改变电枢电压调速仿真； 7. 直流电动机改变励磁电流调速仿真。 要求：编写M文件，在Simulink环境画仿真模型原理图，用二维画图命令画仿真结果图或用示波器观察仿真结果，并加以分析
第二章直流电动机的电力拖动仿真绘制 1）直流电动机的机械特性仿真 clear; U_N=220;P_N=22;I_N=115; n_N=1500;R_a=;R_f=628; Ia_N=I_N-U_N/R_f; C_EPhi_N=(U_N-R_a*Ia_N)/n_N; C_TPhi_N=*C_EPhi_N; Ia=0;Ia_N; n=U_N/C_EPhi_N-R_a/(C_EPhi_N)*Ia; Te=C_TPhi_N*Ia; P1=U_N*Ia+U_N*U_N/R_f; T2_N=9550*P_N/n_N; figure(1); plot(Te,n,'.-'); xlabel('电磁转矩Te/'); ylabel('转矩n/rpm'); ylim([0,1800]); figure(2); plot(Te,n,'rs'); xlabel('电磁转矩Te/'); ylabel('转矩n/rpm');
hold on; R_c=0; for coef=1:;; U=U_N*coef; n=U/C_EPhi_N-(R_a+R_c)/(C_EPhi_N*C_TPhi_N)*Te; plot(Te,n,'k-'); str=strcat('U=',num2str(U),'V'); s_y=1650*coef; text(50,s_y,str); end figure(3); n=U_N/C_EPhi_N-(R_a+R_c)/(C_EPhi_N*C_TPhi_N)*Te; plot(Te,n,'rs'); xlabel('电磁转矩Te/'); ylabel('转矩n/rpm'); hold on; U=U_N;R_c=; for R_c=0::; n=U/C_EPhi_N-(R_a+R_c)/(C_EPhi_N*C_TPhi_N)*Te; plot(Te,n,'k-'); str=strcat('R=',num2str(R_c+R_a),'\Omega'); s_y=400*(4-R_c*; text(120,s_y,str);
刘亚敏1520310052--电机现代设计方法与优化作业
电机现代设计方法与优化作业 电气工程刘亚敏 1520310052 1、所用算法的寻优策略 本篇论文所采用的算法为蚁群算法，又称蚂蚁算法，其定义为：各个蚂蚁在没有事先告诉他们食物在什么地方的前提下开始寻找食物。当一只找到食物以后，它会向环境释放一种挥发性分泌物pheromone (信息素,该物质随着时间的推移会逐渐挥发消失，信息素浓度的大小表征路径的远近)来实现的，吸引其他的蚂蚁过来，这样越来越多的蚂蚁会找到食物。有些蚂蚁并没有像其它蚂蚁一样总重复同样的路，他们会另辟蹊径，如果另开辟的道路比原来的其他道路更短，那么，渐渐地，更多的蚂蚁被吸引到这条较短的路上来。最后，经过一段时间运行，可能会出现一条最短的路径被大多数蚂蚁重复着。 蚁群算法是一种新型的模拟进化算法，该算法通过模拟蚂蚁觅食的方式，使一定数量的蚂蚁在解空间内进行随机搜索，对路径上蚂蚁释放的信息素进行更新，按照转移概率决定前进的方向，最后收敛于全局最优解。蚁群算法具有较强的鲁棒性。相对于其它算法，蚁群算法对初始路线要求不高，即蚁群算法的求解结果不依赖子初始路线的选择，而且在搜索过程中不需要进行人工的调整。其次，蚁群算法的参数数目少，设置简单，易于蚁群算法应用到其它组合优化问题的求解。现在蚁群算法己经在电力网络优化、网络路中分配、函数优化和集成电路布线等领域得到应用。本文将蚁群算法进行了改进，将其用于永磁同步电机的优化设计中。
2、论文对算法的改进 算法与其它智能优化算法相比，存在搜索时间长的缺陷，该算法的复杂度可以反映这一点；而且该算法容易出现停滞现象，即搜索到一定程度后，所有个体发现的解完全一致，不能对解空间进行进一步的搜索，不利于发现更好的解。本文借鉴蚁群算法的进化思想，针对以上提及的两个问题，将算法的数学模型做了三方面的改进。 2.1转移规则的改进 对每只蚂蚁i ，定义其函数值为相应的目标函数值Zi ，并记蚂蚁i 与蚂蚁j 的目标函数值的差值为 蚂蚁j 到蚂蚁i 的转移概率为 式中:———蚂蚁j 邻域内的信息素数量； α和β———算法的权重因子，本文取α=β =1。 2.2信息素更新规则的改进 由于信息素强度Q 是表征蚂蚁所经轨迹数量的一个常数，它影响算法的全局收敛速度[ 5]。蚂蚁之间通过信息素进行交流，因此, 本文针对蚁群算法寻优过程易陷入局部最小的弊端，提出根据算法搜索的情况，动态修改需要增加的信息素的方法。即用时变函数Q(t)
直流电机模糊控制系统的MATLAB-Simulink仿真研究毕业设计
XXXX届毕业设计说明书 直流电机模糊控制系统 的MATLAB／Simulink仿真研究 院、部：电气与信息工程学院 学生姓名：XXX 指导教师：XXXX职称教授 职称 专业：XXXXXXXXXXXXX 班级：XXXXXXXXX 完成时间：20XX.X.X
摘要 在当今控制技术的发展当中，模糊控制技术的发展走在了前列，成为了当今世界上最先进的控制技术之一。模糊控制技术很好的将模糊数学理论应用于控制领域当中, 更加真切地模拟出了人脑的思维方式和判断能力, 以及对产品生产的过程进行筛选和对产品质量上的控制, 从而发展出了基于模糊控制技术的智能化的新技术，为当今控制技术的发展提供了广阔空间。 在本文当中，主要介绍了基于模糊控制理论的直流电机模糊控制系统的原理，以及直流电机模糊控制系统的优点和缺点，并通过使用MATLAB语言中SIMULINK 模块和模糊控制工具箱对直流电机模糊控制系统进行仿真，把控制直流电机调速的实际情况转换成模糊控制规则，再使用这些规则，对过程经过模糊推理和模糊决策所得到的控制量，从而实现在MATLAB语言中SIMULINK模块和模糊控制工具箱对直流电机模糊控制系统的建模与仿真。对仿真结果予以分析，对直流电机模糊控制系统的仿真进行总结。 关键词：MATLAB；SIMULINK；模糊控制；直流电机；电机调速
ABSTRACT Among today’s control technology development, one of the leading enterprises in the development of fuzzy control technology, fuzzy control technology has become one of the most advanced control technology in the world today, it will be a very good fuzzy control technology of fuzzy mathematics theory is applied in control field, the more realistically simulate the human brain’s way of thinking and judgment ability, as well as to the production process of screening and the control on the quality of product, which was developed based on fuzzy intelligent control technology of the new technology, for the development of modern control technology provides a broad expansion of space. in this article, mainly introduced the dc motor based on fuzzy control theory, the principle of fuzzy control system, as well as the advantages and disadvantages of the fuzzy control system for dc motor, and by using the SIMULINK module and the fuzzy control toolbox in MATLAB language for the calculation of the fuzzy control system of dc motor, the control of the actual situation of the dc motor speed control is converted into fuzzy control rules, and then use these rules, the process through fuzzy reasoning and fuzzy decision of control, thus to achieve the SIMULINK module and the fuzzy control toolbox in MATLAB language modeling and simulation of fuzzy control system of a dc motor. And the analysis to the results of simulation and simulation of fuzzy control system of dc motor. Keywordsmatlab；Simulink；fuzzy control；dc motor；motor speed control
MATLAB程序设计教程(第二版)第三章实验报告下载
大学社区网收集整理https://www.360docs.net/doc/4516438309.html, 评分 日期湖南商学院北津学院实验报告 课程名称MATLAB科学计算编程语言 实验名称MATLAB程序设计 专业班级信科1121班 姓名xxx 学号xxx 实验日期2012年11月5日 2012—2013学年度第一学期 一、实验目的 1.掌握利用if语句、switch语句实现选择结构的方法。 2.掌握利用for语句、while语句实现循环结构的方法。 3.熟悉利用向量运算来代替循环操作的方法并理解MATLAB程序设计的特点 4.掌握定义和调用MATLAB函数的方法。
二、实验环境 系统windows7旗舰版 处理器Intel(R)Core(TM)i7-3610M CPU @ 2.30GHz 安装内存 4.00GB （3.07GB 可用）系统类型64位操作系统运行环境 MATLAB 5.3 三、实验基本原理 利用上课所学知识解决以下问题： 1.从键盘输入一个3位数的整数，将它反向输出。如输入639，输出936。 2.输入一个百分制成绩，要求输出成绩等级A、B、C、D、E。其中90~100分为A，80~89分为B，70~79分为C，60~69分为D，60分以下为E。 要求： (1)分别用if 语句和switch 语句实现。 (2)输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。3.输入20个数，求其中最大数和最小数。要求分别用循环结构和调用MATLAB 的max 函数、min 函数来实现。 4.23.0ln )3.0sin(23.03.0a a e e y a a +++?=?，当a 取-3.0、-2.9、-2.8、…、2.8、2.9、 3.0时，求各点的函数值。要求分别用顺序结构和循环结构实现。 5.当n 分别取100、1000、10000时，求下列各式的值： (1)) 6...(n 1...31211122222π=+++++(2)) 2...()12)(12()2)(2(...756653443122π =??? ?????+?????????××????????××????????××n n n n 要求分别用循环结构和向量运算（使用sum 函数）来实现。 6.建立5×6矩阵，要求输出矩阵第n 行元素。当n 值超过矩阵的行数时，自动转为输出矩阵最后一行元素，并给出出错信息。 7已知,o999 ) 20()30()40(f f f y += (1)当)5ln(10)(2 ++=n n n f 时，y 的值是多小。 (2)当+×+×+×=433221)(n f …+)1(+×n n 时，y 的值是多小。 8.先用函数的递归调用定义一个函数文件求 ∑=n i m i 1，然后调用该函数文件求
异步电动机机械特性的MATLAB仿真
辽宁工业大学 实验室开放课题设计（论文） 题目：异步电动机机械特性的MATLAB仿真》 院（系）：电气工程学院 专业班级：自动化 131 学号： 0 ` 学生姓名：徐峰 指导教师：赵丽丽
起止时间：
摘要 异步电动机以其结构简单、运行可靠、效率较高、成本较低等特点，在日常生活中得到广泛的使用。目前，电动机控制系统在追求更高的控制精度的基础上变得越来越复杂，而仿真是对其进行研究的一个重要手段。MATLAB是一个高级的数学分析和运算软件，可用动作系统的建模和仿真。在分析三相异步电动机物理和数学模型的基础上，应用MATLAB软件简历了相对应的仿真模型；在加入相同的三相电压和转矩的条件下，使用实际电机参数，与MALAB给定的电机模型进行了对比仿真。 第一章对异步电机的实验要求做出了相关的描述，第二章对MATLAB仿真软件做了一定的介绍，第三章是对异步电动机的机械特性、启动、制动和正反转进行理论分析和仿真模拟以及仿真结果的分析。 经分析后，表明模型的搭建是合理的。因此，本设计将结合MATLAB的特点，对三相异步电机进行建模和仿真，并通过实际的电动机参数，对建立的模型进行了验证。 关键词：异步电机、数学模型、MATLAB仿真、三相异步电动机
目录 第1章实验任务及要求 (1) 第2章 MATLAB及SIMULINK的介绍 (2) MATLAB介绍 (2) S IMULINK模块的介绍 (3) 第3章仿真实验 (4) 三相异步电动机的机械特性 (4) 三相异步电动机起动的仿真 (6) 三相异步电动机制动仿真 (8) 三相异步电动机正反转仿真 (10) 第4章总结 (12) 参考文献 (13) 附录 (14)
基于Matlab的直流电机速度控制
基于Matlab的直流电机速度控制
系统仿真 课程设计报告 设计题目：基于Matlab的直流电机速度控制 专业：自动化 学生姓名： 班级学号： 指导教师： 开课日期2013年 7 月 1 日至2013年 7 月 13 日南京邮电大学自动化学院
一、课程设计题目 控制系统的执行机构常用直流电机来驱动，电路和原理示意图如下所示 其开环传递函 数 为 ()()0001 .0)15.0)(1.001.0(01 .02+++= +++= s s K R Ls b Js K V θ ，请用时域分析方法设计PID 控制器，使系统满足下列性能指标要求：当仿真输入是单位阶跃信号时，电机输出转速调整时间小于2秒，超调小于5％，稳态误差小于1％。 要求给出详细的设计步骤，matlab 源码及仿真曲线。 二、实验原理 本报告首先介绍了直流电动机的物理模型，并测量计算了它的具体参数。然后根据牛顿第二定律和回路电压法分别列写运动平衡方程式和电机电枢回路方程式，从而通过一些数学变换抽象出了以电压为输入、转速为输出、电流和转速为状态变量的数学模型。借助MATLAB 设计simulink 模块调整PID 模块的各项系数，使系统的阶跃响应达到了设计指标。 1、建立该系统的时域数学模型 由克希霍夫定律得： V=R*i+L +e 直流电机转矩和电枢电流关系为 T=Kt*I 电枢旋转产生反电动势与旋转运动角速度的关系为 e=
由牛顿定律，转子力矩平衡关系为 其中，T：负载转矩，：负载电流 V(s)=R*I(s)+L*sI(s)+E(s) 拉式变换：E=Ke(s) 划去中间变量得： 开环传递函数为： 2、PID控制器的功能 比例环节：Kp增大等价于系统的开环增益增加，会引起系统响应速度加快，稳态误差减少，超调量增加。当Kp过大时，会使闭环系 统不稳定； 积分环节：相当于增加系统积分环节个数，主要作用是消除系统的稳态误差。积分环节作用的强弱取决于积分时间常数Ti，Ti增大， 系统超调量变小，响应速度变慢； 微分环节：主要作用是提高系统的响应速度，同时减少系统超调量，抵消系统惯性环节的相位滞后不良作用，使系统稳定性明显改善。 Td偏大或偏小，都会使超调量增大，调整时间加长。由于该环节所产 生的控制量与信号变化速率有关，故对于信号无变化或变化缓慢的系 统微分环节不起作用。 三、设计步骤 方法1： 搭建simulink模块，利用经验调节法整定PID参数，使整个系统满足调节时间小于2秒，超调小于5%，稳态误差小于1%。 1、搭建的simulink模块图如下：
玩具设计中马达的选用技巧和方法
玩具设计中马达的选用技巧和方法(总2页) -CAL-FENGHAI.-(YICAI)-Company One1 -CAL-本页仅作为文档封面，使用请直接删除
玩具设计中马达的选用技巧和方法 马达是玩具能够保持持续动力的主要动力源。设计人员在玩具生产中很少有机会设计一个马达，主要是选用。玩具马达的选用也有很多技巧和方法。 在玩具制造中常用的马达（电机），有万宝制（Mabuchi）、标准（Standard）等，这些马达生产厂家的产品都有马达性能的参照表，设计人员很少有机会设计一个马达，主要是选用，当然也可能因为某种原因而专门设计一个马达（很多马达厂的工程师都是从Mabuchi出来的，所以Mabuchi的标准基本可以用作玩具业的马达标准）。 用试验找出合适的转速比和扭矩 马达厂的标准是空载的转速，而转速比计算时，一般就是用空载转速的计算。行走类玩具的转速比在80~150之间，因玩具产品并非是一种精密的机械，所以可以通过调节齿轮的传动来得出实际需要的传动比。比如：齿轮传动中某一节的齿数比为20：8，实验发现转速太慢则可以试一下22：6，如果转速太快，又可以换18：10，这种办法非常方便。 马达产品都有一个最大载荷的数据，还有一个载荷与转速的关系。因为玩具产品一般不会很大而且是用电池驱动的较多，所以马达的载荷也不会很大。一般来说，载荷与扭矩越大，马达的体积也越大，所要的电池也越多，而且玩具的重量也越重，所以尽量选用合适的马达。一般保险起见，载荷与转速在曲线上升区的马达就管用。如果要节约成本，在找到一个大功率马达就可以起作用的情况下，换一个小功率的试一下。通过反复试验，找到一个功率较小又能满足要求的马达。 用收音机做干扰测试 马达的工作原理是不断地通过电刷来改变线圈中的电流，从而保证连续的转动。可以想像，马达里的磁铁越大，线圈越长，则扭矩会越大，而这种电流的改变会形成一个电磁波。电磁波的存在会对很多种家用电器造成干扰。如果一个小孩在等飞机的时候玩带有这种马达的玩具，还有可能影响飞机的起降。所以在实际工作中会增加一个抗干扰的垫圈（如DV WASHER , 0.5港币左右），但一般设计时都会在马达的两个电极上并联上一个陶瓷片电容或一个电阻（稳定电路的作用），一方面减少电刷与线圈转换时产生的电流影响，另一方面起稳定作用。实际工作中，用一个小收音机来做干扰测试器材，如不影响所有波段的电台就基本可以了。 马达转动轴一般是铁轴，而铁轴一般直接紧配一个6齿或8齿的0.5模数的齿轮（或一个小带轮或锅杆），当是齿轮或锅杆时，特别是锅杆时，会因为刚起动的瞬间，齿与齿之间会卡位，这时线圈中的电流会加大，而导致铁轴串动，串动就会解卡，所以在固定马达时千万不要把轴在轴的方向上固定死（致少要有0.5mm的余量）。
电机设计matlab程序
%电机设计程序 clear all format short e m1=3;p=2;f=50 %1.额定功率 PN=*10^3 ; %2.额定电压（单位V，三角形接法） UN=380;UN0=380; %3.功电流(单位A) IKW=PN/(m1*UN0) %4.效率eta按照技术条件的规定eta= eta= ; %5.功率因数cos(phi) =,按照技术条件的规定cos(phi)= phi=acos; cos(phi); %6.极对数p=2 p=2; %7.定转子槽数：每极每相槽数取整数。参考类似规格电机取q1=3,则Z1=2m1pq1，再查表10-8选Z2=32，并采用转子斜槽。 q1=3; Z1=2*m1*p*q1 Z2=32 ; %8.定转子每极槽数 Zp1=Z1/(2*p) Zp2=Z2/(2*p) %9.确定电机的主要尺寸；一般可参考类似电机的主要尺寸来确定Di1和lef.现按10-2中的 KE1=*log(PN/1000)*p+ P1=KE1*PN/(eta*cos(phi)) alphap1=;KNm1=;Kdp1=;A1=25000; Bdelta1=;n1=1450; V=(alphap1*KNm1*Kdp1))*(1/(A1*Bdelta1 ))*(P1/n1) D1=; %铁心的有效长度 Di1=; lef =V/((Di1)^2) %气隙的确定 %参考类似产品或由经验公式（10-10a），得 lt=;
delta = lef=lt + 2*delta D2=Di1-2*delta %转子内径先按转轴直径决定（以后再校验转子轭部磁密） Di2= ; %11.极距 tau tau =pi*Di1/(2*p) %12.定子齿距t1 t1=(pi*Di1/Z1) %转子齿距t2 t2=(pi*D2/Z2) bsk=; %15.设计定子绕组 Nphi11=eta*cos(phi)*pi*Di1*A1/(m1*IKW) %取并联支路a1=1,由式（10-15），可得每槽导体数 a1=1; Ns1=47 %16.每相串联导体数Nphi1 Nphi1=Ns1*Z1/(m1*a1) %每相串联匝数N1 N1=Nphi1/2 %17.绕组线规设计 %初选定子电密J11=5.0A/mm^2,由式（10-16），计算导线并绕根数和每根导线面积的乘积。 J11=; %其中定子电流初步估计值 I11=IKW/(eta*cos(phi)) Nt1Ac11=I11/(a1*J11)
经典-同步电机模型的MATLAB仿真h
安徽工业大学工商学院课程设计（论文）同步电机模型的MATLAB仿真 学生姓名：李春笋 学号：111842161 专业班级：气1142 指导教师：范国伟 2013年12月20日
摘要 采用电力电子变频装置实现电压频率协调控制，改变了同步电机历来的恒速运行不能调速的面貌，使它和异步电机一样成为调速电机大家庭的一员。本文针对同步电机中具有代表性的凸极机，在忽略了一部分对误差影响较小而使算法复杂度大大增加的因素（如谐波磁势等），对其内部电流、电压、磁通、磁链及转矩的相互关系进行了一系列定量分析，建立了简化的基于abc三相变量上的数学模型，并将其进行派克变换，转换成易于计算机控制的d/q坐标下的模型。再使用MATLAB中用于仿真模拟系统的SIMULINK 对系统的各个部分进行封装及连接，系统总体分为电源、abc/dq转换器、电机内部模拟、控制反馈四个主要部分，并为其设计了专用的模块，同时对其中的一系列参数进行了配置。系统启动仿真后，在经历了一开始的振荡后，各输出相对于输出时间的响应较稳定。关键词：同步电机 d/q模型 MATLAB SIMULINK 仿真。 The Simulation Platform of Synchronous Machine by MATLAB Abstract: The utilization of transducer realizes the control of voltage’s frequency. It changes the situation that Synchronous Machine is always running with constant speed. Just like Asynchronous Machine, Synchronous machine can also be viewed as a member of the timing machine. This thesis intends to aim at the typical salient pole machine in Synchronous Machine. Some quantitative analysis are made on relations of salient pole machine among current, voltage, flux, flux linkage and torque, under the condition that some factors such as harmonic electric potential are ignored. These factors have less influence on error but greatly increase complexity of arithmetic. Thus, simplified mathematic model is established on the basis of a, b, c three phase variables. By the Park transformation, this model is transformed to d, q model which, is easy to be controlled by computer. Simulink is used to masking and linking all the parts of the system. The system can be divided into four main parts, namely power system, abc/dq transformation, simulation model of the machine and feedback control. Special blocks are designed for the four parts and a series of parameters in these parts are configured. The results of simulation show that each output has a satisfactory response when there is disturbance. Key Words: Synchronous Machine Simulation d/q Model MATLAB SIMULINK
matlab电机实例
unction [sys,x0,str,ts,simStateCompliance] = BLDC_S(t,x,u,flag) %-----------------------------------------------------------------------
% 状态变量：X(1)=ia;X(2)=ib;X(3)=ic;X(4)=SETA;X(5)=OMEGA; % 输入量：u(1)=Ud; u(2)=TL; % 输出量：n, Tem, ia, ib, ic; %------------------------------------------------------------------------ %-----------------电动机参数--------------------------------------------- R = 0.23;L = 0.00498; M = -0.00005478; J = 0.025; P0=2; % 极对数 RR = diag([R R R]); LL = diag([L-M,L-M,L-M]); S = [2,-1,-1; -1,2,-1; -1,-1,2]/3; %----------------------------------------------------------------------- %SFUNTMPL General M-file S-function template % With M-file S-functions, you can define you own ordinary differential % equations (ODEs), discrete system equations, and/or just about % any type of algorithm to be used within a Simulink block diagram. % % The general form of an M-File S-function syntax is: % [SYS,X0,STR,TS,SIMSTATECOMPLIANCE] = SFUNC(T,X,U,FLAG,P1,...,Pn) % % What is returned by SFUNC at a given point in time, T, depends on the % value of the FLAG, the current state vector, X, and the current % input vector, U. % % FLAG RESULT DESCRIPTION % ----- ------ -------------------------------------------- % 0 [SIZES,X0,STR,TS] Initialization, return system sizes in SYS, % initial state in X0, state ordering strings % in STR, and sample times in TS. % 1 DX Return continuous state derivatives in SYS. % 2 DS Update discrete states SYS = X(n+1) % 3 Y Return outputs in SYS. % 4 TNEXT Return next time hit for variable step sample % time in SYS. % 5 Reserved for future (root finding). % 9 [] Termination, perform any cleanup SYS=[]. % % The state vectors, X and X0 consists of continuous states followed % by discrete states. % % Optional parameters, P1,...,Pn can be provided to the S-function and % used during any FLAG operation. % % When SFUNC is called with FLAG = 0, the following information % should be returned: % % SYS(1) = Number of continuous states. % SYS(2) = Number of discrete states. % SYS(3) = Number of outputs. % SYS(4) = Number of inputs. % Any of the first four elements in SYS can be specified % as -1 indicating that they are dynamically sized. The % actual length for all other flags will be equal to the % length of the input, U. % SYS(5) = Reserved for root finding. Must be zero. % SYS(6) = Direct feedthrough flag (1=yes, 0=no). The s-function % has direct feedthrough if U is used during the FLAG=3 % call. Setting this to 0 is akin to making a promise that % U will not be used during FLAG=3. If you break the promise % then unpredictable results will occur. % SYS(7) = Number of sample times. This is the number of rows in TS. % % X0 = Initial state conditions or [] if no states. %
Matlab程序设计教程(第二版)刘卫国课后参考答案解析
第二章 1 求下列表达式的值。 （1） w=sqrt(2)*(1+0.34245*10^(-6)) (2) a=3.5; b=5; c=-9.8; x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/tan(b+c)+a (3) a=3.32; b=-7.9; y=2*pi*a^(2)*[(1-pi/4)*b-(0.8333-pi/4)*a] (4) t=[2,1-3*i;5,-0.65]; z=1/2*exp(2*t)*log(t+sqrt(1+t^(2))) 2 求下列表达式 A=[-1,5,-4;0,7,8;3,61,7]; B=[8,3,-1;2,5,3;-3,2,0]; (1) A+6*B A^2-B+eye (2) A*B A.*B B.*A (3) A/B B\A (4) [A,B] [A([1,3],:);B^2] 3 根据已知，完成下列操作 （1） A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14]; K=find(A>10&A<25); A(K) （2） A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14]; B=A(1:3,:) C=A(:,1:2) D=A(2:4,3:4) E=B*C （3） E f=input('输入一个数：','s'); f(end :-1:1) 2 用if语句 score=input('请输入成绩:'); if score>=90&&score<=100 disp('A'); elseif score>=80&&score<=89 disp('B'); elseif score>=70&&score<=79 disp('C'); elseif score>=60&&score<=69; disp('D'); elseif score<60&&score>=0; disp('E'); else disp('出错'); end 用switch语句 score=input('请输入成绩:'); switch fix(score/10) case {9,10} disp('A'); case {8} disp('B'); case {7} disp('C'); case {6} disp('D'); case {0,1,2,3,4,5} disp('E'); otherwise disp('出错'); end 第四章1题 1) X=0:10; Y=x-x.^3/6; P lot(x,y)