计算机仿真课程实验指导.docVIP

计算机仿真 课程实验指导 河北工业大学电气工程学院 生物医学工程专业 2011年秋季学期  目 录 实验项目一：MATLAB数值计算 3 实验项目二：MATLAB数据的可视化 7 实验项目三：动态电路的仿真 10 实验项目四：模拟电路的仿真 11 实验项目五：组合逻辑电路仿真 13 实验项目六：信源编码和量化 15  实验项目一：MATLAB数值计算 实验的目的和任务：学习使用MATLAB常用命令进行数组和矩阵等数值计算 实验内容： 要求在闭区间上产生具有10个等距采样点的一维数组。试用两种不同的指令实现。（参考提示：用冒号或linspace命令） 已知矩阵，运行指令B1=A.^(0.5), B2=A^(0.5), 可以观察到不同运算方法所得结果不同。 （1）请分别写出根据B1, B2恢复原矩阵A的程序。 （2）用指令检验所得的两个恢复矩阵是否相等。 先运行clear,format long,rand(state,1),A=rand(3,3)，然后根据A写出两个矩阵：一个对角阵B，其相应元素由A的对角元素构成；另一个矩阵C，其对角元素全为0，而其余元素与对应的A阵元素相同。 clear,format long,rand(state,1),A=rand(3,3) A = 0.95278214965662 0.59815852417219 0.83681960067634 0.70406216677500 0.84074319811307 0.51870305972492 0.95387747359223 0.44281884223513 0.02220977857260 命令为： 下面有一段程序，企图用来解决如下计算任务：有矩阵，当依次取10, 9, 8, 7, 6, 5, 4, 3, 2, 1时，计算矩阵“各列元素的和”，并把此求和结果存放为矩阵Sa的第k行。例如时，A阵为，此时它各列元素 的和是一个行数组，并把它保存为Sa的第3行。问题：该段程序的计算结果对吗？假如计算结果不正确，请指出错误发生的根源，并改正之。 原程序为： for k=10:-1:1 A=reshape(1:10*k,k,10); Sa(k,:)=sum(A); end Sa; 修改源程序： 由指令rand(state,0),A=rand(4,6)生成二维数组A，试求该数组中所有大于0.5的元素的位置，分别求出它们的“全下标”和“单下标”。（请参考find指令的用法，利用help find，注：命令(state,0) 为resets the generator to its initial state.?重新恢复初始状态，刚运行matlab时候的状态，更详细的rand 用法请参考help rand） 求解一元七次方程1.25x7+5.66x6+0.71x4+0.68x3+6.21x+5015＝0的根。（使用root指令） 命令为： 7.运行如下指令：（请参考例2.1-3） syms a b x X Y k=sym(‘3’); z=sym(‘c*sqrt(delta)+y*sin(theta)’) EXPR=a*z*X+(b*x^2+k)*Y 在表达式EXPR中，确定1个、2个和3个自由符号变量的命令分别为： 找出EXPR中的全部的自由符号变量的命令为： 8.请编写MATLAB指令，化简。（请参考例2.2-2） 9.请编写MATLAB指令求。（请参考例2.3-1） 10. 请编写MATLAB指令求。（请参考例2.3-9） 11. 产生7阶随机方阵A，其元素为[10,900]区间的随机整数，然后判断A的元素是否能被7整除. 12. 已知矩阵 1 2 0 X＝ 2 5 -1 4 10 -1 计算其特征值V和特征向量D，并验证XV=VD。（提示用eig函数） 13. 用求特征值的方法解方程，3x5-7x4+5x2+2x-18=0。（参考课件例3-3） 14、求方程组的解（参考课件中代数方程组求解） （1） x1+2x2=8 2x1+3x2=13 （2） x1+2x2=1 2x1+3x2=2 3x1+4x2=3 （3） x1+2x2+3x3=1 2x1+3x2+4x3=2 实验项目二：MATLAB数据的可视化 实验的目的和任务：利用MATLAB常用命令进行符号计算和绘图，练习二维图形和三维图形的绘制、符号计算 实验