姓名：孟云飞学号Y080430 任意信道的信道容量C及最佳输入分布P_X源程序
function [P_X,C,N]=contmax(P_YX,e)
% 计算任意信道的信道容量C及最佳输入分布P_X
% P_X -------输入概率分布
% C----------信道容量N----------迭代次数
% P_YX-----DMC信道的转移矩阵e-----------精度
if length(find(P_YX<0)~=0)
error('Not a probable vector.Negtive component');
end
if abs(sum(P_YX')-1)>10e-10
error('Not a probable vector,Component do not add up to "1" ');
end
% 变量初始化
C1=1;C=0;N=0;
r=size(P_YX);P_X=ones(1,r(1))/r(1);
% 调整P_YX的零元素值
Pyx=(P_YX==0).*eps;P_YX=P_YX+Pyx;
% 迭代求解
while (abs(C1-C))>e
P_Y=P_X*P_YX;
I1=sum((P_YX.*log2(P_YX))');
I2=log2(P_Y)*(P_YX'); BETA=exp(I1-I2);
B=P_X*(BETA');
C1=log(B);C=log(max(BETA));
P_X=P_X.*BETA/B;
N=N+1;
end
P_YX =
0.5000 0.3000 0.2000
0.4000 0.3000 0.3000
0.1000 0.6000 0.4000
>> [P_X,C,N]=contmax(P_YX,0.001)
P_X =
0.3599 0.0023 0.6378
C =
0.1596
N =
40
>> [P_X,C,N]=contmax(P_YX,0.0001)
P_X =
0.3620 0.0002 0.6378
C = 0.1591
N = 58
由上面结果可以看出，不同的停止迭代容差限对信道容量没有太大影响
P_YX =
0.3333 0.3333 0.3333
0.3333 0.3333 0.3333
0.3333 0.3333 0.3333
>> [P_X,C,N]=contmax(P_YX,0.00001)
P_X =
0.3333 0.3333 0.3333
C =
N = 1
P_YX =
0.5000 0.3000 0.2000
0.4000 0.3000 0.3000
0.3000 0.3000 0.4000
>> [P_X,C,N]=contmax(P_YX,0.00001)
P_X =
0.5072 0.0002 0.4926
C =
0.0427
N = 181
P_YX =
0.5000 0.3000 0.2000
0.4000 0.3000 0.3000
0.1000 0.6000 0.4000
>> [P_X,C,N]=contmax(P_YX,0.00001)
P_X =
0.3622 0.0000 0.6378
C =
0.1590
N = 76
P_YX =
0.8000 0.1000 0.1000
0.2000 0.7000 0.1000
0.1000 0.3000 0.6000
>> [P_X,C,N]=contmax(P_YX,0.00001)
P_X =
0.4114 0.2412 0.3474
C =
0.4551
N =
26
由上面的结果可以看出信道容量和转移概率矩阵有关，不同的信道参数信道容量不同，当转移概率矩阵的每行的概率越趋于均匀分布，信道容量就越小。
潮流例题：根据给定的参数或工程具体要求（如图），收集和查阅资料；学习相关软件（软件自选：本设计选择Matlab进行设计）。 2.在给定的电力网络上画出等值电路图。 3.运用计算机进行潮流计算。 4.编写设计说明书。 一、设计原理 1．牛顿-拉夫逊原理 牛顿迭代法是取x0 之后，在这个基础上，找到比x0 更接近的方程的跟，一步一步迭代，从而找到更接近方程根的近似跟。牛顿迭代法是求方程根的重要方法之一，其最大优点是在方程f(x) = 0 的单根附近具有平方收敛，而且该法还可以用来求方程的重根、复根。电力系统潮流计算，一般来说，各个母线所供负荷的功率是已知的，各个节点电压是未知的（平衡节点外）可以根据网络结构形成节点导纳矩阵，然后由节点导纳矩阵列写功率方程，由于功率方程里功率是已知的，电压的幅值和相角是未知的，这样潮流计算的问题就转化为求解非线性方程组的问题了。为了便于用迭代法解方程组，需要将上述功率方程改写成功率平衡方程，并对功率平衡方程求偏导，得出对应的雅可比矩阵，给未知节点赋电压初值，一般为额定电压，将初值带入功率平衡方程，得到功率不平衡量，这样由功率不平衡量、雅可比矩阵、节点电压不平衡量（未知的）构成了误差方程，解误差方程，得到节点电压不平衡量，节点电压加上节点电压不平衡量构成新的节点电压初值，将新的初值带入原来的功率平衡方程，并重新形成雅可比矩阵，然后计算新
的电压不平衡量，这样不断迭代，不断修正，一般迭代三到五次就能收敛。 牛顿—拉夫逊迭代法的一般步骤： （1）形成各节点导纳矩阵Y。 （2）设个节点电压的初始值U和相角初始值e 还有迭代次数初值为0。 （3）计算各个节点的功率不平衡量。 （4）根据收敛条件判断是否满足，若不满足则向下进行。 （5）计算雅可比矩阵中的各元素。 （6）修正方程式个节点电压 （7）利用新值自第（3）步开始进入下一次迭代，直至达到精度退出循环。 （8）计算平衡节点输出功率和各线路功率 2．网络节点的优化 １）静态地按最少出线支路数编号 这种方法由称为静态优化法。在编号以前。首先统计电力网络个节点的出线支路数，然后，按出线支路数有少到多的节点顺序编号。当由n 个节点的出线支路相同时，则可以按任意次序对这n 个节点进行编号。这种编号方法的根据是导纳矩阵中，出线支路数最少的节点所对应的行中非零元素也２）动态地按增加出线支路数最少编号在上述的方法中，各节点的出线支路数是按原始网络统计出来的，在编号过程中认为固定不变的，事实上，在节点消去过程中，每消去一个节点以后，与该节点相连的各节点的出线支路数将发生变化(增加，减少或保持不变)。因此，如果每消去一个节点后，立即修正尚未编号节点的出线支路数，然后选其中支路数最少的一个节点进行编号，就可以预期得到更好的效果，动态按最少出线支路数编号方法的特点就是按出线最少原则编号时考虑了消去过程中各节点出线支路数目的变动情况。 3．MATLAB编程应用 Matlab 是“Matrix Laboratory”的缩写，主要包括：一般数值分析，矩阵运算、数字信号处理、建模、系统控制、优化和图形显示等应用程序。由于使用Matlab 编程运算与人进行科学计算的思路和表达方式完全一致，所以不像学习高级语言那样难于掌握，而且编程效率和计算效率极高，还可在计算机上直接输出结果和精美的图形拷贝，所以它的确为一高效的科研助手。 二、设计内容 1.设计流程图
华中科技大学 信息工程学院课程设计报告书题目: 电力系统潮流计算 专业：电气工程及其自动化 班级： 学号： 学生姓名： 指导教师： 2015年 11 月 10 日
2015年11月12日
信息工程学院课程设计成绩评定表
摘要 电力系统稳态分析包括潮流计算和静态安全分析。本文主要运用的事潮流计算，潮流计算是电力网络设计与运行中最基本的运算，对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。本位就是运用潮流计算具体分析，并有MATLAB仿真。 关键词：电力系统潮流计算 MATLAB仿真
Abstract Electric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation. Key words: Power system; Flow calculation; MATLAB simulation
基于MATLAB的计算器制作从入门到高级 项目一、建立GUI框架 1，双击MATLAB图标（本实例使用的是MATLAB R2014），打开软件 2，打开GUI操作界面。有两种方式： a、在命令行窗口输入guide,然后点击回车。 b、先点击新建图标，在选择‘图形用户界面’。 得到如下窗口： 左上方有两个按钮：‘新建GUI’，‘打开现有GUI’。 对于‘新建GUI’：下方的选择框内有你可以选择建立的GUI类型，对于初学者一般选择第一个。下方是用来设置 保存文件位置的。
对于‘打开现有GUI’：初学者可能会在网上下载一些GUI源程序，使用‘打开现有GUI’，再浏览到文件所在位置便可。 3、构建GUI框架： 在选择新建GUI并且确定后，会弹出如下窗口： 在窗口的左边会有很多选项，可以一个个拖动编辑框内； 是不是发现这些选项发挥的作用刚好就是我们平常打开一个网页或者一个软件进行操作时的选择按钮。比较常见的有：按钮、单选按钮、拖动条、弹出式菜单等。看到这里细心地朋友就会发现，MATLAB GUI似乎可以开发出大型
软件哦，只要你有足够的想象力，当你从一个小白进阶到 一个高手之后，你可以做到的。 本次计算器制作只需要用到静态文本（或者是动态文本）用来做显示器，按钮（用来操作）。也可以添加一个面板 用来修饰。 在拖出的选项上双击，会得出该选项的属性栏，比如说 双击一个按钮选项，你会得到： 这个属性框比较复杂，没有必要一个个说清楚（有兴趣的 可以自己研究，其实理解起来也挺简单的） 本次计算器制作，我们只需要知道以下几个：
◆a、BackgroundColor取值为颜色的预定义字符或RGB 数值；缺省（就是默认的意思）值为浅灰色； ◆b、ForegroundColor取值为颜色的预定义字符或RGB 数值，该属性定义控件对象标题字符的颜色；缺省值为黑色； ◆c、String取值为字符串矩阵或块数组，定义控件对 象标题或选项内容； ◆d、FontName取值为控件标题等字体的字库名； ◆e、FontSize取值为数值，控制字体大小； ◆f、Tag取值为字符串，定义了控件的标识值，在任何 程序中都可以通过这个标识值控制该控件对象； ◆g、Style取值可以是pushbutton(缺省值), radiobutton, checkbox, edit, text, slider, frame, popupmenu 或listbox； 项目二、简单加法器的制作 1、框架准备
附录1 使用牛顿拉夫逊法进行潮流计算的Matlab程序代码 % 牛拉法计算潮流程序 %----------------------------------------------------------------------- % B1矩阵：1、支路首端号；2、末端号；3、支路阻抗；4、支路对地电纳 % 5、支路的变比；6、支路首端处于K侧为1，1侧为0 % B2矩阵：1、该节点发电机功率；2、该节点负荷功率；3、节点电压初始值 % 4、PV节点电压V的给定值；5、节点所接的无功补偿设备的容量 % 6、节点分类标号：1为平衡节点（应为1号节点）；2为PQ节点；3为PV节点； %------------------------------------------------------------------------ clear all; format long; n=input('请输入节点数:nodes='); nl=input('请输入支路数:lines='); isb=input('请输入平衡母线节点号:balance='); pr=input('请输入误差精度:precision='); B1=input('请输入由各支路参数形成的矩阵:B1='); B2=input('请输入各节点参数形成的矩阵:B2='); Y=zeros(n);e=zeros(1,n);f=zeros(1,n);V=zeros(1,n);sida=zeros(1,n);S1=zeros(nl); %------------------------------------------------------------------ for i=1:nl %支路数 if B1(i,6)==0 %左节点处于1侧 p=B1(i,1);q=B1(i,2); else %左节点处于K侧 p=B1(i,2);q=B1(i,1); end Y(p,q)=Y(p,q)-1./(B1(i,3)*B1(i,5)); %非对角元 Y(q,p)=Y(p,q); %非对角元 Y(q,q)=Y(q,q)+1./(B1(i,3)*B1(i,5)^2)+B1(i,4); %对角元K侧 Y(p,p)=Y(p,p)+1./B1(i,3)+B1(i,4); %对角元1侧 end %求导纳矩阵 disp('导纳矩阵 Y='); disp(Y) %------------------------------------------------------------------- G=real(Y);B=imag(Y); %分解出导纳阵的实部和虚部 for i=1:n %给定各节点初始电压的实部和虚部
1.需求分析 本次的实验要求是设计一个计算器，主要功能如下： （1）实现基本数学运算（加减乘除等），而且要能进行混合运算 （2）实现部分函数功能，如求平方根、求倒数等 （3）能实现小数运算 界面与标准计算器界面类似 根据要求以及以前的学习情况，决定使用matlab进行编程。Matlab强大的计算功能以及便捷的GUI设计，可以较为简便的实现所要求的功能。按照要求，数据输入和输出支持小数点，支持四则混合运算，决定使用如下几个数据进行分析：（1+3）*5 Sqrt（4） 1/2 Sin4 用以检验是否可以进行加减乘除四则运算、平方根、倒数、正弦的运算。 2.程序设计 M atlab的程序设计较为简便，用GUI设计出一个计算器的模型，然后系统会自动生成一个框架，在框架中，写入每一个按键对应的程序就可以实现功能。 3.调式分析 编程的过程中遇到的问题不是很多，基本就是找要实现各个功能的子程序，通过上网和去图书馆，加上自己的编写，终于实现了实验要求的功能。但是有一点很重要，matlab不支持中文，所以从路径到文件名必须是全英文的，不然就无法识别。此外，给每个按键命名也是很重要的，不然在生成的程序框架里面，就无法识别各个按键的作用，编写程序的时候也就无法做到一一对应。 4.使用说明 程序的使用比较简单，由于是可视化界面，直接打开matlab，然后建立一个GUI 工程，再打开生成的fig文件，就是一个计算器的界面，直接按照市面上卖的计算器的
方法，按键使用即可。 5.测试结果 计算结果为20 4sqrt=2 Sin4结果为 1/2=0.5 经过计算，这些结果均与实际结果相吻合，计算器的功能实现的较为完好。 6.心得体会 本次试验由于不限制语言，于是计算功能强大，操作简便的matlab变成了首选，matlab的GUI设计，操作是较为简单的，首先建立一个GUI工程，然后用可视化界面，
一．简答题(每题5分，共40分) 1．数组运算和矩阵运算的运算符有什么区别？ 答：在加、减运算时数组运算与矩阵运算的运算符相同，乘、除和乘方运算时，在矩阵运算的运算符前加一个点即为数组运算，如a*b为矩阵乘，a.*b为数组乘。 2. 命令文件与函数文件的主要区别是什么？ 答：（1）命令文件是一系列命令的组合，函数文件的第一行必须用function说明； （2）命令文件没有输入参数，也不用返回参数，函数文件可以接受输入参数，也可以返回参数； （3）命令文件处理的变量为工作空间变量，函数文件处理的变量为函数内部的局部变量，也可以处理全局变量。 3. 如何定义全局变量？ 答：用关键字global可以把一个变量定义为全局变量，在M文件中定义全局变量时，如果在当前工作空间已经存在了相同的变量，系统将会给出警告，说明由于将该变量定义为全局变量，可能会使变量的值发生改变，为避免发生这种情况，应该在使用变量前先将其定义为全局变量。 4. 什么是Simulink ? 答：Simulink是一个用来对动态系统进行建模、仿真和分析的软件包，它支持连续的、离散的或二者混合的线性和非线性系统，也支持具有多种采样速率的多速率系统。 5. 如何在Simulink中进行指定仿真时间的操作? 选择Simulation菜单打开参数设置对话框，在Solver页面中第一项Simulation time 中通过输入开始时间和结束时间来确定仿真时间。 6. MATLAB的变量命名有什么规则？能用”NaN”来作为变量名吗，为什么？ 答：变量可以包括数字、字母、下划线，必须以字母开头 可以，因为符合变量名命名规则 7. 简述MATLAB命令窗的主要作用？ 答：命令窗口是MATLAB的主要交互窗口，用于输入命令并显示除图形以外的所有执行结果。是MATLAB提供给用户使用的管理功能的人机界面，其管理功能包括：管理工作空间中的变量、数据的输入输出的方式和方法，开发、调试、管理M 文件和M函数的各种工具。 8. 简述MATLAB桌面集成环境包括的7个窗口及其主要功能？ 答：桌面是MATLAB的主要工作界面。命令窗口是MATLAB的主要交互窗口，用于输入命令并显示除图形以外的所有执行结果。工作空间管理窗口用于显示当前计算机内存中MATLAB变量的名称、数学结构、该变量的字节数及其类型，可对变量进行观察、编辑、保存和删除。命令历史窗口显示用户在命令窗口中所输入的每条命令的历史记录，并标明使用时间，这样可以方便用户查询。当前目录窗口显示当前用户工作所在的路径。Start菜单用于选择菜单中的命令执行MATLAB产品的各种工具，并且可以查阅MATLAB所含的各种资源。编译窗口为用户提供了一个图形界面进行M文件的编写和调试。 二．编程题(每题10分，共60分) 1.编写M文件，使用for循环计算f=5^x，当f（x）>500就终止程序，在命令窗口中调用函数文件，求终止程序时f的值和n运行的值分别为多少? 解： function[y]=f(x) % 2分 for x=1:00 % 4分 f=5^x if f>500 f=z n=x break % 4分 end end 2. 编写程序分别对符号表达式23 3546 f x y x y =+-+的x和y求一阶微分。 解：smys x y % 2分 f=3*x^2+5*y^3-4*x+6*y; % 4分 dfdx=diff(f)% 对x求一阶微分% 2分 dfdy=diff(f,y)% 对y求一阶微分% 2分 3. 写出用黑色画出一条正弦虚线和一条余弦点划线，时间为0到2π并标记横坐标和纵坐标，并限定x轴和y轴的显示范围，并添加图例的步骤。 解： x=0:pi/10:2*pi; y1=sin(x); y2=cos(x); %4分 plot(x,y1,’:b’,x,y2,’-.b’) % 2分 xlabel(‘x’) ylabel(‘y’) %2分 axis([0 2*pi -1 1]) %2分 第1页（共2页）
课程设计任务书 学生姓名：专业班级： 指导教师：工作单位： 题目: 基于MATLAB的图像处理的基本运算 初始条件： 要求完成的主要任务:（包括课程设计工作量及其技术要求，以及说明书撰写等具体要求） （1）能够对图像亮度和对比度变化调整，并比较结果 （2）编写程序通过最近邻插值和双线性插值等算法将用户所选取的图像区域进行放大和缩小整数倍的和旋转操作，并保存，比较几 种插值的效果 （3）图像直方图统计和直方图均衡，要求显示直方图统计，比较直方图均衡后的效果。 （4）对图像加入各种噪声，比较效果。 时间安排： 指导教师签名：年月日 系主任（或责任教师）签名：年月日 目录 摘要.......................................................................................................................... 错误!未定义书签。 1 MATLAB简介 ........................................................................................................ 错误!未定义书签。2图像选择及变换................................................................................................... 错误!未定义书签。 2.1 原始图像选择读取....................................................................................... 错误!未定义书签。 2.1.1 原理图的读入与基本变换 .................................................................... 错误!未定义书签。
信号与系统的MATLAB 仿真 一、信号生成与运算的实现 1.1 实现)3(sin )()(π±== =t t t t S t f a )(sin )sin()sin(sin )()(t c t t t t t t t S t f a '=' '== ==πππ π ππ m11.m t=-3*pi:0.01*pi:3*pi; % 定义时间范围向量t f=sinc(t/pi); % 计算Sa(t)函数 plot(t,f); % 绘制Sa(t)的波形 运行结果： 1.2 实现)10() sin()(sin )(±== =t t t t c t f ππ m12.m t=-10:0.01:10; % 定义时间范围向量t f=sinc(t); % 计算sinc(t)函数 plot(t,f); % 绘制sinc(t)的波形 运行结果： 1.3 信号相加：t t t f ππ20cos 18cos )(+= m13.m syms t; % 定义符号变量t f=cos(18*pi*t)+cos(20*pi*t); % 计算符号函数f(t)=cos(18*pi*t)+cos(20*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果：
1.4 信号的调制：t t t f ππ50cos )4sin 22()(+= m14.m syms t; % 定义符号变量t f=(2+2*sin(4*pi*t))*cos(50*pi*t) % 计算符号函数f(t)=(2+2*sin(4*pi*t))*cos(50*pi*t) ezplot(f,[0 pi]); % 绘制f(t)的波形 运行结果： 1.5 信号相乘：)20cos()(sin )(t t c t f π?= m15.m t=-5:0.01:5; % 定义时间范围向量 f=sinc(t).*cos(20*pi*t); % 计算函数f(t)=sinc(t)*cos(20*pi*t) plot(t,f); % 绘制f(t)的波形 title('sinc(t)*cos(20*pi*t)'); % 加注波形标题 运行结果：
题目：潮流计算与matlab 教学单位电气信息学院姓名 学号 年级 专业电气工程及其自动化指导教师 职称副教授
摘要 电力系统稳态分析包括潮流计算和静态安全分析。本文主要运用的事潮流计算，潮流计算是电力网络设计与运行中最基本的运算，对电力网络的各种设计方案及各种运行方式进行潮流计算，可以得到各种电网各节点的电压，并求得网络的潮流及网络中的各元件的电力损耗，进而求得电能损耗。本位就是运用潮流计算具体分析，并有MATLAB仿真。 关键词：电力系统潮流计算 MATLAB Abstract Electric power system steady flow calculation and analysis of the static safety analysis. This paper, by means of the calculation, flow calculation is the trend of the power network design and operation of the most basic operations of electric power network, various design scheme and the operation ways to tide computation, can get all kinds of each node of the power grid voltage and seek the trend of the network and the network of the components of the power loss, and getting electric power. The standard is to use the power flow calculation and analysis, the specific have MATLAB simulation. Key words: Power system; Flow calculation; MATLAB simulation
目录 计算器的效果图 ..................................... 错误!未定义书签。 一、GUI设计界面： (3) 1.打开GUI (3) 2.添加按钮 (3) 3.根据按钮的作用及视觉效果做一定的修改： (4) 4.保存、添加功能函数 (4) （1）数字键编写 (4) （2）符号键的编写 (4) （3）运算符“=”的编写 (5) （4）按键“←back”的编写 (5) （5）按键“清空”的编写 (5) （6）按键“退出”的编写 (5) （7）按键“二进制数转十进制数”的编写 (5) （8）按键“十进制数转二进制数”的编写 (5) 二、计算器的使用 (6) 除法运算（÷） (6) 平方运算（^2） (6) 函数cos （∏/3）的计算 (7) 函数arctan （∏/3）的计算 (7) 以2为底的对数的计算（log 2） (8) 十进制数转二进制数的计算（调用dec2bin函数） (8) 二进制数转十进制数的计算（调用bin2dec函数） (8) 三、附各按键的程序源代码 (9) 四、问题和解决方法 (15) 五、心得体会 (15) 参考文献 (15)
计算器的效果图:
一、GUI设计界面： 1.打开GUI 输入Guide 回车或者在工具栏上点击图标打开Guide 窗口： 2.添加按钮
3.根据按钮的作用及视觉效果做一定的修改： 双击按钮（Puch Button）进入按键属性修改显示字符串大小、字体和颜色，然后对按钮的位置进行排布，尽量使按钮集中在静态文本框下面。 4.保存、添加功能函数 把做好的按钮及静态文本框保存后自动弹出Editor的M文本，对然后对相应的pushbutton添加功能函数。以下是相应按钮的功能函数。 （1）数字键编写 在function pushbutton1_Callback(hObject, eventdata, handles）下输入：textString = get(handles.text1,'String'); textString =strcat(textString,'0'); set(handles.text1,'String',textString) 这是使用句柄handles指向对象text1，并以字符串形式来存储数据文本框text1的容，并存储数个“0”， 然后由set(handles.text1,'String','textString'在text1中输出。 同理，分别在function pushbutton2~10_Callback(hObject, eventdata, handles）下给1~9数字按键下编写此类程序。 （2）符号键的编写 function pushbutton12_Callback(hObject, eventdata, handles) textString = get(handles.text1,'String'); textString =strcat(textString,'÷'); set(handles.text1,'String',textString) strcat的作用是将两个字符串连接起来，就是在已输入的存储数据textString 后添加“÷”进行运算。 然后执行set(handles.text1,'String',textString)。符号键‘-’、‘*’、‘/’与‘÷’的运算函数类似。“平方运算”,主要是由“^2”功能实现。
摘要 MATLAB是美国mathworks公司开发的大型数学计算软件，它具有强大的矩阵处理和绘图功能，已经广泛地应用于科学研究和工程技术的各个领域，MATLAB 以矩阵和向量为基本数据单元，提供了丰富的矩阵操作和矩阵运算功能，并在这些基本运算基础上提供了可供各种科学研究和工程技术门类使用的工具箱。极大地方便了科学计算和工程问题的求解，使得科技人员从复杂的变成工作中解放出来，专注于数学模型的建立。本文着重介绍MATLAB的基本操作和一些基础应用，并通过例子来阐述说明，熟悉MATLAB函数调用，了解其在电工电子技术中的应用。掌握Matlab基本语法结构及调试方法，熟悉Matlab函数调用，熟练其用于电工电子的一些计算，并学会用Matlab/Simulink进行简单电路的仿真。 一、MALTISM简介 MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。 在MATLAB函数库中，除了基本初等函数外，还有初等矩阵和矩阵变换、线性代数方程组合矩阵特征的求解、多项式运算和求根、函数的插值和数据的多项式拟合、数值积分和常数微分方程数值解、函数求极值、单变量非线性方程求解根、数据分析和傅里叶变换，以及某些特殊的矩阵函数和数学函数，这些函数都可以直接调用。用户可以根据自己的需要任意扩充函数库。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JAVA的支持。可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用. 二、MATLAB的应用。 MATLAB产品族可以用来进行以下各种工作： ●数值分析 ●数值和符号计算
基于MATLAB的电力系统潮流计算 %简单潮流计算的小程序，相关的原始数据数据数据输入格式如下： %B1是支路参数矩阵，第一列和第二列是节点编号。节点编号由小到大编写%对于含有变压器的支路，第一列为低压侧节点编号，第二列为高压侧节点%编号，将变压器的串联阻抗置于低压侧处理。 %第三列为支路的串列阻抗参数。 %第四列为支路的对地导纳参数。 %第五烈为含变压器支路的变压器的变比 %第六列为变压器是否是否含有变压器的参数，其中“1”为含有变压器，%“0”为不含有变压器。 %B2为节点参数矩阵，其中第一列为节点注入发电功率参数；第二列为节点%负荷功率参数；第三列为节点电压参数；第六列为节点类型参数，其中 %“1”为平衡节点，“2”为PQ节点，“3”为PV节点参数。 %X为节点号和对地参数矩阵。其中第一列为节点编号，第二列为节点对地%参数。 n=input('请输入节点数:n='); n1=input('请输入支路数:n1='); isb=input('请输入平衡节点号:isb='); pr=input('请输入误差精度:pr='); B1=input('请输入支路参数:B1='); B2=input('请输入节点参数:B2='); X=input('节点号和对地参数:X='); Y=zeros(n); Times=1; %置迭代次数为初始值 %创建节点导纳矩阵 for i=1:n1 if B1(i,6)==0 %不含变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/B1(i,3); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3)+0.5*B1(i,4); Y(q,q)=Y(q,q)+1/B1(i,3)+0.5*B1(i,4); else %含有变压器的支路 p=B1(i,1); q=B1(i,2); Y(p,q)=Y(p,q)-1/(B1(i,3)*B1(i,5)); Y(q,p)=Y(p,q); Y(p,p)=Y(p,p)+1/B1(i,3);
%*************************电力系统直角坐标系下的牛顿拉夫逊法潮流计算********** clear clc load E:\data\IEEE014_Node.txt Node=IEEE014_Node; weishu=size(Node); nnum=weishu(1,1); %节点总数 load E:\data\IEEE014_Branch.txt branch=IEEE014_Branch; bwei=size(branch); bnum=bwei(1,1); %支路总数 Y=(zeros(nnum)); Sj=100; %********************************节点导纳矩阵******************************* for m=1:bnum; s=branch(m,1); %首节点 e=branch(m,2); %末节点 R=branch(m,3); %支路电阻 X=branch(m,4); %支路电抗 B=branch(m,5); %支路对地电纳 k=branch(m,6); if k==0 %无变压器支路情形 Y(s,e)=-1/(R+j*X); %互导纳 Y(e,s)=Y(s,e); end if k~=0 %有变压器支路情形 Y(s,e)=-(1/((R+j*X)*k)); Y(e,s)=Y(s,e); Y(s,s)=-(1-k)/((R+j*X)*k^2); Y(e,e)=-(k-1)/((R+j*X)*k); %对地导纳 end Y(s,s)=Y(s,s)-j*B/2; Y(e,e)=Y(e,e)-j*B/2; %自导纳的计算情形 end for t=1:nnum; Y(t,t)=-sum(Y(t,:))+Node(t,12)+j*Node(t,13); %求支路自导纳 end G=real(Y); %电导 B=imag(Y); %电纳 %******************节点分类************************************* * pq=0; pv=0; blancenode=0; pqnode=zeros(1,nnum); pvnode=zeros(1,nnum); for m=1:nnum; if Node(m,2)==3 blancenode=m; %平衡节点编号 else if Node(m,2)==0 pq=pq+1; pqnode(1,pq)=m; %PQ 节点编号 else if Node(m,2)==2 pv=pv+1; pvnode(1,pv)=m; %PV 节点编号 end end end end %*****************************设置电压初值********************************** Uoriginal=zeros(1,nnum); %对各节点电压矩阵初始化 for n=1:nnum Uoriginal(1,n)=Node(n,9); %对各点电压赋初值 if Node(n,9)==0;
工程设计报告 设计题目：基于Matlab的简易计算器 学院： 专业： 班级： 学号： 姓名： 电子邮件： 日期：2015年12 月 成绩： 指导教师：
西安电子科技大学 电子工程学院 工 程设计 任务书 学生姓名指导教师职称 学生学号专业 题目基于Matlab 的简易计算器 任务与要求 任务如下： 利用MATLAB GUI 设计实现一个图形用户界面的计算器程序，实现： A.实现十进制数的加、减、乘、除、简单计算。 B. 科学计算函数，包括正弦、余弦、正切、余切、开方、指数等函数运行。 C. 有清除键，能清除操作。 要求如下： A ．熟练掌握MatlabGUI 界面的设计与应用 B ．最终计算器能够实现预期的相关功能 开始日期2015年 11月日完成日期2016年1月日 课程设计所在单位 本表格由电子工程学院网络信息中心编辑录入 https://www.doczj.com/doc/f110887948.html,. …………………………装…………………… … … … … 订 … … … … … … … … … … … …线 … …… …… …… …… …… … …… …… …… …… …… … …
摘要 基于Matlab GUI计算器设计时利用GUI的创建图像用户界面进行计算器设计。设计计算器时，主要是考虑到计算器的易用性、功能的常用程度进行计算器界面与功能的设计。通过调整控件和文本的布局及颜色，使界面简单大方、布局合理，达到界面友好的效果。 计算器设计时主要利用到get和set两个函数进行各个控件属性值的传递和设置。计算器实现的功能有：数字0~9和小数点的输入显示，平方开方和对数的输入显示。进行四则运算、正弦函数、余弦函数、正切函数以及反正弦函数、反余弦函数、反正切函数的计算等等。最后运行调试，实现基于MatlabGUI的计算器的设计。 关键词：MatlabGUI计算器 Abstracts Based on Matlab GUI calculator design using the user interface to create images of GUI calculator design.Design calculator, mainly considering the ease of use, function calculators calculator interface and function of the common level of design.By adjusting the control and the layout of the text and color, make the interface simple and easy, rational layout, to achieve the effect of friendly interface. Calculator design used to get and set two main function for each attribute value transfer and control Settings.Calculator the functions are: 0 ~ 9, according to input and decimal square root and logarithm of input.Arithmetic, sine function and cosine function, tangent function and the arcsine function,arccosine function, the calculation of the arctangent function and so on.Finally running debugging, implementation design based on Matlab GUI calculator. Keywords: Matlab GUI calculator
实验一 MATLAB 及仿真实验（控制系统的时域分析） 一、实验目的 学习利用MATLAB 进行控制系统时域分析，包括典型响应、判断系统稳定性和分析系统的动态特性； 二、预习要点 1、 系统的典型响应有哪些？ 2、 如何判断系统稳定性？ 3、 系统的动态性能指标有哪些？ 三、实验方法 （一） 四种典型响应 1、 阶跃响应： 阶跃响应常用格式： 1、)(sys step ；其中sys 可以为连续系统，也可为离散系统。 2、),(Tn sys step ；表示时间围0---Tn 。 3、),(T sys step ；表示时间围向量T 指定。 4、),(T sys step Y =；可详细了解某段时间的输入、输出情况。 2、 脉冲响应： 脉冲函数在数学上的精确定义：0 ,0)(1)(0 ?==?∞ t x f dx x f 其拉氏变换为：)()()()(1 )(s G s f s G s Y s f === 所以脉冲响应即为传函的反拉氏变换。 脉冲响应函数常用格式： ① )(sys impulse ； ② ); ,();,(T sys impulse Tn sys impulse ③ ),(T sys impulse Y = （二） 分析系统稳定性 有以下三种方法： 1、 利用pzmap 绘制连续系统的零极点图； 2、 利用tf2zp 求出系统零极点； 3、 利用roots 求分母多项式的根来确定系统的极点 （三） 系统的动态特性分析 Matlab 提供了求取连续系统的单位阶跃响应函数step 、单位脉冲响应函数impulse 、零输入响应函数initial 以及任意输入下的仿真函数lsim.
毕业设计（论文）MATLAB下的潮流计算实现-稀疏技术
毕业设计（论文）原创性声明和使用授权说明 原创性声明 本人郑重承诺：所呈交的毕业设计（论文），是我个人在指导教师的指导下进行的研究工作及取得的成果。尽我所知，除文中特别加以标注和致谢的地方外，不包含其他人或组织已经发表或公布过的研究成果，也不包含我为获得及其它教育机构的学位或学历而使用过的材料。对本研究提供过帮助和做出过贡献的个人或集体，均已在文中作了明确的说明并表示了谢意。 作者签名：日期： 指导教师签名：日期： 使用授权说明 本人完全了解大学关于收集、保存、使用毕业设计（论文）的规定，即：按照学校要求提交毕业设计（论文）的印刷本和电子版本；学校有权保存毕业设计（论文）的印刷本和电子版，并提供目录检索与阅览服务；学校可以采用影印、缩印、数字化或其它复制手段保存论文；在不以赢利为目的前提下，学校可以公布论文的部分或全部内容。 作者签名：日期：
摘要 电力系统潮流计算是研究电力系统稳态运行情况的一种计算，它根据给定的运行条件及系统接线情况确定整个电力系统各部分的运行状态：各母线的电压，各元件中流过的功率，系统的功率损耗等等。在电力系统规划的设计和现有电力系统运行方式的研究中，都需要利用潮流计算来定量地分析比较供电方案或运行方式的合理性、可靠性和经济性。因此潮流计算是研究电力系统的一种很重要和很基础的计算。由于电力系统结构及参数的一些特点，并且随着电力系统不断扩大，潮流问题的方程式阶数越来越高，对这样的方程式并不是任何数学方法都能保证给出正确答案的。这种情况成为促使电力系统计算人员不断寻求新的更可靠方法的重要因素。 本文旨在于研究潮流计算的牛顿—拉夫逊法的基本原理，在Matlab环境中实现牛顿—拉夫逊法潮流计算的数学模型，程序流程以及编制相应程序，并在程序中融合了节点优化编号和稀疏技术，以提高计算效率。最后用IEEE-3O节点标准测试系统验证所编程序。 关键词：潮流计算Newtom-Raphson法节点优化稀疏技术Matlab ABSTRACT Power flow calculation is fundanmental of analysis. Network reconfiguration,fault management,state estimator etc also need the data of electrial system power flow.There is important significance to develop power flow calculation in allusion to traits of distribution network. This paper introduces the principle of Newtom-Raphson algorithm, which is developed for calculation of power flow calculation ,where zero sequence network is open.With this algorithm,the three-phase load is resolved into positive/negative sequence power and coupling power,thus,decoupling three phase power flow into sequencet component power flow.The power flow can be obtained by just finding the positive sequence power flow and then finding the negative sequent component from the coupling https://www.doczj.com/doc/f110887948.html,pared with the existing methods,the jacobian matrix with the proposed algorithm is of much lower order,thus substantially reducing the computation burden.The proposed algorithm,together with a reference algorithm,has been simulated on an actual IEEE-30 system using statistic load date.And then it will