MATLAB课程作业练习题

内蒙古工业大学电力学院

《MATLAB在电力系统中的应用》

练习题汇总

编者：顾秀芳

单位：电力学院电力工程系

时间：2012年8月

实验一：《MATLAB操作环境及基本操作》实验指导书

1、实验目的：了解MATLAB的运行环境，界面功能；熟悉MATLAB编程的基

本方法；掌握MATLAB的基本画图功能；熟悉Simulink模块的基本操作及子系统的封装方法。

2、实验内容：操作MATLAB界面的主要菜单并了解其中的功能；练习矩阵的

运算；编写或读懂并运行MATLAB程序；练习并实践MATLAB的画图功能的实现方法。练习Simulink模块的基本操作；熟悉并练习子系统的封装方法；

3、实验要求：通过预习，了解MATLAB的界面的组成，各界面实现的主要功

能；了解MATLAB实现运算、逻辑判断的基本语法、函数的含义等；掌握MATLAB基本绘图功能，并实现简单的绘图；熟悉Simulink的及子系统封装的主要功能及基本操作。

4、实验步骤：

（1）熟悉MATLAB界面的基本功能；

（2）熟悉编程基础，并读通指定程序的含义，运行结果，并分析；

（3）熟悉MATLAB的基本绘图语句，读通指定并读通指定程序的含义，运行结果并分析；实现给定的绘图任务；

（4）熟悉Simulink模块的操作，并完成给定模型的搭建；

（5）熟悉Simulink子系统模型搭建的步骤，并完成给定子系统的搭建及参数设置。

6 、实验题目：

1）MATLAB环境的了解

包括：窗口菜单的主要功能；各窗口实现的功能，包括：命令窗口，workspace/history,command history三个窗口的应用。

2）MATLAB编程基础

内容包括：矩阵的形成及运算（矩阵相乘、点乘）；逻辑运算；函数的应用。 练习内容一：矩阵的形成及操作；

1：形成矩阵1

0.1 2.136016 2.8 3.42.10.10.30.4 1.31.4 5.10 1.10A ⎡⎤⎢⎥-⎢

⎥=⎢⎥-⎢⎥

-⎣⎦ 2:从矩阵中提取子矩阵： 3：矩阵A 中元素A41的值； 4：修改矩阵A 的元素，使A41=3；

5：矩阵A 中最后2行和最后3行交汇的子矩阵的值； 练习内容二：矩阵的运算 1：矩阵定义为：

22111,,12022A B C --⎡⎤⎡⎤⎡⎤

===⎢⎥⎢⎥⎢⎥--⎣⎦⎣⎦⎣⎦

求出下列矩阵运算的结果，并解释部分操作无法实现的原因： （1） R=A+B; （2） R=A*D; （3） R=A.*D; （4） R=A*C; （5） R=A.*C; （6） R=A\B; （7） R=A.\B; （8） R=A.^B.

1.网络如图，图中标注为支路导纳yij.试通过程序完成Y 导纳矩阵的形成

ii ij i j

Y y ∈=∑，j 含对地支路

ii ij Y y =-

20

B

练习内容三：编程语句的理解

1:运行书P20中例题，比较标量循环和矩阵运算的差异；

比较程序：

1.一个窗口画一个图

t=0:0.1:2*pi;

y1=sin(t);

plot(t,y1,'r:')

axis([0,2*pi,-2,2]);

xlabel('huduzhi');ylabel('hanshuzhi');

title('正弦函数和余弦函数');

grid on

gtext({'y1=sin(t)'});

2.一个窗口绘两个图

(1)一条plot语句画

t=0:0.1:2*pi;

y1=sin(t); y2=cos(t);

plot(t,y1,'r:', t,y2,'b-.')

axis([0,2*pi,-2,2]);

xlabel('t');ylabel('y1,y2');

title('正弦函数和余弦函数');

grid on

gtext({'y1=sin(t)'; 'y2=cos(t)'});

(1)分多条plot语句画

t=0:0.1:2*pi;

y1=sin(t); y2=cos(t); y3=tan(t)

plot(t,y1,'r-')

hold

plot(t,y2,'b--')

plot(t,y3,'k-.')

axis([0,2*pi,-2,2]);

xlabel('t');ylabel('y1,y2,y3');

grid on

gtext({'y1=sin(t)'; 'y2=sin(t)';'y3=sin(t)'}); legend('y1','y2','y3')

3、一个窗口分割成三个小窗口

t=0:0.1:2*pi;

y1=sin(t);

subplot(3,1,1)

plot(t,y1,'r:')

axis([0,2*pi,-2,2]);

xlabel('t');ylabel('y1');

title('正弦函数');

grid on

gtext('y1=sin(t)');

subplot(3,1,2)

plot(t,cos(t),'r:')

axis([0,2*pi,-2,2]);

xlabel('t');ylabel('cost');

title('余弦函数');

grid on

gtext('y2=sin(t)');

subplot(3,1,3)

plot(t,tan(t),'r:')

axis([0,2*pi,-2,2]);

xlabel('t');ylabel('tant');

title('正切函数');

grid on

gtext('y3=sin(t)');

作业：

t=0:0.1:2*pi;%

y1=sin(t); y2=cos(t); y3=tan(t);y4=cot(t) subplot(2,2,1)

plot(t,y1,'r-')

axis([0,2*pi,-2,2]);

xlabel('t');ylabel('y1');

subplot(2,2,2)

plot(t,y2,'b--*')

axis([0,2*pi,-2,2]);

xlabel('t');ylabel('y2');

subplot(2,2,3)

plot(t,y3,'k:')

axis([0,2*pi,-2,2]);

xlabel('t');ylabel('y3');

subplot(2,2,4)