matlab-计算机仿真作业.docx

计算机仿真大作业 石家庄铁道大学 数理系数0802 至尊 2011年4月16日 1?典型直流电阻电路的分析计算 用回路电流法求解图示电路中电流I及电压U 解：应用回路电流法，可列出以下方程： Im2-1.4Im3=0 -12.5lm1-2lm2-2.5lm3=14 21.5lm3-2.5lm1-15lm2=0 其中I=-Im3 U=17lm2-2lm1-15lm3 「0 1 -1.4, Im1 0 12.5 2 2.5 Im2 = 14 「2.5 -15 21.5 Jm3 0, 求得 Im1,lm2,lm3 再得 I,U MATLAB语言编程代码： clear; syms Im1; syms Im2; syms Im3; A=[0 1 -1.4;12.5 2 2.5;-2.5 -15 21.5];列出系数矩阵 A B=[lm1;lm2;lm3]; %列出系数矩阵B C=[0;14;0]; %列出系数矩阵C B=A\C; %解出 B I=-B(3); U=-2*B(1)+17*B(2)-15*B(3); I,u 程序运行结果：I二-1.8 U = 15.1 此方法适合多数题目，但列出的方程组必须满足可以有逆矩阵的 条件，即行列式结果不为零。正如程序段“ B=A\C; %解出B”要求 矩阵A必须有逆矩阵。 应用戴维南等效定理求解电路： 如图所示，已知 isi=2.5A,R=15@ ,「2=5q ,「3=4q ,Us2=10V,求该含源一端 口 1-1的戴维南等效电路。 图2应用戴.维南的电路 解：首先，建立数学模型。 利用节点分析法分别求解开路电压 Uoc和短路电流isc及等效电阻 R0,则有: （盒+盒）妆冃讪+瓷 （盒 +盒）妆冃讪+瓷 =藝 isc MATLAB语言编程代码： clear; r1=15;r2=5;r3=4;us2=10;is 仁 2.5; %给定电路参数 a1=1/r1+1/r2; %-端口开路时自导纳 uoc=(us2/r2+is1)/a1; %求解开路电压 a2=1/r1+1/r2+1/r3; %求解短路电流 isc=(us2/r2+is1)/a2/r3; %求解戴维南等效电阻 ro=uoc/isc; ro uoc 结果 ro = 7.7500 uoc =16.8750 此方法适用于戴维南等效电路。但对于不需要化成戴维南等效电路 的题目，还是采用一般方法较为简单，在计算个编程方面都会节省 时间和精力。 时间步长法： 求由曲线(x-1)A2+yA2=4, xA2/9+yA2=1,以及 xA2+(y+1)A2=4 所围成图形 的面积。 画图： t=0:0.01:2*pi; x=sin(t); y=cos(t); plot(2*x+1,2*y,3*x,y,2*x,2*y-1) axis([-3,3,-3,3]) 3 Matlab程序为 x=-3:0.01:3; y=-3:0.01:3; s=0; h=0.01; for i=1:800 for j=1:400 xx=-2+i*h; yy=-2+j*h; if (xx-1)A2+yyA2=4 if xxA29+yyA2=1 if xxA2+(yy+1F2v=4 s=s+hA2; end end end end end s 计算结果是： s = 4.7581 时间步长法使用条件限制比较小，一般在编程时，应注意取用适当 的时间步长，平衡精确性和程序运行步骤，可以更简单直接的求得 事件步长法： 对于一个单服务台的售票处，假设各顾客的到达时间间隔和服务时间均服从负指数分布， 到达率 人入=10人/分钟，服务率 卩=6人/分钟，等待队长 N=20人，仿真时间10分钟。 MATLAB 语言编程代码： clear clc %***************************************** % 初始化顾客源 %***************************************** %总仿真时间 Total_time = 10; %队列最大长度 N = 20; %到达率与服务率 lambda = 10; mu = 6; %平均到达时间与平均服务时间 arr_mean = 1/lambda; ser_mean = 1/mu; %可能到达的最大顾客数 round 四舍五入求整数 arr_num = round(Total_time*lambda*2); %顾客事件表初始化 events = []; %按负指数分布产生各顾客达到时间间隔 events(1,:) = exprnd(arr_mean,1,arr_num); %各顾客的到达时刻等于时间间隔的累积和 events(1,:) = cumsum(events(1,:)); %按负指数分布产生各顾客服