(完整word版)电机设计matlab程序
(完整word版)MATLAB程序设计教程(第二版)课后答案
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+a2（1）n=input('请输入成绩’）；switch ncase num2cell（90：100)p='A'；case num2cell（80:89)p='B'；case num2cell(70：79）p=’C';case num2cell（60:69)p='D'；otherwisep='E';endprice=p(2）n=input(’请输入成绩');if n〉=90&n〈=100p='A’；elseif n>=80&n<=89p='B';elseif n〉=70&n<=79p=’C’;elseif n>=60＆n<=69p='D'；elsep='E'；endprice=p(3）tryn;catchprice='erroe'end3n=［1,5,56，4,3,476，45，6，3，76，45,6,4,3，6，4,23,76,908，6]; a=n(1);b=n(1)；for m=2:20if n（m）>aa=n（m）；elseif n（m)<bb=n（m）;endendmax=amin=b法2n=[1,5，56,4，3，476,45，6,3,76,45，6，4，3,6，4，23，76,908,6］;min=min(n)max=max(n）4b=［—3.0:0.1:3.0];for n=1：61a=b(n）；y(n)=（exp（0.3＊a)-exp(—0。
(完整word版)含答案《MATLAB实用教程》
第二章 MATLAB 语言及应用实验项目实验一 MATLAB 数值计算三、实验内容与步骤1.创建矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321a（1（2）用（3）用（42.矩阵的运算（1）利用矩阵除法解线性方程组。
⎪⎪⎩⎪⎪⎨⎧=+++=-+-=+++=+-12224732258232432143214321421x x x x x x x x x x x x x x x 将方程表示为AX=B ，计算X=A\B 。
（2）利用矩阵的基本运算求解矩阵方程。
已知矩阵A 和B 满足关系式A -1BA=6A+BA ，计算矩阵B 。
其中⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=7/10004/10003/1A ，Ps: format rata=[1/3 0 0;0 1/4 0;0 0 1/7];b=inv(a)*inv(inv(a)-eye(3))*6*a（3）计算矩阵的特征值和特征向量。
已知矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=1104152021X ，计算其特征值和特征向量。
（4）Page:322利用数学函数进行矩阵运算。
已知传递函数G(s)=1/(2s+1),计算幅频特性Lw=-20lg(1)2(2w )和相频特性Fw=-arctan(2w),w 的范围为[0.01，10]，按对数均匀分布。
3.多项式的运算（1）多项式的运算。
已知表达式G(x)=(x-4)(x+5)(x 2-6x+9),展开多项式形式，并计算当x 在[0，20]内变化时G(x)的值，计算出G(x)=0的根。
Page 324（2）多项式的拟合与插值。
将多项式G(x)=x 4-5x 3-17x 2+129x-180，当x 在[0，20]多项式的值上下加上随机数的偏差构成y1，对y1进行拟合。
对G(x)和y1分别进行插值，计算在5.5处的值。
Page 325 四、思考练习题1.使用logspace 函数创建0～4π的行向量，有20个元素，查看其元素分布情况。
Ps: logspace(log10(0),log10(4*pi),20) (2) sort(c,2) %顺序排列 3.1多项式1）f(x)=2x 2+3x+5x+8用向量表示该多项式，并计算f(10)值. 2）根据多项式的根[-0.5 -3+4i -3-4i]创建多项式。
MATLAB程序设计M文件程序控制结构程序调试程序举例
y=input('Please input y=:');
[rho,the]=tran(x,y);
rho
the MATLAB中，函数能够嵌套调用，即一种函数能够调用别旳函数， 甚至调用它本身。一种函数调用它本身称为函数旳递归调用。
例6 分别建立命令文件和函数文件，将华氏温度f转换为 摄氏温度c。
21.1111
x=
21.1111
(5). 函数文件旳其他有关概念
①局部变量和全局变量：
局部变量只存在于单个函数工作空间，全局变量 经过global定义，可穿行于不同函数工作空间， 涉及基本工作空间workspace。
函数文件旳内部变量是局部旳，与其他函数 及MATLAB内存相互隔离； 而假如在若干函数中把某一变量定义为全局 变量，那么这些函数将公用这个变量，全局 变量旳作用域是整个MATLAB旳工作区，即 全程有效，全部函数都能够对其存取和修改。
if (nargin == 1) tol = max(size(x)) * max(s) * eps;
程序部分
end
r = sum(s > tol);
（2）命令m文件建立及其运营
建立 涉及下列环节： 进入m文件编辑器 输入程序 定义文件名，保存程序
命令M文件旳运营方式： 直接在命令窗口输入该文件旳文件名 在m文件编辑器中打开该文件后点击工具条中
② M文件模式
将matlab语句构成旳程序存储成以m为扩展名 旳文件，然后再执行该程序文件，这种工作模式 称为程序文件模式。
程序文件不能在命令窗口下建立，因为命令窗口 只允许一次执行一行上旳一种或几种语句。
开关磁阻电机matlab程序设计
开关磁阻电机设计程序%开关磁阻电机设计程序%%电机给定数据：功率7.5KW，电源电压380V，额定转速1500r/min，额定效率0.88，%调速范围200——2000r/min，4相8/6极结构。
%clear allformat short e%1.功率，转速PN=7.5*10^3;n=1500;%2.相数q=4;%3.定子极数、转子极数Ns=8;Nr=6;%4.绕组端电压(单位V)U=280;%5.主要尺寸选择，电磁功率eta=0.88;Pem=PN*((1+eta)/(2*eta));%6.细长比lambda=1.2;%7.电磁负荷初选值A1=28000;Bdelta1=0.4;%8.转子外径Daki=0.5;km=0.8;Da=6.1/(Bdelta1*A1)*(ki/km)*(Pem/n)/(1.05*lambda);Da=Da^(1/3);%9.铁芯叠长la=0.1355;%10.定子外径Ds=0.21;%11.气隙g=0.0004;%12.定转子极弧betas=21*pi/180;betar=23*pi/180;%13.定转子极宽bps=(Da+2*g)*sin(betas/2);bpr=Da*sin(betar/2);%14.第二气隙，设%%%%%%%%%%%% 此处可以输入数值 %%%%%%%%%%%%%%%%%%%例如：gi=input('输入第二气隙gi=')gi=0.01615;%15.定转轭高hcs=1.3*bps/2;hcr=1.4*bpr/2;%16.轴径。
%%%%%%%%%%%% 此处可以输入数值 %%%%%%%%%%Di=0.05;%17.定子槽深ds=(Ds-Da-2*g-2*hcs)/2;%19.有效铁芯长度kFe=0.93;lFe=kFe*la;%20.转子极距taur=pi*Da/Nr;%21.控制参数。
开通脚，关断角，导通角thetaon=0;thetau=0;thetaoff=(1/2)*(2*pi/Nr-betar)+thetau;thetahr=thetaoff;thetac=thetaoff-thetaon;%22.每相绕组串联匝数，取Bps=1.6T，重新校验BdeltaBps=1.6;ldelta=1.05*la;Bdelta=0.805*bps*Bps/taur;Nph=3.04*Nr*U*thetac/(n*Bdelta*Da*ldelta);%%%%%%%%%%%% 此处可以输入数值 %%%%%%%%%%%23.磁化曲线计算，选择DR510-50硅钢片冲片，delta=0.0005，计算机计算得：%不对齐位置电感Lu=0.00918; %%%%%%%%%%%% 此处可以输入数值 %%%%%%%%%%%25.额定电磁转矩omega=2*pi*n/60;Tem=Pem/omega;%26.额定电磁转矩Tem所需的理想方波电流幅值Ims%在图6-17中，取Ims=28A时，W1=13.61J，Ims=28;W1=13.61;Tem1=Ns*Nr*W1/4/pi; %%%%%%%%%%%% 此处可以输入数值 %%%%%%%%%%%27.绕组电流有效值I=Ims/(2^(1/2));%28.实际电流峰值ki=0.5;im=I/ki;%29.定子极间窗口面积SW=(1/2)*(pi/Ns)*((Ds/2-hcs)^2-(Da/2+g)^2)-bps*ds/2;%30.导线规格、导线截面积Sa和槽满率ks%预取电流密度J1=5A/mm^2,J1=5;Sa1=I/J1;Sa=3.94; %%%%%%%%%%%%%%%%%%%%此处可以输入数值%%%%%%%%%%%每槽导体净截面积SCu=Nph*Sa/2;%槽满率ks=SCu/SW;%31.电流密度J=5; %%%%%%%%%%%%%%%%%%%%此处可以输入数值%%%%%%%%%%%32.绕组平均匝长lavbw=(Da+2*g)*sin(pi/Ns)-bps/2;a=bps+bw;b=la+2*5*10^(-3)+bw;r=5*10^(-3)+bw/2;lav=2*la+2*(bps-2*5*10^(-3))+2*pi*r;%33.每相绕组导线总长l=Nph*lav;%34.每相绕组电阻Rp=0.0217*l/Sa;%参数计算%35.铜重GCu=q*Sa*l*10^(3)*8.9*10^(-6);%36.定子铁芯体积VSFe=((1/2)*(pi/Ns)*((Ds/2)^2-(Ds/2-hcs)^2)+bps*ds)*2*Ns*lFe;%37.转子铁芯体积VRFe=((1/2)*(pi/Nr)*((Di/2+hcr)^2-(Di/2)^2)+bpr*(gi-g)/2)*2*Nr*lF e;%38.铁芯重GFe=(VSFe+VRFe)*10^(9)*7.8*10^(-6);%39.电负荷Dsi=Da+2*g;A=q*Nph*I/(pi*Dsi);%40.铜耗PCu=q*I^2*Rp;%41.电机利用系数TN=PN/(2*n*pi/60);K=TN/(Da^2*la);。
电力系统分析潮流计算课程序设计及其MATLAB程序设计-范本模板
电力系统分析潮流计算程序设计报告题目:13节点配电网潮流计算学院电气工程学院专业班级学生姓名学号班内序号指导教师房大中提交日期 2015年05月04日目录一、程序设计目的 (1)二、程序设计要求 (3)三、13节点配网潮流计算 (3)3.1主要流程................................................................................................... 错误!未定义书签。
3。
1.1第一步的前推公式如下（1—1）-（1—5)： ................................. 错误!未定义书签。
3。
1.2第二步的回代公式如下（1-6)—（1-9)： ..................................... 错误!未定义书签。
3.2配网前推后代潮流计算的原理 (7)3。
3配网前推后代潮流计算迭代过程 (7)3.3计算原理 (8)四、计算框图流程 (9)五、确定前推回代支路次序.......................................................................................... 错误!未定义书签。
六、前推回代计算输入文件 (10)主程序： (10)输入文件清单： (11)计算结果: (12)数据分析: (12)七、配电网潮流计算的要点 (13)八、自我总结 (13)九、参考文献 (14)附录一 MATLAB的简介 (14)一、程序设计目的开式网络潮流计算：配电网的结构特点呈辐射状，在正常运行时是开环的；配电网的潮流计算采用的方法是前推回代法，本程序利用前推回代法的基本原理、收敛性。
(1)在电网规划阶段，通过潮流计算，合理规划电源容量及接入点，合理规划网架，选择无功补偿方案，满足规划水平年的大、小方式下潮流交换控制、调峰、调相、调压的要求。
MATLAB用于电机电磁计算的计算机编程
MATLAB用于电机电磁计算的计算机编程原文发表于《微特电机》2007年第6期MATLAB是一种功能极其强大的解释性语言，它提供的运算符及数据结构之丰富，图形功能之完美、功能箱之丰富，远超过FORTRAN、C等语言；而它的灵活简便又超过BASIC等语言，MATLAB语法限制不严格，特别适用于电机电磁计算的计算机程序，程序的自由度和可移植性好，还可通过Notebook将程序移到Word环境中享用。
以下就程序编制的一些心得与读者进行交流。
使用MATLAB编写电机设计程序，通常编写一个主程序文件和一个辅助文件夹（内放置需要的各个子程序，即专为设计服务的各种曲线和函数）。
通过path（路线）设置使其链接能通过主程序方便加以引用。
编写主程序的过程可以边输入、边验算，及时检验程序。
MATLAB没有FORTRAN、C等语言等要求的严格的语法规则，因此没有多少计算机语言基础的人都能很快掌握和随心所欲地使用。
1 曲线和图表的处理电机中的曲线和图表可能有两个来源：由复杂的解析式计算而来和根据试验——经验数据绘制。
对于由解析式计算来的曲线，一部分可以将其还原成原始公式直接代入。
对找不到原公式的曲线和图表，MATLAB提供了以下十分方便的方法：1.1 一元插值函数的应用MATLAB具有最强大和最完善的矩阵功能，在程序编制中可以直接利用MATLAB提供的一元插值函数准确查取曲线和图表数据。
现举例示范如下：例磁路计算中导磁材料的磁化曲线的查取编制一函数Hts=aw(Bts,FE1)，计算出电枢齿部磁密Bts 后，便自动查出对应单位磁路长度的励磁安匝，继续后面的计算，其中FE1是材料代号，需要采取的步骤是：（1）建立新文件“aw”，编制函数：function y=aw(xi,xj)，其中xi和xj对应于计算中的Bts和FE1；注：文件中编入了三种材料：21—DW540；315—DW315 ；235—Q235（欲直接用牌号代入应加引号，参见MATLAB文字计算部分）；下文中括号中为方便读者阅读的解释和说明，不能用此格式代入程序。
(完整word版)Matlab解决电路问题
如下图所示的电桥电路, 其中I1是16V 的电压源, I2是1A 的电流源,R1为8 , 电桥的四个臂分别为R2, R3, R4, R5电阻值如图所示, 求流过R4的电流I 的大小？解法一: 利用戴维南定理进行求解:解题思路：将A.B 两点断开, 求A.B 两点之间的等效电阻与等效电压, 等效之后的图形 如下图所示:I=? ABAB其中R6是等效电阻, I3是等效电压。
①求解等效电阻:求解等效电阻时把所有的电流源开路, 电压源短路, 得到如下所示的电路：AB则AB两端的电阻值即等效电阻R6=(R2+R3)//R1+R5②求解等效电压可以利用叠加法求解AB 两端的电压值, 先不看电压源（即电压源相当于短路）, 计算电流源对AB 两端的电压值, 再不看电流源（即电流源相当于断路）, 再计算AB 两端的电压值, 然后将俩种情况下的电压值叠加即得到AB 两端的等效电压。
不看电压源的电路图如下:则UCB+I2*R5+I2*（R1+R2）//R3=0 可以得到:UCB =-[I2*R5+I2*（R1+R2）//R3]U AB1 =-I2*R5-I2*3213)21(R R R R R R ++•++I2*R2*3213R R R R ++不看电流源的电路图如下:ABC很容易的知道AB 两端的电压值为:U AB2=321)32(*1R R R R R I +++所以UAB=UAB1+UAB2则经过戴维南等效之后的电路图如下:可以很简单的求解出II=64R R U AB+ABMatlab求解程序如下:(程序代码如下)R1=8;R2=4;R3=20;R4=3;R5=3;I1=16;I2=1;R6=R5+(R2+R3)*R1/(R1+R2+R3);UAB1=-I2*R5-I2*(R1+R2)*R3/(R1+R2+R3)+I2*R2*R3/(R1+R2+R3); UAB2=I1*(R2+R3)/(R1+R2+R3);UAB=UAB1+UAB2;I=UAB/(R4+R6);解法二: 运用叠加定理直接求解①先考虑电压源对AB两点的电流影响, 此时不看电流源, 电流源相当于断路, 电路图如下：根据电路图, 容易知道: AB 之间的电流I1 为I 1=543232)54//()32(11R R R R R R R R R R R I ++++•+++②再考虑电流源对AB 端电流源的影响, 此时不看电压源, 即将电压源短路, 电路图如下所示:根据电路图, 分析容易知道: 可以根据三角形与Y 形电路之间的转换, 将三角形电阻ACD 转换为Y 形电阻, 公式为：ABI 1BCD形电阻之和相邻电阻的乘积形电阻∆∆Y转换之后的电路图如下:可以得到:R12=32121R R R R R ++•R13=32131R R R R R ++•由于是电流源, 电流一定, 可以忽略与电流源串联的电阻R23 所以I 2=-I2*541312513R R R R R R ++++综上知道:I=I 1+I 2Matlab 求解程序如下: (程序代码如下) R1=8 R2=4;I 2R3=20; R4=3; R5=3; I1=16; I2=1;i1=[(R2+R3)/(R2+R3+R4+R5)]*I1/[R1+(R2+R3)*(R4+R5)/(R2+R3+R4+R5)];R12=R1*R2/(R1+R2+R3); R13=R1*R3/(R1+R2+R3);i2=-I2*(R13+R5)/(R12+R13+R4+R5); I=i1+i2解法三: 利用回路电流法进行求解 实验电路图如下:将无伴电流源的支路作为一个回路电流, 可以有电路图结合回路电i1i2流法列出如下方程:i1=I2I*（R2+R3+R4+R5）+i1*(R3+R5)-i2*（R2+R3）=0 -I*(R2+R3)-i1*R3+i2*(R1+R2+R3)=I1解方程可以很容易解的I 的值。
完整word版,BOOST电路设计及matlab仿真
Boost升压电路及MATLAB仿真一、设计要求1.输入电压（VIN）:12V2.输出电压（VO）：18V3.输出电流（IN）:5A4.电压纹波：0.1V5.开关频率设置为50KHz需设计一个闭环控制电路，输入电压在10—14V或负载电流在2—5A范围变化时，稳态输出能够保持在18V 。
根据设计要求很显然是要设计一个升压电路即Boost电路。
Boost电路又称为升压型电路,是一种开关直流升压电路，它可以是输出电压比输入电压高。
其工作过程包括电路启动时的瞬态工作过程和电路稳定后的稳态工作过程。
二、主电路设计图1主电路2.1 Boost电路的工作原理Boost升压电路电感的作用：是将电能和磁场能相互转换的能量转换器件，当MOS开关管闭合后，电感将电能转换为磁场能储存起来，当MOS断开后电感将储存的磁场能转换为电场能，且这个能量在和输入电源电压叠加后通过二极管和电容的滤波后得到平滑的直流电压提供给负载，由于这个电压是输入电源电压和电感的磁场能转换为电能的叠加后形成的，所以输出电压高于输入电压，既升压过程的完成。
Boost升压电路的肖特基二极管主要起隔离作用，即在MOS开关管闭合时，肖特基二极管的正极电压比负极的电压低，此时二极管反向截止，使此电感的储能过程不影响输出端电容对负载的正常供电；因在MOS管断开时，两种叠加后的能量通过二极向负载供电，此时二极管正向导通，要求其正向压降越小越好，尽量使更多的能量供给到负载端。
闭合开关会引起通过电感的电流增加。
打开开关会促使电流通过二极管流向输出电容因储存来自电感的电流，多个开关周期以后输出电容的电压升高，结果输出电压高于输入电压。
接下来分两部分对Boost电路作具体介绍即充电过程和放电过程。
充电过程在充电过程中，开关闭合（三极管导通），等效电路如图二，开关（三极管）处用导线代替。
这时，输入电压流过电感。
二极管防止电容对地放电。
由于输入是直流电，所以电感上的电流以一定的比率线性增加，这个比率跟电感大小有关。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
%电机设计程序clear allformat short em1=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.875eta=0.855 ;%5.功率因数cos(phi) =0.84,按照技术条件的规定cos(phi)=0.84phi=acos(0.84);cos(phi);%6.极对数p=2p=2;%7.定转子槽数：每极每相槽数取整数。
参考类似规格电机取q1=3,则Z1=2m1pq1，再查表10-8选Z2=32，并采用转子斜槽。
q1=3;Z1=2*m1*p*q1Z2=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.931P1=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.0004lef=lt + 2*deltaD2=Di1-2*delta%转子内径先按转轴直径决定（以后再校验转子轭部磁密）Di2=0.048 ;%11.极距tautau =pi*Di1/(2*p)%12.定子齿距t1t1=(pi*Di1/Z1)%转子齿距t2t2=(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.每相串联导体数Nphi1Nphi1=Ns1*Z1/(m1*a1)%每相串联匝数N1N1=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-h11r21=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层压板，槽绝缘占面积Aideltai=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%每相有效串联导体数NefNef=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.5mmb02=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/JB1AR11=(a+b)/2*(h02+h12+h22)*10^6AR1=input('输入端环面积AR1=')%其中端环电密JR1=0.6,JB1=2.1A/mm^2.按照工艺要求有所需面积确定端环内径及厚度如图10-25b,得端环面积AE=40*10^(-6)%(二)磁路计算%22KE1=0.918;epsilonl1=1-0.918;E1=(1-epsilonl1)*UN%23Ks(1)=1.2;KNm(1)=1.095;precision=1;while precision >0.01for n=1:3KNm(n)=input('输入波形系数KNm(n)=')alphap(n)=input('输入极弧系数alphap(n)=')phi(n)=E1/(4*KNm(n)*Kdp1*f*N1)%24At1=KFe*lt*bt1*Zp1At2=KFe*lt*bt213*Zp2%25hs1=h21+h01+h11+r21;hs2=h02+h12+h22;hj11=(D1-Di1)/2-hs1+r21/3hj21=(D2-Di2)/2-hs2Aj1=KFe*lt*hj11Aj2=KFe*lt*hj21%26Adelta=tau*lef%27Fs(n)=input('输入Fs(n)=')%28Bdelta(n)=Fs(n)*phi(n)/Adelta%29Bt1=Fs(n)*phi(n)/At1%30Bt2=Fs(n)*phi(n)/At2%31Ht1(n)=input('输入Ht1(n)=');Ht2(n)=input('输入Ht2(n)=');%32Kdelta1=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*Kdelta2deltaef=Kdelta*delta%33Lt1=(h11+h21)+1/3*r21Lt2=(h12+h22)%34Lj11=pi*(D1-hj11)/(4*p)Lj21=pi*(Di2+hj21)/(4*p)%35mu0=4*pi*10^(-7);Fdelta(n)=Kdelta*Bdelta(n)*delta/mu0%36Ft1(n)=Ht1(n)*Lt1*100Ft2(n)=Ht2(n)*Lt2*100%37Ks(2*n)=(Fdelta(n)+Ft1(n)+Ft2(n))/Fdelta(n)Ks(2*n+1)=Ks(2*n)-(Ks(2*n)-Ks(2*n-1))/3precision=abs((Ks(2*n-1)-Ks(2*n))/Ks(2*n))endendKs=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);%迭代的出来的结论%38Bj1=0.5*phi/Aj1%39Bj2=0.5*phi/Aj2%40Hj1=input('输入Hj1=')Hj2=input('输入Hj2=')%41hj11/tauCj1=input('根据Bj1和hj11/tau输入Cj1=')Fj1=Cj1*Hj1*Lj11*10^2hj21/tauCj2=input('根据Bj2和hj21/tau输入Cj2=')Fj2=Cj2*Hj2*Lj21*10^2%42F0=Fdelta+Ft1+Ft2+Fj1+Fj2%43Im=2*p*F0/(0.9*m1*N1*Kdp1)%44Ims=Im/IKW%45Xms=(4*f* mu0*m1/pi)*((N1*Kdp1)^2/(Ks*p))*lef*(tau/deltaef)Xmss=Xms*IKW/UN%（三)参数计算%46beta=(2*8+1*7)/(3*9)tauy=pi*(Di1+2*(h01+h11)+h21+r21)*beta/(2*p)d1=0.015;lB=lt+2*d1K0=1.2;l0=lB+K0*tauy%47lE=2*d1+K0*tauy%48Cx=4*pi*f*mu0*(N1*Kdp1)^2*lef*PN/(m1*p*UN^2)%49KU1=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%50Xs1s=2*m1*p/(Z1*Kdp1^2)*(lt/lef)*lambdas1*Cx%51sigmas=0.0129;Xdelta1s=(m1/(pi^2)*(tau/deltaef)*(sigmas/((Kdp1^2)*Ks)))*Cx %52XE1s=0.47*(lE-0.64*tauy)*Cx/(lef*Kdp1^2)%53Xsigma1s=Xs1s+Xdelta1s+XE1s%54%55lambdaU2=h02/b02hb2=h22/b22,b1b2=b12/b22lambdaL=input('输入lambdaL=');lambdaL2=2*h12/(b02+b12)+lambdaLlambdas2=lambdaU2+lambdaL2%56Xs2s=2*m1*p/Z2*(lt/lef)*lambdas2%57Z22p=Z2/(2*p)sigmaR=input('根据Z22p输入sigmaR=')Xdelta2s=m1*tau*sigmaR/(pi^2*deltaef*Ks)%58DR=0.107;XE2s=(0.2523*Z2*DR/(2*p*lef*2*p))*(2*m1*p/Z2)%59Xsks=0.5*(bsk/t2)^2*Xdelta2s%60Xsigma2s=(Xs2s+Xdelta2s+Xsks+XE2s)*Cx%61Xsigmas=Xsigma1s+Xsigma2s%62rhow=0.0217*10^(-6);R1=rhow*(2*N1*l0)/(Nt1*Ac11*a1)%63R1s=R1*IKW/UN%64C=1.05;rhoCu=8.9*10^3;GCu=C*l0*Ns1*Z1*Ac11*Nt1*rhoCurhoFe=7.8*10^3;GFe=KFe*lt*(D1+delta)^2*rhoFe%65 AR1代替ARKB=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/UNRR1=rhow1*Z2*DR/(2*pi*p^2*AR1*10^(-6))*(4*m1*(N1*Kdp1)^2/Z2)RRs=RR1*IKW/UNR2s=R21*IKW/UN%66eta(1)=0.855;I1Ps(1)=1/eta(1);precisioneta=1;while precisioneta>0.005for m=1:2I1Ps(m)=1/eta(2*m-1)%67sigma1=1+Xsigma1s/XmssIXs(m)=sigma1*Xsigmas*(I1Ps(m))^2*(1+(sigma1*Xsigmas*I1Ps(m))^2)%68I1Qs(m)=Ims+IXs(m)%69epsilonL(m)=I1Ps(m)*R1s+I1Qs(m)*Xsigma1sKE=1-(I1Ps(m)*R1s+I1Qs(m)*Xsigma1s)%70epsilon0=Ims*Xsigma1s1-epsilon0%71Bt10(m)=(1-epsilon0)/(1-epsilonL(m))*Bt1Ht10(m)=input('输入Ht10(m)=')Bt20(m)=(1-epsilon0)/(1-epsilonL(m))*Bt2Ht20(m)=input('输入Ht20(m)=')%73Bj10(m)=(1-epsilon0)/(1-epsilonL(m))*Bj1Hj10(m)=input('输入Hj10(m)=')%74Bj20(m)=(1-epsilon0)/(1-epsilonL(m))*Bj2Hj20(m)=input('输入Hj20(m)=')%75Bdelta0(m)=(1-epsilon0)/(1-epsilonL(m))*Bdelta%76Ft10(m)=Ht10(m)*Lt1*100%77Ft20(m)=Ht20(m)*Lt2*100%78hj11/tau,Bj10Cj1(m)=input('输入Cj1(m)=')Fj10(m)=Cj1(m)*Hj10(m)*Lj11*100%79hj21/tau,Bj20Cj2(m)=input('输入Cj2(m)=')Fj20(m)=Cj2(m)*Hj20(m)*Lj21*100%80Fdelta0(m)=Kdelta*delta*Bdelta0(m)/mu0%81F00(m)=Fdelta0(m)+Fj20(m)+Fj10(m)+Ft20(m)+Ft10(m) %82Im0=2*p*F00(m)/(0.9*m1*N1*Kdp1)%(四)工作性能计算%83I1s(m)=(I1Ps(m)^2+I1Qs(m)^2)^(0.5)I1(m)=I1s(m)*IKW%84J1(m)=(I1(m)/(a1*Nt1*Ac11))/10^6%85A1(m)=m1*Nphi1*I1(m)/(pi*Di1)%86I2s(m)=(I1Ps(m)^2+IXs(m)^2)^(0.5)I2(m)=I2s(m)*IKW*m1*Nphi1*Kdp1/Z2IR(m)=I2(m)*Z2/(2*pi*p)%87JB(m)=I2(m)/(AB*10^6)JR(m)=IR(m)/(AR1)%88pCu1s(m)=I1s(m)^2*R1spCu1(m)=pCu1s(m)*PN%89pA12s(m)=I2s(m)^2*R2spA12(m)=pA12s(m)*PN%90pss=0.02;ps=pss*PN%91pfw=(3/p)^2*(D1)^4*10^4pfws=pfw/PN%92Bj10phej(m)=input('输入phej(m)=')k2=2;Gj=4*p*Aj1*Lj11*rhoFepFej(m)=k2*phej(m)*GjBt10phet(m)=input('输入phet(m)=')Gt=2*p*At1*Lt1*rhoFek1=2;pFet(m)=k1*phet(m)*GtpFe(m)=pFej(m)+pFet(m)pFes(m)=pFe(m)/PN%93sigmaps(m)=pCu1s(m)+pA12s(m)+pss+pfws+pFes(m) %94PN1s(m)=1+sigmaps(m)%95eta(2*m)=1-sigmaps(m)/PN1s(m)eta(2*m+1)=eta(2*m)+(eta(2*m)-eta(2*m-1))/5precisioneta=abs((eta(2*m)-eta(2*m-1))/eta(2*m))endend%96I1s=I1s(2);I1Ps=I1Ps(2);phi=acos(I1Ps/I1s);cos(phi)%97pFers=((1-1/2)*pFej+(1-1/2.5)*pFet)*10^(-3)/PN*1000sN=pA12s/(1+pA12s+pFers+pss+pfws)%98nN=60*f/p*(1-sN)%99Tms=(1-sN)/(2*(R1s+(R1s^2+Xsigmas^2)^(0.5)))%(五)启动性能%100K(1)=3;%101precisionist=1;K(1)=3;Ist(1)=K(1)*Tms*IKWwhile precisionist>0.03for k=1:3Fst(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)%102KZ(k)=input('输入启动漏抗饱和KZ（k）=')%103cs1(k)=(t1-b01)*(1-KZ(k))cs2(k)=(t2-b02)*(1-KZ(k))%105deltalambdaU1(k)=(h01+0.58*h11)/b01*(cs1(k)/(cs1(k)+1.5*b01)) lambdas1st(k)=KU1*(lambdaU1-deltalambdaU1(k))+KL1*lambdaL1 %106Xs1sts(k)=lambdas1st(k)/lambdas1*Xs1s%107Xdelta1sts(k)=KZ(k)*Xdelta1s%108Xsigma1sts(k)=Xs1sts(k)+Xdelta1sts(k)+XE1s%109hB=h22;xi=1.987*10^(-3)*hB*(f/rhow1)^(0.5)%110b1b2=b22/b12;KF(k)=input('输入电阻增加系数KF(k)=')Kx(k)=input('输入电抗增加系数Kx(k)=')%111deltalambdaU2(k)=h02/b02*(cs2(k)/(cs2(k)+b02))lambdas2st(k)=(lambdaU2-deltalambdaU2(k))+Kx(k)*lambdaL2%112Xs2sts(k)=lambdas2st(k)*Xs2s/lambdas2%113Xdelta2sts(k)=KZ(k)*Xdelta2s%114Xksts(k)=KZ(k)*Xsks%115Xsigma2sts(k)=(Xs2sts(k)+Xdelta2sts(k)+Xksts(k)+XE2s)*Cx%116Xsigmasts(k)=Xsigma1sts(k)+Xsigma2sts(k)%117R2sts(k)=KF(k)*RBs+RRs%118Rsts(k)=R1s+R2sts(k)Zsts(k)=(Rsts(k)^2+(Xsigmasts(k))^2)^(0.5)%120Ist(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)=')endend%121ist=Ist(6)/I1(2)%122Tsts=R2sts(3)/Zsts(3)^2*(1-sN)。