哈工大机械设计电算实验(matlab2010)
哈工大 机电产品现代设计方法实验报告
1 实验目的（1）掌握典型机电产品多学科协同优化设计软件环境组成，包括建模软件、分析软件、协同平台；（2）自主设计产品模型、分析过程、优化目标；(3) 对得到的优化结果进行定性分析，解释结果的合理性，编写上机实验报告。
2 实验内容(1) 轴或负载台的有限元分析(2) 基于Adams的运动学分析与仿真3实验相关情况介绍（包含使用软件或实验设备等情况）网络协同设计环境，如图1所示：包括产品CAD建模、有限元分析FEM、动力学仿真ADAMS和控制仿真MATLAB。
计算机网络硬件环境和相应软件环境。
图形工作站和路由器，安装协同设计仿真软件。
型图1 协同设计仿真平台组成典型机电产品协同设计仿真工作流程如下图2所示。
1）利用CAD建模工具，建立产品模型；2）利用ADAMS建立产品运动学模型；3）根据CAD和ADAMS传过来的结构模型和边界条件分析零件应力场和应变场；4）用ADAMS分析得到的运动参数（位移、速度）。
图2 协同设计仿真平台组成SysML语言是UML语言(Unified Modeling Language，统一建模语言，一种面向对象的标准建模语言，用于软件系统的可视化建模)在系统工程应用领域的延续和扩展，是近年提出的用于系统体系结构设计的多用途建模语言，用于对由软硬件、数据和人综合而成的复杂系统的集成体系结构进行可视化的说明、分析、设计及校验。
在这里我们绘制参数图如下。
在下面的参数图中，我们确定了系统中各部件的相互约束情况。
图４产品初步结构与ＳｙｓＭＬ图4实验结果（含操作过程说明、结果记录及分析和实验总结等，可附页）（一）底座转台关键件有限元分析：1，在CAD中打开零件的三维模型图，导出为IGES格式模型文件（*.igs）,在Ansys中运行file->import->IGSE...导入该模型; 或者按照以下步骤创建零件模型。
运行Preprocessor->Modeling->V olumns->Cylinder->Solide Cylinder，弹出如下对话框，在对话框中输入相应数值，。
哈尔滨工业大学机械设计大作业_带传动电算
H a r b i n I n s t i t u t e o f T e c h n o l o g y上机电算说明书课程名称：机械设计电算题目：普通V带传动院系：机电工程学院班级：设计者：学号：指导教师：设计时间：2015.11.11-2015.12.1哈尔滨工业大学目录一、普通V带传动的内容 (1)二、变量标识符 (1)三、程序框图 (2)四、V带设计C程序 (3)五、程序运行截图 (10)参考文献 (11)一、普通V带传动的内容给定原始数据：传递的功率P，小带轮转速n1 传动比i及工作条件设计内容：带型号，基准长度Ld，根数Z，传动中心距a，带轮基准直径dd1、dd2，带轮轮缘宽度B，初拉力F0和压轴力Q。
二、变量标识符为了使程序具有较好的可读性易用性，应采用统一的变量标识符，如表1所示。
表1变量标识符表。
表1 变量标识符表三、程序框图四、V带设计c程序#include<stdio.h>#include<math.h>#include<string.h>#define PAI 3.14int TYPE_TO_NUM(char type); /*将输入的字符（不论大小写）转换为数字方便使用*/float Select_Ki(float i); /*查表选择Ki的程序*/float Select_KA(float H,int GZJ,int YDJ); /*查表选择KA的程序*/float Select_KL(float Ld,int TYPE); /*查表选择KL的程序*/float Select_M(int TYPE); /*查表选择m的程序*/float Select_dd1(int TYPE); /*查表选择小轮基准直径dd1的程序*/ float Select_dd2(int dd1,float i); /*查表选择大轮直径dd2的程序*/ float Compute_P0(float V,int TYPE,float dd1); /*计算P0的程序*/float Compute_DIFP0(int TYPE,float Ki,int n1); /*计算DIFP0的程序*/float Compute_VMAX(int TYPE); /*计算VMAX的程序*/float Compute_KALF(float ALF1); /*计算KALF的程序*/float Compute_B(int TYPE,int z); /*计算带宽B的程序*/float* Compute_LAK(float dd1,float dd2,int TYPE); /*计算Ld，a，KL的程序*/ main(){floatP,H,i,n1,KA,Ki,dd1,dd2,V,P0,DIFP0,Pd,VMAX,*LAK,m,Ld,KALF,a,KL,z,F0,ALF1, Q,B;int YDJ,GZJ,TYPE,ANS;char type,ans;printf(" V带传动设计程序\n");printf(" 程序设计人:×××\n 班号:123456678\n 学号:1234567896\n"); START: printf("请输入原始参数:\n");printf("传递功率P(KW):");scanf("%f",&P);printf("小带轮转速n1(r/min)：");scanf("%f",&n1);printf("传动比i:");scanf("%f",&i);printf("每天工作时间H(h):");scanf("%f",&H);printf("原动机类型(1或2):");scanf("%d",&YDJ);printf("工作机载荷类型:\n1,载荷平稳\n2,载荷变动较小\n3,载荷变动较大\n4,载荷变动很大\n");scanf("%d",&GZJ);DX: printf("带型:");scanf(" %c",&type);TYPE=TYPE_TO_NUM(type);KA= Select_KA(H,GZJ,YDJ);Pd=KA*P;VMAX=Compute_VMAX(TYPE);DD1:dd1=Select_dd1(TYPE);V=PAI*dd1*n1/60000;while(V>VMAX){printf("所选小轮基准直径过小!请重新选择!\n");dd1=Select_dd1(TYPE);V=PAI*dd1*n1/60000;}dd2=Select_dd2(dd1,i);P0=Compute_P0(V,TYPE,dd1);Ki=Select_Ki(i);DIFP0=Compute_DIFP0(TYPE,Ki,n1);LAK=Compute_LAK(dd1,dd2,TYPE);Ld=LAK[0];a=LAK[1];KL=LAK[2];ALF1=180-57.3*(dd2-dd1)/a;KALF=Compute_KALF(ALF1);z=ceil(Pd/(KALF*KL*(P0+DIFP0)));if(z>=10)goto DX;m=Select_M(TYPE);F0=500*Pd*(2.5-KALF)/(z*V*KALF)+m*V*V;Q=2*z*F0*sin(ALF1*PAI/360);B= Compute_B(TYPE,z);printf(" 计算结果\n");printf("小轮直径:%f mm\n大轮直径:%f mm\n中心距:%f mm\n带长:%f mm\n带宽:%f mm\n带的根数:%f\n初拉力:%f N\n轴压力:%f N\n",dd1,dd2,a,Ld,B,z,F0,Q); X: printf("计算下一组带轮直径按'Y',重新开始按'S',结束按'N':");scanf(" %c",&ans);ANS=TYPE_TO_NUM(ans);if(ANS==6)goto DD1;else if(ANS==7)goto START;else if(ANS==8)printf("程序结束");else{printf("输入错误，请重新输入");goto X;}}float Select_KA(float H,int GZJ,int YDJ){float ka1[4][3]={{1.0,1.1,1.2},{1.1,1.2,1.3},{1.2,1.3,1.4},{1.3,1.4,1.5}};float ka2[4][3]={{1.1,1.2,1.3},{1.2,1.3,1.4},{1.4,1.5,1.6},{1.5,1.6,1.8}};float KA;if(YDJ==1){if(H<10)KA=ka1[GZJ-1][0];else if(H>16)KA=ka1[GZJ-1][2];elseKA=ka1[GZJ-1][1];}if(YDJ==2){if(H<10)KA=ka2[GZJ-1][0];else if(H>16)KA=ka2[GZJ-1][2];elseKA=ka2[GZJ-1][1];}return KA;}float Select_Ki(float i){float m;floatKi[10]={1.0000,1.0136,1.0276,1.0419,1.0567,1.0719,1.0875,1.1036,1.1202,1.1 373};float I[10]={1.00,1.02,1.05,1.09,1.13,1.19,1.25,1.35,1.52,2.00};int j=9;for(j=9;j>=0;j--){m=i-I[j];if(m>=0)break;}return Ki[j];}int TYPE_TO_NUM(char type){int x;if(type<91)type+=32;switch(type){case 'z':x=0;break;case 'a':x=1;break;case 'b':x=2;break;case 'c':x=3;break;case 'd':x=4;break;case 'e':x=5;break;case 'y':x=6;break;case 's':x=7;break;case 'n':x=8;break;}return x;}float Compute_P0(float V,int TYPE,float dd1){float P0;float K1[6]={0.246,0.449,0.794,1.48,3.15,4.57};float K2[6]={7.44,19.02,50.6,143.2,507.3,951.5};float K3[6]={0.441e-4,0.765e-4,1.31e-4,2.34e-4,4.77e-4,7.06e-4};P0=(K1[TYPE]*pow(V,-0.09)-K2[TYPE]/dd1-K3[TYPE]*V*V)*V;return P0;}float Compute_DIFP0(int TYPE,float Ki,int n1){float Kb,DIFP0;float KB[6]={0.2925e-3,0.7725e-3,1.9875e-3,5.625e-3,19.95e-3,37.35e-3};Kb=KB[TYPE];DIFP0=Kb*n1*(1-1/Ki);return DIFP0;}float Select_dd1(int TYPE){int i;float dd1;float DD1[6][4]={{50,63,71,80},{75,90,100,125},{125,140,160,180},{200,250,315,400},{355,400,450,500},{500,560,630,710}};printf("该带型推荐的小带轮基准直径有(mm):\n");for(i=0;i<4;i++){printf("%.2f ",DD1[TYPE][i]);}printf("\n请选择小带轮直径dd1:");scanf("%f",&dd1);return dd1;}float Select_dd2(int dd1,float i){float DD2[]={50,56,63,71,80,90,100,112,125,140,150,160,180,200,224, 250,280,315,355,400,425,450,500,560,600,630,710,800};float dd2,temp,t;int j=0;temp=i*dd1;for(j=0;j<28;j++){t=temp-DD2[j];if(t<=0)break;}if((DD2[j]-temp)<(temp-DD2[j-1]))dd2=DD2[j];elsedd2=DD2[j-1];return dd2;}float* Compute_LAK(float dd1,float dd2,int TYPE){int j=0;float a1,a2,a0,Ld0,temp,KALF,Ld,a,KL;float LAK[3];float kl[23][6]={{400,0.87},{450,0.89},{500,0.91},{560,0.94},{630,0.96,0.81},{710,0.99,0.83},{800,1.00,0.85,0.82},{900,1.03,0.8 7,0.84,0.83},{1000,1.06,0.89,0.86,0.86,0.83},{1120,1.08,0.91,0.88,0.88,0.86},{1250,1.10,0.93,0.90,0.91,0.89},{1400,1.14, 0.96,0.92,0.93,0.91},{1600,1.16,0.99,0.95,0.95,0.93},{1800,1.18,1.01,0.98,0.97,0.96},{2000,0,1.03,1.00,0.99,0.98},{2240,0,1.06,1 .03,1.02,1.00},{2500,0,1.09,1.05,1.04,1.03},{2800,0,1.11,1.07,1.07,1.06},{3150,0,1.13,1.09,1.09,1.08},{3550,0,1.17,1.13 ,1.12,1.11},{4000,0,1.19,1.15,1.15,1.14},{4500,0,0,1.18,1.18,1.17},{5000,0,0,0,1. 21,1.20}};a1=0.7*(dd1+dd2);a2=2*(dd1+dd2);printf("请输入初估中心距a0(范围为%.2fmm~%.2fmm):",a1,a2);scanf("%f",&a0);Ld0=2*a0+(dd1+dd2)*PAI/2+(dd2-dd1)*(dd2-dd1)/a0/4;for(j=0;j<23;j++){temp=Ld0-kl[j][0];if(temp<=0)break;}if((Ld0-kl[j-1][0])>(kl[j][0]-Ld0)){ Ld=kl[j][0];KL=kl[j][TYPE+1];}else{ Ld=kl[j-1][0];KL=kl[j-1][TYPE+1];}a=a0+(Ld-Ld0)/2;LAK[0]=Ld;LAK[1]=a;LAK[2]=KL;return LAK;}float Compute_VMAX(int TYPE){float VMAX;switch(TYPE){case 0:case 1:case 2:case 3:VMAX=25;break;case 4:case 5:VMAX=30;break;}return VMAX;}float Select_M(int TYPE){float m;float M[6]={0.06,0.1,0.17,0.3,0.6,0.9};m=M[TYPE];return m;}float Compute_B(int TYPE,int z){float f[6]={8,10,12.5,17,23,29};float e[6]={12,15,19,25.5,37,44.5};float B;B=(z-1)*e[TYPE]+2*f[TYPE];return B;}float Compute_KALF(float ALF1){floatkalf[2][14]={220,210,200,190,180,170,160,150,140,130,120,110,100,90,1.20,1. 15,1.10,1.05,1.00,0.98,0.95,0.92,0.89,0.86,0.82,0.78,0.73,0.68},KALF;int i;for(i=0;i<14;i++){if(ALF1<kalf[0][i]&&ALF1>kalf[0][i+1]){KALF=kalf[1][i]+(kalf[1][i]-kalf[1][i+1])*(ALF1-kalf[0][i])/10;break;}}return KALF;}五、程序运行截图电动机功率4kw,,转速960r/min,传动比2，工作时间24h。
机械设计matlab课程设计
机械设计matlab课程设计一、课程目标知识目标：1. 理解并掌握机械设计中常用的MATLAB命令和功能；2. 学会运用MATLAB进行机械系统建模、仿真和优化；3. 掌握运用MATLAB解决机械设计中的实际问题的方法和步骤。
技能目标：1. 能够运用MATLAB进行数据分析和处理，提高数据处理能力；2. 能够熟练使用MATLAB进行机械设计绘图，提高绘图技能；3. 能够独立完成机械设计相关MATLAB程序编写和调试，提升编程能力。
情感态度价值观目标：1. 培养学生对机械设计MATLAB课程的学习兴趣，激发学习热情；2. 培养学生严谨的科学态度，强调实际操作与理论知识的相结合；3. 增强学生的团队协作意识，培养合作解决问题的能力。
课程性质：本课程为实践性较强的学科，结合机械设计与MATLAB软件应用，注重培养学生的实际操作能力和解决问题的能力。
学生特点：学生具备一定的机械设计基础知识，对MATLAB软件有一定了解，但实际操作能力有待提高。
教学要求：结合课程性质和学生特点，采用任务驱动、案例教学等方法，使学生在实践中掌握知识，提高技能，培养情感态度价值观。
通过分解课程目标为具体学习成果，为教学设计和评估提供明确依据。
二、教学内容1. MATLAB基础命令与操作：介绍MATLAB的基本命令、数据类型、矩阵运算等，为学生后续学习打下基础。
（对应教材第一章）2. MATLAB绘图与可视化：讲解MATLAB绘图功能，如二维、三维图形绘制，让学生能够将机械设计结果可视化。
（对应教材第二章）3. 机械系统建模与仿真：结合教材第三章内容，教授学生如何使用MATLAB 进行机械系统的建模、仿真和结果分析。
4. 机械设计优化方法：介绍教材第四章内容，使学生掌握运用MATLAB进行机械设计优化的基本方法。
5. MATLAB在机械设计中的应用案例：分析教材第五章案例，让学生了解MATLAB在机械设计领域的实际应用。
6. 实践操作与项目实战：组织学生进行实践操作，结合教材内容开展项目实战，巩固所学知识，提高实际操作能力。
哈工大《电子系统》实验报告模板
H a r b i n I n s t i t u t e o f T e c h n o l o g y《电子系统》实验报告院系：电信学院班级：设计者：学号：指导教师：孙思博实验一连续波雷达测速实验一、实验目的：1、掌握雷达测速原理。
2、了解连续波雷达测速实验仪器原理及其使用。
3、使用Matlab对实验数据进行分析，得到回波多普勒频率和目标速度。
二、实验原理：1、多普勒测速原理：由于运动目标相对辐射源的运动而引起发射信号的中心频率发生多普勒频移的现象称为多普勒效应。
目标运动方向的不同决定了多普勒频移的正负。
（如图1所示）图1.多普勒效应假设发射的是重复频率为错误！未找到引用源。
的脉冲串，雷达发射信号的波长为错误！未找到引用源。
时，设目标的速度为错误！未找到引用源。
，多普勒频率为错误！未找到引用源。
,以目标接近雷达为例，错误！未找到引用源。
错误！未找到引用源。
为接收脉冲串频率新频率错误！未找到引用源。
错误！未找到引用源。
为雷达发射信号的载频则：错误！未找到引用源。
，当|错误！未找到引用源。
|<<c时（1）2、多普勒信息的提取：在连续波工作状态时，利用相干检波器可以得到和错误！未找到引用源。
相关的一系列频谱分量，回波分量中的错误！未找到引用源。
、错误！未找到引用源。
、2错误！未找到引用源。
等高频分量被多普勒滤波器滤除，则最后获得就是多普勒分量，利用公式(1) 可以求得目标的速度。
本实验中发射波长为3cm,采样率是2048HZ。
三、实验仪器：实验装置如下：5402DSP测速传感器混频器连续波发射机传感器输出信号放大滤波AD 串行接口PC 机FFT图3-2 连续波雷达测速实验仪器原理框图图3. 测速雷达传感器三、 实验内容与步骤：1、 利用给定装置，使用一挡光板作为目标物体，移动该物体，则通过测速雷达传感器（如图3）能够获得回波数据，并被DSP 芯片采样，采样频率为2048HZ 。
2、 通过示波器观察波形，选择一高频干扰少的波形，利用软件获得其2048个数据，并存储在计算机中。
哈工大机械设计电算实验(matlab2010)
机械设计电算实验一：普通V带传动设计内容和任务1、普通V带传动设计内容给定原始数据：传递的功率P，小带轮转速n，传动比i及工作条件。
设计内容：带型号，基准长度Ld，根数Z，传动中心距a，带轮基准直径dd1、dd2，带轮轮缘宽度B，初拉力F0，和压轴力FQ。
2、CAD任务：（1）编制V带传动设计程序框图。
（2）编制V带传动设计原程序。
（3）按习题或作业中数据运行程序，要求对每一组数据各按三种V带型号计算，对每一种带型号选三种小带轮直径进行计算并输出所有结果。
二、变量标识符序号名称代号程序符号单位1 名义功率P p kwpd kw2 设计功率Pdkka3 工作情况系数KA4 小带轮转速nn r/min15 工作机载荷特性载荷类型6 原动机类型原动机类型7 工作年限a年b班8 带型号V带带型dd1 mm9 小带轮基准直径dd1dd2 mm10 大带论基准直径dd211 带速v v m/s12 最大带速vvmax m/smaxa0 mm13 预估中心距ald mm14 基准长度Ld15 实际中心距 a aa mmr 度16 小带轮包角α117 V带根数Z zpo kw18 基本额定功率P19 功率增量ΔPpp0 kwki20 传动比系数Ki21 长度系数KklLka22 包角系数Kαkb23 弯曲影响系数Kb24 初拉力Ff0 N25 压轴力Q fq N26 轮缘宽度 B bb mm三、程序框图开始输入原始数据P,n,i并选择电机计算设计功率由带型号确定最大带速，小带轮直径系，带的相关参数m，f，e 相关系数Kb确定带型号选择小带轮直径验算带速v<vmax?计算大带轮直径给定范围，预估中心距AA计算Ld并取标准值，计算实际中心距并圆整计算小带轮包角计算选择相关系数：Kα，Ki 。
插值计算P0，计算ΔP0计算带根数ZZ<10?计算相关参数结束下一型号带下一带轮直径B B四、源程序与其说明程序说明：本程序用Matlab2010b软件编制，主要针对机械设计大作业上的题型设计。
哈工大 电机学 MATLAB 仿真 实验报告
基于MATLAB的电机学计算机辅助分析与仿真实验报告班级：学号：姓名：完成时间：一、实验内容1.1单相变压器不同负载性质的相量图通过MATLAB 画出单相变压器带感性，阻性，容性三种不同性质负载的变压器向量图1.2感应电机的S T -曲线通过MATLAB 画出三相感应电动机的转矩转差率曲线二、实验要求2.1单相变压器不同负载性质的相量图根据给定的仿真实例画出负载相位角30,0,302-=j 三种情况下得向量图，观察电压大小与相位的关系，了解总结负载性质不同对向量图的影响2.2感应电机的S T -曲线根据给定的实例，画出3.1~3.1-=s 的S T -曲线，了解感应电机临界转差率的大小和稳定工作区间的大小，给出定性分析三、实验方法3.1单相变压器不同负载性质的相量图1.单相变压器不同负载性质的相量图 (1)先画出负载电压'2U 的相量；(2)根据负载的性质和阻抗角画出二次电流（规算值）的相量(3)在2U 上加上一个与电流方向相同的压降，其大小为二次电流规算值'2I 与二次漏电阻规算值'2R 之积；再加上一个超前电流方向︒90的压降，其大小为二次电流'2I 规算值与二次漏电抗规算值'2χ之积； (4)根据上一步结果连线，得出'2E ； (5)超前'2E 方向︒90画出m Φ；(6)根据励磁电阻与电抗的大小得出励磁阻抗角，并超前m Φ一个励磁阻抗角的大小得出m I 的方向；(7)根据平行四边形法则，做出'2I -与m I 的和，即为1I ； (8)根据'21E E =得出1E ，并得出1E -。
(9)在1E -上加上一个与电流方向相同的压降，其大小为一次电流1I 与一次漏电阻1R 之积；再加上一个超前电流方向︒90的压降，其大小为一次电流1I 与一次漏电抗1χ之积；(10) 根据上一步结果连线，得出1U ；3.2感应电机的S T -曲线实验采用matlab 对转矩转差率曲线进行仿真。
哈工大机械设计大作业三-V带传动电算编程-推荐下载
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术关，系电，力通根保1据护过生高管产中线工资0不艺料仅高试可中卷以资配解料置决试技吊卷术顶要是层求指配，机置对组不电在规气进范设行高备继中进电资行保料空护试载高卷与中问带资题负料2荷试2，下卷而高总且中体可资配保料置障试时各卷，类调需管控要路试在习验最题；大到对限位设度。备内在进来管行确路调保敷整机设使组过其高程在中1正资中常料，工试要况卷加下安强与全看过，25度并52工且22作尽护下可1关都能于可地管以缩路正小高常故中工障资作高料；中试对资卷于料连继试接电卷管保破口护坏处进范理行围高整，中核或资对者料定对试值某卷，些弯审异扁核常度与高固校中定对资盒图料位纸试置，卷.编工保写况护复进层杂行防设自腐备动跨与处接装理地置，线高尤弯中其曲资要半料避径试免标卷错高调误等试高，方中要案资求，料技编试术5写、卷交重电保底要气护。设设装管备备置线4高调、动敷中试电作设资高气，技料中课并术3试资件且、中卷料中拒管包试试调绝路含验卷试动敷线方技作设槽案术，技、以来术管及避架系免等统不多启必项动要方方高式案中，；资为对料解整试决套卷高启突中动然语过停文程机电中。气高因课中此件资，中料电管试力壁卷高薄电中、气资接设料口备试不进卷严行保等调护问试装题工置，作调合并试理且技利进术用行，管过要线关求敷运电设行力技高保术中护。资装线料置缆试做敷卷到设技准原术确则指灵：导活在。。分对对线于于盒调差处试动，过保当程护不中装同高置电中高压资中回料资路试料交卷试叉技卷时术调，问试应题技采，术用作是金为指属调发隔试电板人机进员一行，变隔需压开要器处在组理事在；前发同掌生一握内线图部槽纸故内资障，料时强、，电设需回备要路制进须造行同厂外时家部切出电断具源习高高题中中电资资源料料，试试线卷卷缆试切敷验除设报从完告而毕与采，相用要关高进技中行术资检资料查料试和，卷检并主测且要处了保理解护。现装场置设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
matlab课程设计任务书2010测
课程设计任务书发布人：梁丽发布时间：2010年12月19日1、上机安排：（19周周一、三下午5,6节B307）注：按时到实验室上机，不得迟到早退2、课程设计说明书要求：1)写清楚整个设计过程，包括系统数学模型的分析、建立，控制器的设计、不同参数之间的控制效果比较，仿真过程等。
2)字数5000字左右；3)统一封皮，设计任务书放在最前面（封皮之后）；4)按照教务网站上毕业设计说明书格式要求打印设计说明书；5)设计说明书如有雷同，两人都不及格；3、题目分配：一班：学号尾数1-12作第1题学号尾数13-23作第2题学号尾数24-33作第3题（必须到图书馆三楼借参考资料2，复印相关内容）学号尾数34-43作第4题另外有两个选做题，感兴趣的可选择内蒙古科技大学课程设计任务书6、完成5000字左右的课程设计报告（包括设计原理、设计过程及结果分析几部分），要求给出设计的模型图和仿真曲线图。
四、进度安排1、根据给定的参数或工程具体要求，收集和查阅资料（一天）2、Matlab/Simulink建模（一天）3、控制系统设计与优化（两天）4、编写课程设计说明书（一天）五、评分标准1）工作态度（占10%）；2）基本技能的掌握程度（占20%）；3）方案的设计是否可行和优化(40%)；4）课程设计技术设计书编写水平(占30%)。
分为优、良、中、合格、不合格五个等级。
六、建议参考资料1、李国勇，控制系统数字仿真与CAD[M]，北京：电子工业出版社，2003，92、王丹力，MATLAB控制系统设计仿真应用[M]，北京：中国电力出版社，2007，92、薛定宇，控制系统仿真与计算机辅助设计[M]，北京：机械工业出版社，2005，13、金以慧，过程控制[M]，北京：清华大学出版社，2003，6内蒙古科技大学课程设计任务书课程名称控制系统仿真设计题目基于模糊控制器的自动着舰控制系统指导教师梁丽时间2010.12.27~2010.12.31一、教学要求1、学会收集和查阅资料，学会针对指定控制系统建立数学模型的方法；2、学会使用Matlab/Simulink建模和仿真的方法；3、掌握控制器的设计方法，以及控制器参数整定和优化的方法。
实验二《机械工程控制基础》Matlab仿真实验报告单
红河学院工学院实验报告单《机械工程控制基础》Matlab仿真实验报告单课程名称：机械工程控制基础姓名：徐海洋日期：2010-11-26 成绩：年级专业：07机械学号：5 实验场地：文514实验序号：实验二实验名称：二阶线性系统时间响应与时域性能指标的Matlab仿真实验一、实验内容。
1、二阶线性系统216 416s s++单位脉冲响应、单位阶跃响应、单位正弦输入响应的Matlab 仿真。
2、二阶线性系统2424s s++单位阶跃响应的时域性能指标。
3、二阶线性系统236 1236s sξ++，当0.1,0.2,0.4,0.7,1.0,2.0ξ=时，单位阶跃响应的Matlab仿真。
二、实验目的。
三、相关Matlab仿真程序。
程序一：>> clearnum=[16];den=[1 4 16];Gs=tf(num,den)t=0:0.1:20;figure(1);impulse(Gs,t)xlabel('时间')ylabel('输出')title('二阶系统单位脉冲响应') figure(2);step(Gs,t)xlabel('时间')ylabel('输出')title('二阶系统单位阶跃响应') figure(3);r=sin(t);lsim(Gs,r,t)xlabel('时间')ylabel('输出')title('二阶系统单位正弦输入响应')Transfer function:16--------------s^2 + 4 s + 16程序二：>> clearnum=[4];den=[1 2 4];disp('二阶系统传递函数')Gs=tf(num,den)disp('无阻尼固有频率与阻尼比')[Wn Xita ]=damp(Gs)Wn=Wn(1)Xita= Xita(1)disp('二阶系统性能指标如下')tr=(pi-atan(sqrt(1-Xita^2)/ Xita))/Wn/sqrt(1-Xita^2) tp=pi/Wn/sqrt(1-Xita^2)ts= 3.5/Xita/WnMp= exp(-pi*Xita/sqrt(1-Xita^2))*100二阶系统传递函数Transfer function:4-------------s^2 + 2 s + 4无阻尼固有频率与阻尼比Wn =2.00002.0000Xita =0.50000.5000Wn =2.0000Xita =0.5000二阶系统性能指标如下tr =1.2092tp =1.8138ts =3.5000Mp =16.3034程序三：clearclcnum=1;Wn=6;xita1=0.1;den1=[1/Wn^2 2*xita1/Wn 1]; xita2=0.2;den2=[1/Wn^2 2*xita2/Wn 1]; xita3=0.4;den3=[1/Wn^2 2*xita3/Wn 1]; xita4=0.7;den4=[1/Wn^2 2*xita4/Wn 1]; xita5=1.0;den5=[1/Wn^2 2*xita5/Wn 1]; xita6=2.0;den6=[1/Wn^2 2*xita6/Wn 1];t=0:0.05:4;[y1 x t]=step(num,den1,t); [y2 x t]=step(num,den2,t);[y3 x t]=step(num,den3,t);[y4 x t]=step(num,den4,t);[y5 x t]=step(num,den5,t);[y6 x t]=step(num,den6,t);figure(1)plot(t,y1,'b-',t,y2,'r--',t,y3,'k',t,y4,'g:',t,y5,'m-.',t,y6,'b.');xlabel('t(s)'),ylabel('y(t)');title('阻尼比 xita=0.1, 0.2, 0.4, 0.7, 1.0, 2.0')legend('xita1=0.1','xita2=0.2','xita3=0.4','xita4=0.7','xita5=1','xita6=2') 四、相关Matlab仿真图形。
哈工大电力系统分析短路及潮流计算实验上机程序
fprintf(fid,'+%8.4fjn',im);
end
end
end end %%%============
fclose(fid);
%关闭文件
附：'node5.txt'
2、程序输出结果 节点导纳阵：
节点阻抗阵：
短路电流：
If= 0.0001+ 0.1082j 节点电压：
各支路电流：
for j=1:1:4
k=I(i,j);
re=real(k);
im=imag(k);
if(re~=0||im~=0)
fprintf(fid,'Branch%d--%d:I%d%d=',i,j,i,j);
fprintf(fid,'%8.4f',re);
if im<0
fprintf(fid,'%8.4fjn',im);
二、实验内容
1、编写数据输入、形成导纳阵程序
2、电力系统短路计算实用公式的计算程序及编制和调试
三、实验过程
1、程序代码及说明
%清屏 clear %读数据
fid=fopen('node5.txt','r');A=fscanf(fid,'%d',2);B=fscanf(fid,'%f',[5,A(2,1)]);fclose(fid);B=B';%求节点导纳矩阵
%===========================================
for j=1:1:A(1,1)-1
%===========================i~=j
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
哈工大机械设计电算实验(matlab2010)机械设计电算实验一：普通V带传动设计内容和任务1、普通V带传动设计内容给定原始数据：传递的功率P，小带轮转速n，传动比i及工作条件。
设计内容：带型号，基准长度Ld，根数Z，传动中心距a，带轮基准直径dd1、dd2，带轮轮缘宽度B，初拉力F0，和压轴力FQ。
2、CAD任务：（1）编制V带传动设计程序框图。
（2）编制V带传动设计原程序。
（3）按习题或作业中数据运行程序，要求对每一组数据各按三种V带型号计算，对每一种带型号选三种小带轮直径进行计算并输出所有结果。
二、变量标识符kka3 工作情况系数KAn r/min4 小带轮转速n15 工作机载荷特性载荷类型6 原动机类型原动机类型7 工作年限a年b班8 带型号V带带型dd1 mm9 小带轮基准直径dd1dd2 mm10 大带论基准直径dd211 带速v v m/svmax m/s12 最大带速vmaxa0 mm13 预估中心距ald mm14 基准长度Ld15 实际中心距 a aa mmr 度16 小带轮包角α117 V带根数Z z18 基本额定功率Ppo kwpp0 kw19 功率增量ΔPki20 传动比系数Ki三、程序框图四、源程序与其说明程序说明：本程序用Matlab2010b软件编制，主要针对机械设计大作业上的题型设计。
使用时只要打开m 文件，并点击运行，按照提示进行即可。
首先输入原始数据，然后根据自己的需要选择带型，中心距即可得到设计结果，无需再查找资料，方便高效,计算过程如有错误会进行提示，并返回到输入处进行改正。
而且该程序可以直接计算下一带轮直径或者计算下一带型，比较方便。
源程序如下(先复制到记事本，再新建一个m文件，粘贴)clear all;disp('欢迎使用本程序，请输入V带传动设计的原始数据');p=input('电动机工作功率(kw) P=');n=input('电动机满载转速(r/min) nm=');i=input('第一级传动比 i1=');a=input('请输入最短工作工作年限 a年b班 a=');b=input(' b=');disp('是否反复起动、正反转频繁或工作条件恶劣');ka1=input('是请输入1，否请输入0。
请输入：');disp('原动机类型:');disp('I类原动机包括普通笼型交流电动机，同步电动机，');disp(' 直流电动机（并激），n>=600r/min的内燃机')disp('II类原动机包括交流电动机(双笼型、滑环式、单相、大转差率)，');disp(' 直流电动机（复激、串激），单缸发动机，n<=600r/min的内燃机')d1=input('请选择原动机的类型，输入1或2。
请输入:');disp('载荷类型：');disp('1、载荷平稳；2、载荷变动小；3、载荷变动较大；4、载荷变动很大。
');e1=input('请选择载荷类型请输入:');c=[3 710 38 80;2.2 940 28 60;3 960 38 80;4 960 38 80];if p==3&&n==710d=38;e=80;disp('根据您的数据为您选择的电动机为 Y132M-8'); elseif p==2.2&&n==940d=28;e=60;disp('根据您的数据为您选择的电动机为 Y112M-6'); elseif p==3&&n==960d=38;e=80;disp('根据您的数据为您选择的电动机为 Y132S-6'); elseif p==4&&n==960d=38;e=80;disp('根据您的数据为您选择的电动机为 Y132M1-8'); elsedisp('欢迎继续使用本程序');endka2=[1.0 1.1 1.2;1.1 1.2 1.3;1.2 1.3 1.4;1.3 1.4 1.5]; ka3=[1.1 1.2 1.3;1.2 1.3 1.4;1.4 1.5 1.6;1.5 1.6 1.8]; if d1<=1k=ka2;elsek=ka3;endif ka1>=1kka=1.1*k(e1,b);elsekka=k(e1,b);endpd=kka*p;t2=1;while t2==1z=10;while z>=10disp('V带带型：1、Z型:dd1min=50mm；2、A型:dd1min=75mm；3、B型:dd1min=125mm'); x=input('请输入（1、2、3）来选择带型。
请选择：');x1=[0.246 0.449 0.794;7.44 19.02 50.6;0.441 0.765 1.31;0.2925 0.7725 1.9875]; k1=x1(1,x);k2=x1(2,x);k3=(10^(-4))*x1(3,x);kb=10^(-3)*x1(4,x);vmax=25;disp('优选带轮直径系列： 50 56 63 71 80 90 100')disp(' 112 125 140 150 160 180 200')disp(' 224 250 280 315 355 400 425')disp(' 450 500 560 600 630 710 800')dd1=input('请参考dd1min以及优选系列，输入小带轮基准直径(mm) dd1=');v=pi*dd1*n/60/1000;while v>=vmaxdisp('请减小小带轮基准直径dd1');dd1=input('请参考dd1min,重新输入小带轮基准直径(mm) dd1=');v=pi*dd1*n/60/1000;enddd2=round(i*dd1);a1=0.7*(dd1+dd2);a2=2*(dd1+dd2);disp('a1=');disp(a1);disp('a2=');disp(a2);a0=input('已知中心距a1<=a0<=a2,请输入a0的值（尽量为整数）。
a0=');while a0<a1||a0>a2disp('wrong!您输入的数值不在a1-a2之间，请重新输入');a0=input('请重新输入a0的值（尽量为整数）。
a0=');endld=2*a0+pi/2*(dd1+dd2)+(dd1-dd2)^2/4/a0;kk=[400 0.87 0 0;450 0.89 0 0;500 0.91 0 0;560 0.94 0 0;630 0.96...0.81 0;710 0.9 0.83 0;800 1.00 0.85 0.82;900 1.03 0.87 0.84;...1000 1.06 0.89 0.86;1120 1.08 0.91 0.88;1250 1.10 0.93 0.90;...1400 1.14 0.96 0.92;1600 1.16 0.99 0.95;1800 1.18 1.01 0.98;...2000 0 1.03 1.00;2240 0 1.06 1.03;2500 0 1.09 1.05;2800 0 ...1.11 1.07;3150 0 1.13 1.09;3550 0 1.17 1.13;4000 0 1.19 1.15;... 4500 0 0 1.18];for g=1:22if ld>=kk(g,1)&&ld<=kk(g+1,1)if (ld-kk(g,1))>=(kk(g+1,1)-ld)ld1=kk(g+1,1);elseld1=kk(g,1);endbreakendendkl=kk(g,x+1);aa=round(a0+(ld-ld1)/2);r=180-(dd2-dd1)/aa*57.3;rr=[220 1.20;210 1.15;200 1.10;190 1.05;180 1.00;170 0.98;160 0.95;... 150 0.92;140 0.89;130 0.86;120 0.82;110 0.78;100 0.73;90 0.68];for ii=1:13if r<=rr(ii,1)&&r>=rr(ii+1,1)ka=rr(ii+1,2)+(r-rr(ii+1,1))/10*(rr(ii,2)-rr(ii+1,2));breakendendp0=(k1*v^(-0.09)-k2/dd1-k3*v^2)*v;if i==1.8ki=1.1202;elseif i>=2ki=1.1372;endm=[0.06 0.1 0.17];pp0=kb*n*(1-1/ki);z=ceil(pd/(p0+pp0)/kl/ka);if z<10f0=500*pd/z/v*(2.5-ka)/ka+m(x)*v^2; fq=2*z*f0*sin(r/2/180*pi);bb=round((z-1)*15+20);disp('所选原动机种类：')if d1==1disp('I类原动机');elsedisp('II类原动机');enddisp('所选带型为：');if x==1disp('Z型V带');elseif x==2disp('A型V带');elsedisp('B型V带');enddisp('小带轮直径(mm) dd1=')disp(dd1);disp('大带轮直径(mm) dd2=')disp(dd2);disp('实际中心距(mm) a0=');disp(aa);disp('小带轮包角（角度）α=');disp(r);disp('V带根数 Z=');disp(z);disp('V带的基准长度(mm) Ld=');disp(ld);disp('初拉力(N) F0=');disp(f0);disp('作用在轴上的压力(N) FQ=');disp(fq);disp('带轮轮缘宽度(mm) B=');disp(bb);endendt1=input('是否计算下一种带轮直径？ 1、yes；2、no。