Matlab实验2
MATLAB实 验 报 告2
实验项目名称：晶闸管三相桥式整流器的仿真实验时间：
实验室指导教师：
一、实验目的
1、掌握晶闸管三相桥式整流器构成
2、掌握晶闸管三相桥式整流器的仿真模型
二、预习与思考
1、三相电源如何设置？
2、6脉冲触发器模块在哪个库中？
3、整流桥模块在哪个库中？
三、仿真步骤
通用桥臂模块参数设置
Number of bridge arms：桥臂数量，可以选择1、2、3相桥臂，构成不同形式的整流器。Port configuration：端口形式设置。Snubber resistance Rs(ohms):缓冲电阻Rs。Snubber capacitance Cs(F):缓冲电容Cs。Resistance Ron（ohms）：晶闸管的内电阻Ron，单位为Ω。Inductance Lon（H）：晶闸管的内电感Lon，单位为H，电感不能设置为0。Forward voltage Vf（v）:晶闸管元件的正向管压降Vf，单位为V。Measurements：测量可以选择5中形式，即无（None），装置电压（Device voltages），装置电流（Device UBC UCA UDC）或所有电压电流（All voltages and currents）
（2）同步脉冲触发器
同步脉冲触发器用于触发三相全控整流桥的6个晶闸管，同步6脉冲触发器可以给出双脉冲,双脉冲间隔为60°,触发器输出的1～6号脉冲依次送给三相全控整流桥对应编号的6个晶闸管.同步脉冲触发器包括同步电源和六脉冲触发器两个部分alpha_deg：此端子为脉冲触发角控制信号输入；AB, BC, CA：三相电源的三相线电压输入即Vab, Vbc, and VBlock：触发器控制端，输入为“0”时开放触发器，输入大于零时封锁触发器；
实验二MATLAB绘制图形
grid on %在所画出的图形坐标中加入栅格
绘制图形如下
50
10
1
0.8
40
10
0.6
0.4
30
10
0.2
0
1020
-0.2
-0.4
1010
-0.6
-0.8
0
10
-1
-2
0
2
-2
0
2
10
10
10
10
10
10
如果在图中不加栅格
程序如下：
clear x=logspace(-1,2);%在10^(-1)到10^2之间产生50个 对数等分的行向量 subplot(121); loglog(x,10*exp(x),'-p') subplot(122); semilogx(x,cos(10.^x))
（2）plot(x,y)： 基本格式，x和y可为向量或矩阵. 1. 如果x，y是同维向量，以x元素为横坐标，以y元素 为纵坐标绘图. 2. 如果x是向量，y是有一维与x元素数量相等的矩阵， 则以x为共同横坐标， y元素为纵坐标绘图，曲线数目 为y的另一维数. 3. 如果x，y是同维矩阵，则按列以x，y对应列元素为 横、纵坐标绘图，曲线数目等于矩阵列数.
y=2*exp(-0.5*x).*cos(4*pi*x);
2
plot(x,y)
1.5
1
0.5
0
-0.5
-1
-1.5
-2
0
1
2
3
4
5
6
7
例4 绘制曲线
t=(0:0.1:2*pi);
x=t.*sin(3*t);
y=t.*sin(t).*sin(t);
实验2利用MATLAB分析信号频谱及系统的频率特性
实验2利用MATLAB分析信号频谱及系统的频率特性引言：在信号处理和通信领域中，频谱分析是一项非常重要的技术。
频谱分析可以帮助我们了解信号的频率特性，包括频率成分和幅度。
MATLAB是一款功能强大的数学软件，提供了多种工具和函数用于信号处理和频谱分析。
本实验旨在通过MATLAB分析信号频谱及系统的频率特性，深入理解信号处理和频域分析的原理和应用。
实验步骤：1.生成一个信号并绘制其时域波形。
首先，我们可以使用MATLAB提供的函数生成一个信号。
例如，我们可以生成一个用正弦函数表示的周期信号。
```matlabt=0:0.001:1;%时间范围为0到1秒，采样率为1000Hzf=10;%信号频率为10Hzx = sin(2*pi*f*t); % 生成正弦信号plot(t,x) % 绘制信号的时域波形图title('Time domain waveform') % 添加标题```2.计算信号的频谱并绘制频谱图。
使用MATLAB中的FFT函数可以计算信号的频谱。
FFT函数将信号从时域转换为频域。
```matlabFs=1000;%采样率为1000HzL = length(x); % 信号长度NFFT = 2^nextpow2(L); % FFT长度X = fft(x,NFFT)/L; % 计算X(k)f = Fs/2*linspace(0,1,NFFT/2+1); % 计算频率轴plot(f,2*abs(X(1:NFFT/2+1))) % 绘制频谱图title('Frequency spectrum') % 添加标题```3.使用MATLAB分析系统的频率特性。
MATLAB提供了Signal Processing Toolbox，其中包含了分析系统频率特性的函数和工具。
```matlabHd = designfilt('lowpassfir', 'FilterOrder', 6,'CutoffFrequency', 0.3, 'SampleRate', Fs); % 设计一个低通滤波器fvtool(Hd) % 显示滤波器的频率响应``````matlab[W,F] = freqz(Hd); % 计算滤波器的频率响应plot(F,abs(W)) % 绘制滤波器的振幅响应title('Frequency response of lowpass filter') % 添加标题```实验结果：运行上述代码后，我们可以得到如下结果：1.时域波形图2.频谱图3.滤波器频率响应讨论与结论：本实验通过MATLAB分析信号频谱及系统的频率特性，深入理解了信号处理和频域分析的原理和应用。
matlab实验二
北京工业大学Matlab实验报告**: ***学号: ************: **实验二、Matlab 的基本计算（一）实验目的1.掌握建立矩阵的方法。
2.掌握Matlab 各种表达式的书写规则以及常用函数的使用。
3.能用Matlab 进行基本的数组、矩阵运算。
4.掌握矩阵分析的方法以及能用矩阵运算或求逆法解线性方程组。
5.掌握Matlab 中的关系运算与逻辑运算。
(二）实验环境1．计算机2．MATLAB7.0集成环境（三）实验内容及要求1、熟练操作MATLAB7.0运行环境；2、自主编写程序，必要时参考相关资料；3、实验前应写出程序大致框架或完整的程序代码；4、完成实验报告。
（四）实验程序设计1.利用diag 等函数产生下列矩阵。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=032570800a ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=804050702b2.利用reshape 函数将1题中的a 和b 变换成行向量。
3.产生一个均匀分布在（-5，5）之间的随机矩阵（10×2），要求精确到小数点后一位。
4.已知：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡-=76538773443412A ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=731203321B求下列表达式的值：(1) B A K *611+=和I B A K +-=12（其中I 为单位矩阵）(2) B A K *21=和B A K *.22=(3) 331^A K =和3.32^A K =(4) B A K /41=和A B K \42=(5) ],[51B A K =和]2:);],3,1([[52^B A K = 5.下面是一个线性方程组：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡52.067.095.03216/15/14/15/14/13/14/13/12/1x x x(1)求方程的解（矩阵除法和求逆法）(2)将方程右边向量元素3b 改为0.53，再求解，并比较3b 的变化和解的相对变化。
MATLAB实验二
3. 三维图形
• 三、视点
➢view(az,el)-------------设置观察点方向
az为方位角，el为仰角。 • 方位角为视点位置在XY平面上的投影与X轴形成的角度，
正值表示逆时针，负值表示顺时针。 • 仰角为XY平面的上仰或下俯角，正值表示视点在XY平
面上方，负值表示视点在XY平面下方。
figure grid gtext hold subplot text title xlabel ylabel
创建图形窗口 放置坐标网格线 用鼠标放置文本 保持当前图形窗口内容 创建子图 放置文本 放置图形标题 放置X轴坐标标记 放置Y轴坐标标记
3. 三维图形 • 一、 plot3函数
• 将二维函数plot的功能扩展到三维空间，绘制三维图形。 • 函数格式：plot3(x1,y1,z1,c1,x2,y2,z2,c2,…)
theta=[0:0.01:2*pi];
rho=sin(2*theta).*cos(2*theta);
polar(theta,rho); 绘制极坐标图命令
title('polar plot');
例9、程序：
theta=linspace(0, 2*pi);
r=cos(4*theta);
polar(theta, r);
所组成的画面。
4. 动画设计
• 【例14】 播放一个不断变化的眼球程序段。
• m=moviein(20); 建立一个20个列向量组成的矩阵
• for j=1:20
•
plot(fft(eye(j+10))) %绘制出每一幅眼球图并保存到m矩阵中
• m(:,j)=getframe;
MATLAB实验2 MATLAB的矩阵操作(1)
6、求特征值的方法求解方程
4x4+6x+3=0
(对比eig和roots、分析它们的不同和各自的应用范围)
7、P55思考练习3 , 5
实验过程中，请结合MATLAB的帮助系统。
四、实验操作过程与结果（写明使用到的函数及结算结果即可）
五、学生收获
MATLAB实验2MATLAB的矩阵操作(1)
班级专业：姓名：学号：日期：
一、实验目的
ห้องสมุดไป่ตู้1.熟悉MATLAB基本命令与操作；
2.熟悉MATLAB的矩阵运算；
3.了解MATLAB的多项式运算；
二、实验准备
通读书本第二章---MATLAB矩阵及其运算
三、实验内容
1.在命令窗口中键入表达式 ，并求 时 的值。
2. P54二.1
3.已知 ， ，在MATLAB命令窗口中建立A、B矩阵并对其进行以下操作：
(1)计算矩阵A的行列式的值
(2)分别计算下列各式：
、 和 、 、 、 、
4、产生3阶随机方阵，其元素为[10，60]区间的随机整数，然后判断A的元素是否能被5整除。（分析fix、rand、rem函数的使用方法）
实验二MATLAB的矩阵操作_参考答案
1
5
>> A(k)
ans =
23
10
（2）取出A前3行构成矩阵B,前两列构成矩阵C,右下角 子矩阵构成矩阵D,B与C的乘积构成矩阵E.
>> B=A([1,2,3],:)
B =
23.0000 10.0000 -0.7780 0
41.0000 -45.0000 65.0000 5.0000
32.0000 5.0000 0 32.0000
y =
-128.4271
2.已知 ,
求下列表达式的值：
（1） ， （其中I为单位矩阵）
A=[-1,5,-4;0,7,8;3,61,7]
B=[8,3,-1;2,5,3;-3,2,0]
>> A+6*B
ans =
47 23 -10
12 37 26
-15 73 7
&
ans =
-1.2768 -0.4743 0.2411
2.1229 1.3173 -0.2924
3.已知
完成下列操作
（1）输出A在[10,25]范围内的全部元素
A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14]
>> k=find(A>=10&A<=25)
(2)
(3)
2.已知 ,
求下列表达式的值：
（1） ， （其中I为单位矩阵）
(2)A*B、A.*B和B*A
(3)A/B及B/A
3.已知
完成下列操作
(1)输出A在[10,25]范围内的全部元素
(2)取出A前3行构成矩阵B,前两列构成矩阵C,右下角 子矩阵构成矩阵D,B与C的乘积构成矩阵E.
实验二MATLAB程序设计含实验报告
实验二 MATLAB 程序设计一、 实验目的1.掌握利用if 语句实现选择结构的方法。
2.掌握利用switch 语句实现多分支选择结构的方法。
3.掌握利用for 语句实现循环结构的方法。
4.掌握利用while 语句实现循环结构的方法。
5.掌握MATLAB 函数的编写及调试方法。
二、 实验的设备及条件计算机一台（带有MATLAB7.0以上的软件环境）。
M 文件的编写：启动MATLAB 后，点击File|New|M-File ，启动MATLAB 的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save 保存程序，注意文件名最好用英文字符。
点击Debug|Run 运行程序，在命令窗口查看运行结果，程序如有错误则改正三、 实验内容1.编写求解方程02=++c bx ax 的根的函数（这个方程不一定为一元二次方程，因c b a 、、的不同取值而定），这里应根据c b a 、、的不同取值分别处理，有输入参数提示，当0~,0,0===c b a 时应提示“为恒不等式!”。
并输入几组典型值加以检验。
（提示：提示输入使用input 函数）2.输入一个百分制成绩，要求输出成绩等级A+、A 、B 、C 、D 、E 。
其中100分为A+，90分～99分为A ，80分～89分为B ，70分～79分为C ，60分～69分为D ，60分以下为E 。
要求：（1）用switch 语句实现。
（2）输入百分制成绩后要判断该成绩的合理性，对不合理的成绩应输出出错信息。
（提示：注意单元矩阵的用法）3.数论中一个有趣的题目：任意一个正整数，若为偶数，则用2除之，若为奇数，则与3相乘再加上1。
重复此过程，最终得到的结果为1。
如：2?13?10?5?16?8?4?2?16?3?10?5?16?8?4?2?1运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。
请为关键的Matlab 语句填写上相关注释，说明其含义或功能。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
运行结果: Desired length = 10 Sampling frequency = 2000
2.编写MATLAB程序,利用rand产生一个均匀分 布的随机信号,并加到一实指数序列 e 0.8 n 上,计算此受干扰序列的自相关序列,证明受 干扰序列的自相关序列在零延时处呈现一个尖峰.
N = input('Length of sequence = '); n = 0:N-1; x = exp(-0.8*n); y = randn(1,N)+x; n1 = length(x)-1; r = conv(y,fliplr(y)); k = (-n1):n1; stem(k,r); xlabel('Lag index');ylabel('Amplitude'); gtext('r_{yy}[n]');
�
[y,fs,nbits] = wavread('oilyrag.wav'); len=length(y); x=1:len; plot(x,y); title('orial speech');
估计x1]=max(y1); xmax1=3100+xmax1-1; y2=y(3250:3350); [ymax2,xmax2]=max(y2); xmax2=3250+xmax2-1; Tvoice=xmax2-xmax1; T=Tvoice/fs; fvoice=fs/Tvoice;
1.产生下列序列,并用stem绘出其图形: ①单位取样序列 δ (n ) ②单位阶跃序列 (n ) ③斜变序列 n (n )
所需输入的数据是产生序列的长度L和抽样频率FT.
L = input('Desired length = '); n = 1:L; FT = input('Sampling frequency = ');T = 1/FT; imp = [1 zeros(1,L-1)];step = ones(1,L); ramp = (n-1).*step; subplot(3,1,1); stem(n-1,imp); xlabel(['Time in ',num2str(T), ' sec']);ylabel('Amplitude'); subplot(3,1,2); stem(n-1,step); xlabel(['Time in ',num2str(T), ' sec']);ylabel('Amplitude'); subplot(3,1,3); stem(n-1,ramp); xlabel(['Time in ',num2str(T), ' sec']);ylabel('Amplitude');
Length of sequence = 50
四,本实验用到的matlab命令 本实验用到的 命令
diap abs cos conv clf subplot filter impz wavread soundc
五,思考题: 思考题:
查阅语音的相关资料,回答什么是清音, 什么是浊音?读入一段语音,分别画出清 音和浊音的波形,并且根据波形,估测浊 音的基音.
一,实验目的:
1,学习MATLAB; 2,加深对常用离散信号及其基本运算的 理解; 3,常见离散信号的MATLAB产生和图形 显示.
二,实验内容: 实验内容:
1,典型序列的产生与显示; 2,序列的简单运算; 3,复杂信号的产生与显示; 4,语音信号的读入,回放,波形显示.
三,实验例程: 实验例程: