实验三 Matlab图像处理基本操作及摄像机标定(DLT)

实验三 Matlab图像处理基本操作及摄像机标定（D L T）

1、实验目的

通过应用Matlab的图像处理基本函数，学习图像处理中的一些基础操作和处理。

理解摄像机标定（DLT）方法的原理，并利用程序实现摄像机内参数和外参数的估计。

2、实验内容:

1)读取一幅图像并显示。

2)检查内存（数组）中的图像。

3)实现图像直方图均衡化。

4)读取图像中像素点的坐标值。

5)保存图像。

6)检查新生成文件的信息。

7)使用阈值操作将图像转换为二值图像。

8)根据RGB图像创建一幅灰度图像。

9)调节图像的对比度。

10)在同一个窗口内显示两幅图像。

11)掌握matlab命令及函数，获取标定块图像的特征点坐标。

12)根据摄像机标定（DLT）方法原理，编写Matlab程序，估计摄像机内参数和

外参数。

3、实验要求：

1)选取一幅图像，根据实验内容1）—10）给出结果。

2)根据给定的标定块图像及实验内容11），12）进行编程实验。

3)书写实验报告

4、实验设备

1)微机。

2)Matlab软件。

5、实验原理

DLT变换：

Abdal-Aziz和Karara于70年代初提出了直接线性变换像机定标的方法，他们从摄影测量学的角度深入的研究了像机图像和环境物体之间的关系，建立了像机成像几何的线性模型，这种线性模型参数的估计完全可以由线性方程的求解来实现。

直接线性变换是将像点和物点的成像几何关系在齐次坐标下写成透视投影矩阵的形式：

其中 为图像坐标系下的点的齐次坐标，

为世界坐标系下的空间点的欧氏坐标， P 为3*4的透视投影矩阵， 为未知尺度因子。

消去S ，可以得到方程组:

当已知N 个空间点和对应的图像上的点时，可以得到一个含有2*N 个方程的方程组：

其中A 为(2N*12)的矩阵， L 为透视投影矩阵元素组成的向量：

像机定标的任务就是寻找合适的L ，使得为 最小，即

给出约束：

L ‘为L 的前11个元素组成的向量， C 为A 前

11列组成的矩阵， B 为A 前12列组成的向量。

6、程序代码

I=imread ('DSCN0831.JPG'); imshow(I); whos

figure,

I1=rgb2gray(I);

||

|| min AL L

??

??

??

??????=????

???????1143w w

w Z Y X P v u s ()1,,v u ()w w w Z Y X ,,0034333231142322213433323114131211=----+++=----+++u p uZ p uY p uX p p Z p Y p X p u p uZ p uY p uX p p Z p Y p X p w w w w w w w w w w w w 0=AL []T p p p p p p p p p p p p 343332312423222114131211,,,,,,,,,,,||||AL 134=p B

C C C L T

T 1')(--=

figure,imshow(I1)

imhist(I1);

I2=histeq(I1)

figure,imshow(I2)

figure,imhist(I2)

a=imread('DSCN0831.JPG')

imwrite(I1,'DSgray.jpg');

imwrite(I2,'DSgrayeq.jpg');

inf=imfinfo('DSCN0831.JPG') ˉ

level=graythresh(I2);

bw=im2bw(I2,level);

figure,imshow(bw)

whos

I3=imadjust(I2,stretchlim(I2),[0 1]);

figure,imshow(I3);

a=imread('DSCN0831.JPG');

b=imread('DSgray.jpg');

subplot(1,2,1),imshow(a);

subplot(1,2,2),imshow(b);

clear;

xpot=[];

ypot=[];

imshow('DSCN0831.JPG');

hold on;

n=input('Please input the number of pot:');

for i=1:n

[x,y]=ginput(1);

plot(x,y,'or');

text(x+1,y+1,num2str(i));

xpot=[xpot,x];

ypot=[ypot,y];

end

[xpot;ypot]

[xw;yw;zw]

a=[xw',yw',zw',ones([100 1]),zeros([100

4]),(-1)*xpot'.*xw',(-1)*xpot'.*yw',(-1)*xpot'.*zw',-1*xpot';zeros([1

00 4]),xw',yw',zw',ones([100

1]),(-1)*ypot'.*xw',(-1)*ypot'.*yw',(-1)*ypot'.*zw',-1*ypot'];

c=a(:,1:11);

b=a(:,12);

l=(-1)*(c'*c)^(-1)*c'*b

7、运行结果

图像信息：

Filename: 'DSCN0831.JPG'

FileModDate: '21-十月-2011 21:29:34'

FileSize: 160501

Format: 'jpg'

FormatVersion: ''

Width: 1024

Height: 768

BitDepth: 24

ColorType: 'truecolor'

FormatSignature: ''

NumberOfSamples: 3

CodingMethod: 'Huffman'

CodingProcess: 'Sequential'

Comment: {}

ImageDescription: ' '

Make: 'NIKON '

Model: 'E990 '

Orientation: 1

XResolution: 300

YResolution: 300

ResolutionUnit: 'Inch'

Software: 'E990v1.0 '

DateTime: '0000:00:00 00:00:00 '

YCbCrPositioning: 'Co-sited'

DigitalCamera: [1x1 struct]

检查内存中的图像：

Name Size Bytes Class Attributes

I 768x1024x3 2359296 uint8 I1 768x1024 786432 uint8 I2 768x1024 786432 uint8 a 768x1024x3 2359296 uint8 bw 768x1024 786432 logical inf 1x1 12144 struct level 1x1 8 double 图像

图一直方图

图二均衡后的直方图

图三调节对比度后的图像

图四二进制图

图五灰度图

图六显示两幅图像

图七标定图

定块图像的特征点坐标

Columns 1 through 9

-0.0818 -0.3738 0.3232 0.8692 0.5916 0.1596 -0.5593 -1.1503 -0.9188

1.0132 1.2610 1.2829 0.8808 0.7317 0.8048 0.0446 0.3487 0.8531

Columns 10 through 18

-0.1365 -0.7079 -1.4787 -1.4430 -1.2126 -1.8520 -2.8819 -3.5305 -4.3289

0.4174 -0.1265 0.0022 0.5168 0.6806 1.8721 2.7697 3.0680 3.6265

Columns 19 through 27

-4.9902 -5.4902 -5.9902 -6.4902 -6.9672 -7.4556 -7.9556 -8.4556 -8.9556

4.1557 4.6557

5.1557 5.6557

6.1703 6.6703

7.1703 7.6703

8.1703

Columns 28 through 36

-9.3865 -9.3796 -9.2874 -9.2702 -9.2702 -9.2702 -9.2702 -9.2759 -9.2817

8.7069 9.2727 9.7580 10.2580 10.7580 11.2580 11.7580 12.2654 12.7873

Columns 37 through 45

-9.2874 -9.2874 -9.2874 -9.2874 -9.2874 -9.2874 -9.2874 -9.2874 -9.2874

13.2946 13.7946 14.2946 14.7946 15.2946 15.7946 16.2946 16.7946 17.2946

Columns 46 through 54

-9.2874 -9.2874 -9.2874 -9.2874 -9.1319 -9.0801 -8.9706 -8.5455 -8.8796

17.7946 18.2946 18.7946 19.2946 19.2902 19.2317 19.0855 18.8516 18.2054

Columns 55 through 63

-9.1734 -9.7483 -9.6031 -9.2978 -8.8381 -8.7298 -9.0697 -9.6849 -10.2137

17.8297 18.0095 18.8896 18.8487 18.4320 17.8589 17.4174 16.9174 16.4247

Columns 64 through 72

-10.8174 -11.5421 -12.8646 -13.8704 -14.8589 -15.5248 -16.0467 -16.1262 -15.5904

16.0197 16.1221 16.4803 16.9262 17.7332 18.6382 19.6089 20.4700 20.4437

Columns 73 through 81

-15.1653 -15.1492 -15.6503 -16.4153 -16.9153 -16.4591 -15.9591 -15.5167 -16.5593

20.0855 19.7273 19.2712 19.3341 19.3633 19.3560 19.3560 19.3487 19.7727

Columns 82 through 90

-17.5651 -18.4268 -18.9672 -18.9315 -18.6146 -19.7033 -20.3520 -20.8520 -20.3727

20.1776 19.8224 19.0665 18.3253 17.9905 17.8414 17.5636 17.5636 17.5636

Columns 91 through 99

-19.7402 -18.9061 -17.9349 -17.4706 -17.1377 -17.3796 -17.6734 -18.5132 -19.4614

17.5782 17.6294 17.9437 18.5285 19.6893 20.7478 21.7332 22.1849 22.3268

Column 100

-19.8347

22.1148

ans =

1.0e+003 *

Columns 1 through 9

0.5020 0.5050 0.5065 0.5125 0.5080 0.5155 0.5125 0.5140 0.4480

0.7212 0.6492 0.5683 0.4797 0.3957 0.2982 0.2022 0.1047 0.6943

Columns 10 through 18

0.4465 0.4465 0.4465 0.4465 0.4465 0.4465 0.4450 0.3880 0.3370

0.6178 0.5383 0.4557 0.3717 0.2832 0.1857 0.0957 0.0867 0.0837

Columns 19 through 27

0.2875 0.2455 0.2035 0.1705 0.1735 0.2065 0.2455 0.2890 0.3385

0.0747 0.0702 0.0627 0.0642 0.1332 0.1422 0.1482 0.1572 0.1662

Columns 28 through 36

0.3880 0.3910 0.3400 0.2935 0.2515 0.2110 0.1795 0.1780 0.2155

0.1797 0.2652 0.2532 0.2382 0.2292 0.2172 0.2082 0.2757 0.2892

Columns 37 through 45

0.2530 0.2950 0.3415 0.3910 0.3925 0.3415 0.2965 0.2560 0.2170

0.3027 0.3207 0.3357 0.3552 0.4347 0.4152 0.3957 0.3792 0.3582

Columns 46 through 54

0.1825 0.1885 0.2200 0.2590 0.2980 0.3445 0.3955 0.3940 0.3475

0.3447 0.4092 0.4302 0.4437 0.4662 0.4902 0.5112 0.5893 0.5623

Columns 55 through 63

0.3025 0.2635 0.2230 0.1900 0.2275 0.2680 0.3040 0.3475 0.3955

0.5367 0.5172 0.4962 0.4767 0.5577 0.5817 0.6042 0.6312 0.6597

Columns 64 through 72

0.5665 0.5725 0.5725 0.5740 0.5755 0.5800 0.5815 0.5860 0.6490

0.6943 0.6252 0.5473 0.4617 0.3762 0.2862 0.1932 0.0987 0.0927

Columns 73 through 81

0.6475 0.6430 0.6400 0.6370 0.6340 0.6280 0.6265 0.6820 0.6850

0.1827 0.2727 0.3627 0.4422 0.5232 0.5998 0.6687 0.6478 0.5757

Columns 82 through 90

0.6880 0.6940 0.6970 0.7015 0.7045 0.7135

1.1485 0.7285 0.7360

0.5037 0.4227 0.3447 0.2622 0.1767 0.0867

-0.2478 0.6222 0.5533

Columns 91 through 99

0.7390 0.7465 0.7495 0.7555 0.8065 0.8560

0.8350 0.8335 0.8455

0.4842 0.4122 0.3327 0.2517 0.2442 0.1587

0.3807 0.4452 0.3102

Column 100

0.8215

0.5787

转换为世界坐标

ans =

Columns 1 through 17

60 60 60 60 100 100 100 100 100 120 140 140 140 140 140 140 140

0 0 0 0 0 0 0 0 0 0 0 0 20 20 20 20 60

60 80 100 120 100 80 60 40 20 80 120 60 120 100 80 40 80

Columns 18 through 34

140 140 140 140 140 140 140 140 20 20 20 20 60 60 60 60 60

60 60 100 100 120 120 120 140 0 40 80 120 120 100 80 60 40

60 40 60 120 120 100 60 40 140 140 140 140 140 140 140 140 140

Columns 35 through 50

100 100 120 120 120 120 140 140 140 140 140

100 80 60 40 20

80 0 100 80 60 40 20 40 60 100 120 140 140 140 140 140

140 140 140 140 140 140 140 140 140 140 140 140 140 140 140 140

实验一 MATLAB基本操作及运算(含实验报告).

实验一 MATLAB 基本操作及运算 一、 实验目的 1、 理解Matlab 数据对象的特点； 2、 掌握基本Matlab 运算规则； 3、 掌握Matlab 帮助的使用方法； 二、 实验的设备及条件 计算机一台（带有MATLAB7.0以上的软件环境）。 三、 实验内容 要求建立一个名为experiment01.m 的，把与实验内容1-7相关的实验命令都放入该文件中，题与题之间用相应注释分割。注意对实验中出现的相关函数或变量，请使用help 或doc 查询相关帮助文档，学习函数的用法。 1、 建立以下标量： 1） a=10 2） b=2.5×1023 3） c=2+3i ，（i 为虚数单位） 4） d=3/2πj e ,（j 为虚数单位，这里要用到exp ，pi ） 2、 建立以下向量： 1） aVec=[3.14 15 9 26] 2） bVec=????? ???????18228871.2 3） cVec=[5 4.8 … -4.8 -5 ] （向量中的数值从5到-5，步长为-0.2） 4） dVec=[100 100.01 … 100.99 101] （产生1到10之间的等对数间隔向量，参考logspace ，注意向量的长度） 3、 建立以下矩阵： 1）???? ??????=2222 aMat aMat 一个9×9的矩阵，其元素全为2；（参考ones 或zeros ）

2）??????? ?????????=1000005000001 bMat bMat 是一个9×9的矩阵，除主对角上的元素为[1 2 3 4 5 4 3 2 1]外，其余元素均为0。（参考diag ）。 3）100 20109212291111 =cMat cMat 为一个10×10的矩阵，可有1:100的向量来产生（参考reshape ） 4）???? ??????=NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN dMat dMat 为3×4的NaN 矩阵，（参考nan ） 5）?? ????---=8710225113eMat 6）产生一个5×3随机整数矩阵fMat ，其值的范围在-3到3之间。（参考rand 和floor 或ceil ） 4、 使用题1中的变量计算下列等式的x,y,z 的值： 1） ) 6/)15((11--+=a e x 2） g g h h b a y /121,)(=+=提示π,参考sqrt 。 3） c c a d c d c R z ))3/sin()]))([(log(π-+= ，其中R 表示取括号内复数的实数部分，c 表示c 的共轭复数，log 是自然对数。（参考real ，conj ，log ） 5、 使用题2中的向量求解一下等式： 1）)25.2/(22 25.221 cVec e xVec -=π， 其中cVec 指的是题2 中定义的向量cVec ，一下雷同。 2）22)(bVec aVec yVec T +=,T aVec 表示aVec 的转置 3） )/1(log 10dVec zVec =，10log 表示已10为底的对数，参考log10 6、 使用题2和题3中所产生的向量和矩阵计算以下等式，注意本题的操作

基于MATLAB平台的可视化图像处理系统设计方法

[收稿日期]20050228 [作者简介]黄书先(1962),女,1983年大学毕业,硕士,副教授,现主要从事石油勘探开发研究与科研管理工作。 基于MATLAB 平台的可视化图像处理系统 设计方法 黄书先 (长江大学科学技术处,湖北荆州434023) [摘要]以MATLAB 为工作语言和开发环境,开发了一个在M AT LAB 平台下的可视化图像处理系统, 可实现包括对测井图像在内的一般图像的精细处理,并能和用户开发的程序接口。为M AT LAB 的再开发 和可视化系统的设计作了有益的探索。 [关键词]MATLAB;图像处理;可视化GUI 界面 [中图分类号]TP 311111 [文献标识码]A [文章编号]16731409(2005)04015803 MAT LAB 的图像处理工具为自然科学各学科领域的学者、研究人员和工程师提供了一个直观的灵活的环境,用以解决复杂的图像处理问题。用MAT LAB 语言开发的图像处理算法可以在所有支持MAT LAB 的平台上共享。也可以将m 语言算法和现存的C 程序集成在一起或者将MATLAB 开发的m 语言算法和GU Is 编译为C/C++代码,供其他程序调用,或者发布为一个独立的应用程序。下面笔者提出一个基于MATLAB 平台的可视化图像处理系统设计方法,可用于包括测井图像资料在内的一般图像的精细处理。 1 系统总体设计 在MAT LAB 中有个重要的图像处理工具包[1],该工具包是由一系列支持图像处理操作的函数组成的,按功能可以分为以下几类:图像显示;图像文件输入与输出;几何操作;像素值统计;图像增强;图像识别;图像滤波;图像变换;邻域和块操作;二值图像操作;颜色映射和颜色空间转换;图像格式转换等。和其他工具包一样,用户还可以根据需要书写自己的函数,以满足特定的需要;也可以将这个工具包和信号处理工具包或小波工具包等其他工具包联合起来使用。 MAT LAB 提供了交互式的GU I 开发环境[2~4],用户只需要设置各个对象相应的属性,系统自动生成与之对应的界面,大大减少了开发的难度。本设计的MAT LAB 图像处理系统由封面、主界面和各个子功能界面组成,其框图如图1所示。程序总流程图如图2所示。 2 封面界面设计 首先用MAT LAB 编辑封面的脚本文件,生成系统的封面界面,然后再进入处理的主界面。 在设计封面时,要做好封面的总体布局,力求完美。要插入背景,首先要找到所用的函数(这在一般的书中很难找到);其次要注意设置axes 的属性。在显示背景图像时,不能用imshow ()函数,如果用此函数,背景图像只能显示在封面的一部分,不能整屏覆盖;而要用imagesc ()函数。MATLAB 提供修改文本中文字的字体,它支持华文中宋、华文彩云、华文仿宋、楷书、黑体等,功能较强。 在设计封面的过程中,最主要的技术是要解决时间的显示,如果只用MATLAB 中的clock ()函数,则显示的是静态时间,而不会显示和电脑同步的动态时间。为了解决这个问题,需做一个循环判断语句:while find (get (0,'c hildren'))==h0。成立的条件是,只要是当前窗口循环,否则停止,这样可以减少CPU 的负担。同时,要实时提取clock ()函数,可用fix (clock),使提取的时间更美观。#158#长江大学学报(自科版)2005年4月第2卷第4期/理工卷第2卷第2期 Jour nal of Yangtze University (Nat Sci Edit)Apr 12005Vol 12No 14/Sci &Eng V,Vol 12No 12

实验1_基于MATLAB的图像基本操作

第1次实验基于MATLAB的图像基本操作 二、实验内容和要求： 1．实现图像Baboon.bmp（MATLAB自带）的读入（可使用imread）和显示（可使用imshow）操作，代码加上足够的注释，需要建立一个M文件实现。 I=imread('F:\标准图像\Baboon.bmp');//读入图像 imshow(I);//显示图像 2．编程实现将一幅RGB图像转换为二值图像，并在一个窗口同时显示处理过程中得到的每一个图像和原图像，同时需要给图像加上标题。（原始数据可以是任意的RGB图像）。需要新建一个M文件实现。 figure,subplot(1,3,1),imshow(I(:,:,1)),title('R'); subplot(1,3,2),imshow(I(:,:,2)),title('G'); subplot(1,3,3),imshow(I(:,:,3)),title('B'); 3.计算图象统计参数： 读取图像（文件名为‘cameraman.tif’）; 最大值 最小值 均值 K=imread('cameraman.tif'); d_max=max(K(:)) d_min=min(K(:)) d_mean=mean(K(:)) 4.利用帮助系统了解im2double,imresize,image函数的作用和语法，并利用这些函数处理已知图像pout.tif（MATLAB自带）并显示处理前后效果。 J=imread('pout.tif'); J1=im2double(J); figure,subplot(1,2,1),imshow(J),title('Before') subplot(1,2,2),imshow(J1),title('After') J2=imresize(J,0.3); figure,subplot(1,2,1),imshow(J),title('Before') subplot(1,2,2),imshow(J2),title('After') figure,subplot(1,2,1),imshow(J),title('Before') subplot(1,2,2),image(J);title('After') 1

实验室常用器材使用方法及注意事项

实验室常用器材使用方法及注意事项

实验室常见仪器使用方法及注意事项 一、常见的仪器 （一）初中化学实验常见仪器 反应容器可直接受热的：试管、蒸发皿、燃烧匙、坩埚等能间接受热的：烧杯、烧瓶、锥形瓶（加热时，需加石棉网） 常存放药品的仪器：广口瓶（固体）、细口瓶（液体）、滴瓶 （少量液体）、集气瓶（气体） 用加热仪器：酒精灯 计量仪器：托盘天平（称固体质量）、量筒（量液体体积） 仪分离仪器：漏斗 取用仪器：药匙（粉末或小晶粒状）、镊子（块状或较大颗粒）、胶头滴管（少量液体） 器夹持仪器：试管夹、铁架台（带铁夹、铁圈）、坩埚钳其它仪器：长颈漏斗、石棉网、玻璃棒、试管刷、水槽 不能加热：量筒、集气瓶、漏斗、温度计、滴瓶、表面皿、广口瓶、细口瓶等 1、试管 （1）、用途： a、在常温或加热时，用作少量试剂的反应容器。 b、溶解少量固体。 c、收集少量气体的容器 d、用于装置成小型气体的发生

器。 （2）、注意事项： a、加热时外壁必须干燥，不能骤热骤冷，一般要先均匀受热，然后才能集中受热， 防止试管受热不均而破裂。 b、加热时，试管要先用铁夹夹持固定在铁架台上（短时间加热也可用试管夹夹持）。 试管夹应夹在的中上部（或铁夹应夹在离试管口的1/3处）。c、加热固体时，试管口要略向下倾斜，且未冷前试管不能直立，避免管口冷凝水倒流 使试管炸裂。 d、加热液体时，盛液量一般不超过试管容积的1/3（防止液体受热溢出），使试管与桌面 约成45°的角度（增大受热面积，防止暴沸），管口不能对着自己或别人（防止液体喷出伤人）。反应时试管内的液体不超过试管容积的1/2。 2、烧杯用途：①溶解固体物质、配制溶液，以及溶液的稀释、浓缩 ②也可用做较大量的物质间的反应 注意事项：受热时外壁要干燥，并放在石棉网上使其受热均匀（防止受热不均使烧杯炸裂）， 加液量一般不超过容积的1/3（防止加热沸腾使液体外溢）。

实验二 Matlab程序设计基本方法1

实验二Matlab程序设计基本方法 覃照乘自092 电气工程学院 一、实验目的： 1、熟悉MATLAB 程序编辑与设计环境 2、掌握各种编程语句语法规则及程序设计方法 3、函数文件的编写和设计 4、了解和熟悉跨空间变量传递和赋值 二、实验基本知识： ◆for循环结构 语法：for i＝初值：增量：终值 语句1 …… 语句n end 说明：1．i＝初值：终值，则增量为1。 2．初值、增量、终值可正可负，可以是整数，也可以是小数，只须符合数学逻辑。 ◆while 循环结构 语法：while 逻辑表达式 循环体语句 end 说明：1、whiIe结构依据逻辑表达式的值判断是否执行循环体语勾。若表达式的值为真，执行循环体语句一次、在反复执行时，每次都要进行判断。若表达 式的值为假，则程序执行end之后的语句。 2、为了避免因逻辑上的失误，而陷入死循环，建议在循环体语句的适当位 置加break语句、以便程序能正常执行。（执行循环体的次数不确定； 每一次执行循环体后，一定会改变while后面所跟关系式的值。） 3、while循环也可以嵌套、其结构如下：

while逻辑表达式1 循环体语句1 while逻辑表达式2 循环体语句2 end 循环体语句3 end ◆if-else-end分支结构 if 表达式1 语句1 else if 表达式2(可选) 语句2 else(可选) 语句3 end end 说明：1.if结构是一个条件分支语句，若满足表达式的条件，则往下执行；若不满足，则跳出if结构。 2．else if表达式2与else为可选项，这两条语句可依据具体情况取舍。 3．注意：每一个if都对应一个end，即有几个if，记就应有几个end。 ◆switch-case结构 语法：switch表达式 case常量表达式1 语句组1 case常量表达式2 语句组2 …… otherwise 语句组n end

实验一 Matlab基本操作

实验一Matlab基本操作 题目： 1．利用基本矩阵产生 3x3 和15x8 的单位阵，全1 阵，全0 阵，均匀分布的随 机阵（[-1，1]之间），正态分布随机阵（方差4，均值1） 2．利用diag()函数和rot90()产生下列矩阵： 然后求解a 阵的逆矩阵aa 及b 阵的特征值和对应特征向量，并利用reshape 将 aa 阵变换成行向量。 3．产生一均匀分布在(-5,5)随机阵(50x2)，精确到小数点后一位。 4．编程实现当α∈[-π，π]，间隔为1o 时，求解正弦和余弦的值，并利用plot() 函数绘制正弦，余弦曲线。 5．利用rand 函数产生(0,1)间均匀分布的10x10 随机矩阵a,然后统计a 中大于等于0.6 的元素个数。 6．利用randn 函数产生均值为0，方差为1 的10x10 正态分布随机阵，然后统计其中大于-0.5，小于0.5 的元素个数。 7．编程实现下表功能： 8．有一矩阵a,找出矩阵中其值大于1 的元素，并将他们重新排列成列向量b。 9．在一保定市区9 月份平均气温变化测量矩阵temp_Baoding_sep 中(48x30)，存在有奇异值（大于42o C，小于0o C）,编程实现删除奇异值所在的行。 10．在给定的100x100 矩阵中，删除整行内容全为0 的行，删除整列内容全为0 的列。 程序： 1. %3X3矩阵 a1=eye(3) a2=ones(3) a3=zeros(3) a4=1-2*rand(3) a5=2*randn(3)+1 %15X8矩阵 b1=eye(15,8) b2=ones(15,8) b3=zeros(15,8) b4=1-2*rand(15,8) b5=2*randn(15,8)+1 运行结果：

常用化学仪器及使用方法、化学实验基本操作

常用化学仪器及使用方法、化学实验基本操作 知识分析： 包括的内容有1. 常用化学仪器及使用方法；2. 化学实验基本操作； 3. 物质的分离、提纯； 4. 常见气体的制备； 5. 物质的检验； 6. 常见试剂的保存； 7. 综合实验—实验方案的设计和评价。 在实验中复习对于基本实验要熟练掌握原理和基本操作，在此基础上，顺利地完成实验方案的设计。 常用化学仪器及使用方法 （一）1. 能直接加热的仪器 仪器图形与名称主要用途使用方法和注意事项 用于蒸发溶剂或浓缩溶液可直接加热，但不能骤冷。蒸发溶液时不可加得太满，液面应距边缘1厘米处。 常用作反应器，也可收集少量气体可直接加热，拿取试管时，用中指、食指、拇指拿住试管口占全长的1／3处，加热时管口不能对着人。放在试管内的液体不超过容积的1／2，加热的不超过1／3。加热时要用试管夹，并使试管跟桌面成45°的角度，先给液体全部加热，然后在液体底部加热，并不断摇动。给固体加热时，试管要横放，管口略向下倾。 用于灼烧固体，使其反应（如分解）可直接加热至高温。灼烧时应放于泥三角上，应用坩埚钳夹取。应避免聚冷。 燃烧少量固体物质可直接用于加热，遇能与Cu、Fe反应的物质时要在匙内铺细砂或垫石棉绒。 2. 能间接加热（需垫石棉网） 仪器图形和名称主要用途使用方法和注意事项 （分为50、100、250、500、1000ml等规格）用作配制、浓缩、稀释溶液。也可用 作反应器和给试管水浴加热等。 加热时应垫石棉网 根据液体体积选用不同规格烧杯 用作反应器（特别是不需加热的） 不能直接加热，加热时要垫石棉网。不 适于长时间加热，当瓶内液体过少时， 加热容易使之破裂。 用作在加热条件下进行的反应器 不能直接加热，应垫石棉网加热。所装 液体的量不应超过其容积1／2。 用于蒸馏与分馏，也可用作气体发生 器 加热时要垫石棉网。也可用其他热浴。 用作接受器 用作反应器，常用于滴定操作 一般放在石棉网上加热。在滴定操作中 液体不易溅出。 3. 不能加热的仪器 仪器图形与名称主要用途使用方法及注意事项 用于收集和贮存少量气体上口为平面磨砂，内侧不磨砂，玻璃片 要涂凡士林油，以免漏气，如果在其中 进行燃烧反应且有固体生成时，应在底 部加少量水或细砂。 分装各种试剂，需要避光保存时用 棕色瓶。广口瓶盛放固体，细口瓶 盛放液体。 瓶口内侧磨砂，且与瓶塞一一对应，切 不可盖错。玻璃塞不可盛放强碱，滴瓶 内不可久置强氧化剂等。 制取某些气体的反应器固体＋液 体 固体为块状，气体溶解性小反应无强热 放出，旋转导气管活塞控制反应进行或 停止。 （二）计量仪器 仪器图形与名称主要用途使用方法及注意事项

matlab实验二

实验2 MATLAB数值计算、符号运算功能 一、实验目的 1、掌握建立矩阵、矩阵分析与处理的方法。 2、掌握线性方程组的求解方法。 3、掌握数据统计和分析方法、多项式的常用运算。 4、掌握求数值导数和数值积分、常微分方程数值求解、非线性代数方程数值求解的方法。 5、掌握定义符号对象的方法、符号表达式的运算法则及符号矩阵运算、符号函数极限及导数、符号函数定积分和不定积分的方法。 二、预习要求 （1）复习4、5、6章所讲内容； （2）熟悉MATLAB中的数值计算和符号运算的实现方法和主要函数。 三、实验内容 1、已知 29618 20512 885 A -?? ?? =?? ?? - ?? ，求A的特征值及特征向量，并分析其数学意义。 >> A=[-29,6,18;20,5,12;-8,8,5]; >> [V,D]=eig(A) V = 0.7130 0.2803 0.2733 -0.6084 -0.7867 0.8725 0.3487 0.5501 0.4050 D = -25.3169 0 0 0 -10.5182 0 0 0 16.8351 V为A的特征向量，D为A的特征值，3个特征值是-25.3169、10.5182和16.8351。 >> A*V ans = -18.0503 -2.9487 4.6007 15.4017 8.2743 14.6886 -8.8273 -5.7857 6.8190 >> V*D

ans = -18.0503 -2.9487 4.6007 15.4017 8.2743 14.6886 -8.8273 -5.7857 6.8190 经过计算，A*V=V*D 。 2、 不用rot90函数，实现方阵左旋90°或右旋90°的功能。例如，原矩阵为A ，A 左旋后得到B ，右旋后得到C 。 147102581136912A ????=??????，101112789456123B ??????=??????，321654987121110B ??????=?????? 提示：先将A 转置，再作上下翻转，则完成左旋90°；如将A 转置后作左右翻转，则完成右旋转90°，可用flipud 、fliplr 函数。 >> a=[1 4 7 10;2 5 8 11;3 6 9 12] a= 1 4 7 10 2 5 8 11 3 6 9 12 >> B=rot90(a) B = 10 11 12 7 8 9 4 5 6 1 2 3 >>C= rot90(s,3) C= 3 2 1 6 5 4 9 8 7 12 11 10

matlab操作实验报告

实验一matlab基本操作 一、实验目的 熟悉matlab的安装与启动；熟悉matlab用户界面；熟悉matlab功能、建模元素；熟悉matlab优化建模过程。 二、实验设备与工具 1.计算机 2.matlab软件 三、实验步骤 1. 了解matlab的硬件和软件必备环境； 2. 启动matlab； 3. 学习优化建模过程。 四、实验报告要求 1. 写出matlab系统界面的各个构成；以及系统布局区的组成；以及每一部 分的功能； 2. 优化建模过程应用举例 五、实验内容 （一）、Matlab操作界面 1.命令窗口（command window） 2.命令历史窗口（command history） 3.工作空间管理窗口（workspace） 4.当前路径窗口（current directory） （二）、优化建模过程应用举例 1、简单矩阵 123 456 789 A ?? ?? =?? ?? ?? 的输入步骤。 （1）在键盘上输入下列内容 A = [1,1,3; 4,5,6; 7,8,9] （2）按【Enter】键，指令被执行。 （3）在指令执行后，MATLAB指令窗中将显示以下结果： A = 1 2 3 4 5 6 7 8 9

2、矩阵的分行输入。 A=[1,2,3 4,5,6 7,8,9] A = 1 2 3 4 5 6 7 8 9 3、指令的续行输入 S=1-1/2+1/3-1/4+ ... 1/5-1/6+1/7-1/8 S = 0.6345 4、画出衰减振荡曲线t e y t 3sin 3-=及其它的包络线3 0t e y -=。t 的取值范围是]4,0[π。 t=0:pi/50:4*pi; y0=exp(-t/3); y=exp(-t/3).*sin(3*t); plot(t,y,'-r',t,y0,':b',t,-y0,':b') grid 5、画出2222) sin(y x y x z ++=所表示的三维曲面。y x ,的取值范围是]8,8[-。 clear;x=-8:0.5:8; y=x';

化学实验室基本操作

化学实验室基本操作 化学实验室基本操作2010-07-24 11：54一、常用仪器的主要用途和使用 方法 反应容器：试管、燃烧匙、烧杯、锥形瓶、集气瓶 存放容器：集气瓶(气体)、细口瓶(液体)、广口瓶(固体)、滴瓶(少量液体) 计量仪器：托盘天平(称固体质量)、量筒(量液体体积) 取用仪器：镊子(块状或较大颗粒)、药匙或纸槽(粉末或小颗粒)、胶头滴 管(少量液体) 夹持容器：试管夹、坩埚钳、铁架台(带铁圈、铁夹) 其它仪器：漏斗、长颈漏斗、分液漏斗、石棉网、玻璃棒、水槽、试管刷 可直接加热的：试管、蒸发皿、燃烧匙 能间接加热的(需垫石棉网)：烧杯、烧瓶、锥形瓶 加热仪器：酒精灯 1.烧杯圆柱状玻璃容器，杯口有便于倒出液体的嘴。 常用的有25mL、50mL、100 mL、250 mL、500 mL等 (1)用于大量物质的溶解和配制溶液或者进行化学反应的容器，也常用于接 过滤后的液体。 (2)实验时盛放液体的量不超过烧杯容积的1/2，以防搅拌时溅出。 (3)向烧杯中注入液体的时候，应沿烧杯内壁或玻璃棒引流。

(4)加热时要垫石棉网，也防受热不均而使其破裂。烧杯不能用作加热固体试剂。 2.试管 (1)用于少量物质的溶解或发生化学反应的仪器，也常用于制取或收集少量气体。 (2)振荡试管的方法：手持试管、手腕摆动。 3)实验时盛放液体量不能超过试管容积的1/3，以防振荡或加热时溅出。可直接加热。 (4)用试管夹或者铁夹固定时，要从试管底部向上套，夹持在试管的中上部(或离管口1/3的部位)。 3.蒸发皿 (1)用于溶液的蒸发、结晶 2)蒸发过程中需用玻璃棒不断搅拌，防止液体由于局部温度过高而飞溅 3)当溶液的量减少只有大量晶体析出时，停止加热并放至石棉网上，以防晶体飞溅 (4)取放蒸发皿，要用坩埚钳夹持 4.集气瓶 (1)用于收集气体、短时间贮存气体、用做物质在气体中的燃烧的反应器 (2)在收集气体或贮存气体时，要用毛玻璃片盖住瓶口。 5、试剂瓶 试剂瓶包括滴瓶、细口瓶、广口瓶等。分为无色和棕色两种。

MATLAB实验报告实验二

实验二 MATLAB矩阵及其运算 学号：3121003104 姓名：刘艳琳专业：电子信息工程1班日期：2014.9.20 一实验目的 1、掌握Matlab数据对象的特点以及数据的运算规则。 2、掌握Matlab中建立矩阵的方法以及矩阵处理的方法。 3、掌握Matlab分析的方法。 二实验环境 PC_Windows 7旗舰版、MATLAB 7.10 三实验内容 4、1. （1）新建一个.m文件，验证书本第15页例2-1； （2）用命令方式查看和保存代码中的所有变量；

（3）用命令方式删除所有变量； （4）用命令方式载入变量z。 2. 将x=[4/3 1.2345e-6]在以下格式符下输出：短格式、短格式e方式、长格式、长格式e方式、银行格式、十六进制格式、+格式。 短格式 短格式e 长格式

长格式e方式 银行格式 十六进制格式 3.计算下列表达式的值 （1）w=sqrt(2)*(1+0.34245*10^(-6)) （2）x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a) a=3.5;b=5;c=-9.8; （3）y=2*pi*a^2*((1-pi/4)*b-(0.8333-pi/4)*a) a=3.32;b=-7.9; （4）z=0.5*exp(2*t)*log(t+sqrt(1+t*t)) t=[2,1-3i;5,-0.65];

4. 已知A=[1 2 3 4 5 ;6 7 8 9 10;11 12 13 14 15;16 17 18 19 20]，对其进行如下操作：（1）输出A在[ 7, 10]范围内的全部元素； （2）取出A的第2，4行和第1，3，5列； （3）对矩阵A变换成向量B，B=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20]； （4）删除A的第2，3，4行元素； (1) (2)

MATLAB基本操作实验报告

MATLAB基本操作 实验报告 课程名称： 院系： 专业班级： 学号： 学生姓名： 指导教师： 开课时间：至学年第学期

一、学生撰写要求 按照实验课程培养方案的要求，每门实验课程中的每一个实验项目完成后，每位参加实验的学生均须在实验教师规定的时间内独立完成一份实验报告，不得抄袭，不得缺交。 学生撰写实验报告时应严格按照本实验报告规定的内容和要求填写。字迹工整，文字简练，数据齐全，图表规范，计算正确，分析充分、具体、定量。 二、教师评阅与装订要求 1.实验报告批改要深入细致，批改过程中要发现和纠正学生实验报告中的问题，给出评语和实验报告成绩，签名并注明批改日期。实验报告批改完成后，应采用适当的形式将学生实验报告中存在的问题及时反馈给学生。 2.实验报告成绩用百分制评定，并给出成绩评定的依据或评分标准（附于实验报告成绩登记表后）。对迟交实验报告的学生要酌情扣分，对缺交和抄袭实验报告的学生应及时批评教育，并对该次实验报告的分数以零分处理。对单独设课的实验课程，如学生抄袭或缺交实验报告达该课程全学期实验报告总次数三分之一以上，不得同意其参加本课程的考核。 3.各实验项目的实验报告成绩登记在实验报告成绩登记表中。本学期实验项目全部完成后，给定实验报告综合成绩。 4.实验报告综合成绩应按课程教学大纲规定比例（一般为10-15%）计入实验课总评成绩；实验总评成绩原则上应包括考勤、实验报告、考核（操作、理论）等多方面成绩； 5.实验教师每学期负责对拟存档的学生实验报告按课程、学生收齐并装订，按如下顺序装订成册：实验报告封面、实验报告成绩登记表、实验报告成绩评定依据、实验报告（按教学进度表规定的实验项目顺序排序）。装订时统一靠左侧按“两钉三等分”原则装订。

MATLAB实验二(修改)

实验二 信号的表示及其基本运算 一、实验目的 1、掌握连续信号及其MATLAB 实现方法； 2、掌握离散信号及其MA TLAB 实现方法 3、掌握离散信号的基本运算方法，以及MA TLAB 实现 4 熟悉应用MATLAB 实现求解系统响应的方法 4、了解离散傅里叶变换的MA TLAB 实现 5、了解IIR 数字滤波器设计 6、了解FIR 数字滤波器设计1 二、实验设备 计算机，Matlab 软件 三、实验内容 （一）、 连续信号及其MATLAB 实现 1、 单位冲激信号 ()0,0()1,0 t t t dt ε ε δδε-?=≠??=?>??? 例1.1：单位冲击信号的MATLAB 实现程序如下： t1=-4; t2=4; t0=0; dt=0.01; t=t1:dt:t2; n=length(t); x=zeros(1,n); x(1,(-t0-t1)/dt+1)=1/dt; stairs(t,x); axis([t1,t2,0,1.2/dt]); 2、 任意函数 ()()()f t f t d τδττ+∞ -∞ =-? 例1.2：用MA TLAB 画出如下表达式的脉冲序列 ()0.4(2)0.8(1) 1.2() 1.5(1) 1.0(2)0.7(3)f n n n n n n n δδδδδδ=-+-+++++++ 3 单位阶跃函数 1,0()0, t u t t ?≥?=?

t=-0.5:0.001:1; t0=0; u=stepfun(t,t0); plot(t,u) axis([-0.5 1 -0.2 1.2]) 4 斜坡函数 0()()g t B t t =- 例1.4：用MA TLAB 实现g(t)=3(t-1) clear all; t=0:0.01:3; B=3; t0=1; u=stepfun(t,t0); n=length(t); for i=1:n u(i)=B*u(i)*(t(i)-t0); end plot(t,u) axis([-0.2 3.1 -0.2 6.2]) 5 抽样信号 抽样信号Sa(t)=sin(t)/t 在MATLAB 中用 sinc 函数表示。 定义为 )/(sin )(πt c t Sa = t=-3*pi:pi/100:3*pi; ft=sinc(t/pi); plot(t,ft); grid on; axis([-10,10,-0.5,1.2]); %定义画图范围，横轴，纵轴 title('抽样信号') %定义图的标题名字 6 指数函数 ()at f t Ae = 例1.5：用MA TLAB 实现0.5()3t f t e = 7 正弦函数 2()cos( )t f t A T π?=+ 例1.6：用MA TLAB 实现正弦函数f(t)=3cos(10πt+1) 8 虚指数信号 例 虚指数信号 调用格式是f=exp((j*w)*t) t=0:0.01:15;

实验室常用检验仪器操作规范

检验仪器操作规范 1. 仪器分类和作业规范 理化检验仪器（序号从4.1-4.30） 4.1原子吸收仪 4.2离子色谱仪 4.3 pH/电导率仪（或其他pH计及电导率仪） 4.4电光分析天平 4.5电子天平/电子分析天平 4.6浊度仪 4.7糖度计 4.8余氯测定仪 4.9分光光度计 4.10阿贝折射仪 4.11低速台式离心机 4.12定氮仪 4.13密度/比重/浓度计 4.14比色管 4.15电热恒温干燥箱 生化检验仪器（序号从4.31-4.40） 4.31自动立式压力蒸汽灭菌锅 4.32生物显微镜 4.33生化培养箱 4.34霉菌培养箱 其他（序号从4.61-4.70） 4.61激光粒子计数器

备注：以下仪器操作规范按仪器分类顺序编写 4.1原子吸收仪 4.1.1. 准备工作： a）接通仪器主机电源，再接通计算机及打印机电源。打开电脑，待 Win dows95 屏幕 左下角显示Start将箭头指向Start，并点击。联机正常后出示 Aawinlab，点击 打开。 b）将空白溶液置于自动进样器位置1上，标准溶液于2上，试样溶液 在其他编码位置上。 c）打开氩气钢瓶，并调节为300-450Kpa d）接通石墨炉冷却水系统电源。 4.1.2. 编制程序： 工作前将下列主要工作页的必要参数一一输入。分别为仪器页、校正页、石墨炉页。 4.1.3. 进行自动分析： 箭头指向Workspace并点击它，出现下图所示的对话框，然后点击 Calibrate 进行校正曲线，完毕点击An alyze Samples 进行试样的测 ^定。 备注：作结束后如需存储文件，点击File上的Save,显示Save Method As窗口，在Method Name^栏上打上方法名字，并点击0K 4.1.4. 关机： 工作结束后点击Automated Analysis Control 窗口的Flush Sampler ，冲洗 进样系统；退出软件，关掉主机就、计算机、打印机、冷却水系统和石墨炉电 源；关紧氩气钢瓶。 4.2离子色谱仪 4.2.1.仪器 离子色谱仪、移液管、滤纸（0.45 卩m、 4.2.2.试剂 阴离子淋洗液、阳离子淋洗液、甲烷磺酸 4.2.3.操作方法 a）洗液的配置: 1. 阴离子淋洗液的配置：

化学实验基本操作

化学实验基本操作 一、药品的取用 1、实验室所用的药品，很多是易燃、易爆、有腐蚀性或有毒的，为保证安全，在取用药品时，不能用手接触药品、不要把鼻孔凑到容器口闻药品，不得尝任何药品的味道。 2、取药品应按实验规定的用量取用，没有说明用量，一半应该按最少量取用：液体1-2ml，固体只需盖满试管底部。 3、实验剩余的药品既不能放回原瓶，也不要随意丢弃，更不要拿出实验室，要放入指定的容器内。 4、取用固体药品一般用广口瓶，块状的可用镊子夹取。用过的药匙或者镊子要立刻用干净得纸擦拭干净，以备下次使用。 5、往试管里装入固体粉末时，为避免药品沾在管口和管壁上，应做到一斜二送三直立，即先使试管倾斜，把药匙小心送到试管底部，然后使试管直立起来，让药品全部落到底部。 6、把块状的药品或者密度较大的金属颗粒放入玻璃容器时应做到一横二放三慢竖，即应该先把容器横放，把药品放入容器口以后，再把容器慢慢地竖立起来，使药品或者金属颗粒缓缓地滑到容器的底部，以免打破容器。 7、液体药品通常盛在细口瓶里。取用药液时，先拿下瓶塞，倒放在桌面上，然后拿起瓶子，瓶口要紧挨着试管口，使液体缓慢地倒入试管。倾倒时，标签要朝向手心，以防止药品腐蚀标签。 8、取用一定量的液体药品，也可以用量筒量出体积，量液时量筒必须放平，视线要跟量筒内液体凹液面的最低处保持水平，如果俯视读数时所量体积大于（大于或小于）实际体积，仰视读数时所量体积小于（大于或小于）实际体积。 9、滴管是用来吸取和滴加少量试剂的一种仪器，取液后滴管滴加时垂直、悬空放在烧杯正上方，试剂便滴入烧杯中。 10、实验中要特别注意保护眼睛，万一眼睛里进入药液，要立即用水冲洗（切不可用手揉眼睛），洗的时候要眨眼睛，必要时请医生治疗，提倡使用防护眼镜。 二、托盘天平的使用 1、托盘天平能准确到0.1g。 2、化学实验称量的药品，常是一些粉末状或是易潮解的，有腐蚀性的药品，称量时应在两个托盘上各放一张质量相同的纸，如果称量易潮解、有腐蚀性的药品，应放在玻璃器皿如烧杯、表面皿里进行称量。 三、连接仪器装置 把玻璃管插入橡皮塞和胶皮管时，先要把玻璃管口用水润湿，在容器口塞橡皮塞时，切不可把容器放在桌上再使劲塞进塞子，因为这样做容易压破

Matlab实验

MATLAB实验报告 学校：湖北文理学院 学院：物理与电子工程学院 专业：电子信息工程 学号： 2013128182 姓名：张冲 指导教师：宋立新

实验一 MATLAB环境的熟悉与基本运算 一、实验目的： 1．熟悉MATLAB开发环境 2．掌握矩阵、变量、表达式的各种基本运算 二、实验内容 1、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明， 学习使用指令eye（其它不会用的指令，依照此方法类推） 2、学习使用clc、clear，观察command window、command history和workspace 等窗口的变化结果。 3、初步程序的编写练习，新建M-file，保存（自己设定文件名，例如exerc1、 exerc2、exerc3……），学习使用MATLAB的基本运算符。 三、练习 1)help rand，然后随机生成一个2×6的数组，观察command window、 command history和workspace等窗口的变化结果。 2)学习使用clc、clear，了解其功能和作用。 3)用逻辑表达式求下列分段函数的值 4)求[100,999]之间能被21整除的数的个数。（提示：rem,sum的用法） 四、实验结果 1）

2）clc:清除命令窗口所有内容，数值不变；clear:初始化变量的值。3） 4）

实验二 MATLAB数值运算 一、实验目的 1、掌握矩阵的基本运算 2、掌握矩阵的数组运算 二、实验内容 1)输入C=1:2:20，则C（i）表示什么？其中i=1,2,3, (10) 2)输入A=[7 1 5；2 5 6；3 1 5]，B=[1 1 1; 2 2 2; 3 3 3]，在命令窗 口中执行下列表达式，掌握其含义： A(2, 3) A(:,2) A(3,:) A(:,1:2:3) A(:,3).*B(:,2) A(:,3)*B(2,:) A*B A.*B A^2 A.^2 B/A B./A 3)二维数组的创建和寻访，创建一个二维数组（4×8）A，查询数组A第2 行、第3列的元素，查询数组A第2行的所有元素，查询数组A第6列的所有 元素。 4)两种运算指令形式和实质内涵的比较。设有3个二维数组A 2×4，B 2×4 ，C 2×2 ， 写出所有由2个数组参与的合法的数组运算和矩阵指令。 5)学习使用表4列的常用函数（通过help方法） 6)学习使用表5数组操作函数。 7)生成一个3行3列的随机矩阵，并逆时针旋转90°，左右翻转，上下翻转。 8)已知a=[1 2 3]，b=[4 5 6]，求a.\b和a./ b 9)用reshape指令生成下列矩阵，并取出方框内的数组元素。 三、实验结果 1）C(i)表示C中的第i个的数值；

实验一MATLAB基本操作及运算.doc

实验一MATLAB基本操作及运算 一、实验目的 二、实验的设备及条件 三、实验内容 1、建立以下标量： 1） a=3 2） b=5+ 3 j，（j为虚数单位） 3） c=e j 2 / 3 2、建立以下向量： 2.71 38 1） Vb= 28 82 2） Vc=[4 3.8-3.8 -4 ] （向量中的数值从 4 到 -4，步长为 -0.2） 3、建立以下矩阵： 1） 3 L 3 Ma M O M 3 L 3 Ma 为一个 7×7的矩阵，其元素全为 3. 2） 1 11 L91 2 12 O92 Mb M M O M 10 20 L100

Mb 为一个 10× 10的矩阵 . 3） 1 14 5 Mc 2 5 17 3 23 8 4、使用题 1 中的变量计算下列等式的x,y,z的值： 1） x1 1 1 e( (a 15)/6) 2）x2 ( a 15 x1) 3）x3 ln( R [( b c)(b c)]sin( a / 3)) ，其中R表示复数实部。 5、求解函数值 y e ct2/(2.252)，其中 c 取值见题 1， t 的取值范围为题 2 中行 向量 Vc。 6、使用题 1 和题 3 中所产生的标量和矩阵计算等式 Mx a Mc (Mc) 1 (Mc )T 其中 * 为矩阵所对应行列式的值，参考det 。 7、函数的使用和矩阵的访问。 1）计算矩阵 Mb 每一列的和，结果应为行向量形式。 2）计算整个矩阵 Mb 的平均值。 3）用向量 [1 1 1] 替换 Mb 的最上一行的值 4）将矩阵 Mb 的第 2~5 行，第 3 到 9 列的元素所构成的矩阵赋值给矩阵SubMb。 5）删除矩阵 Mb 的第一行； 6）使用函数 rand 产生一个 1× 10的向量 r ，并将 r 中值小于 0.5 的元素设置为 0。 8、已知 CellA（1, 1）＝‘中国’， CellA（ 1， 2）＝‘北京’， CellA（ 2， 1）是一个 3 乘 3 的单位阵， CellA（ 2, 2）＝[1 2 3]，试用 MATLAB创建一个 2 ×2 的细胞数组 CellA。 9、已知结构数组student 中信息包含有姓名，学号，性别，年龄和班级，试用 MATLAB创建相应的结构数组 student。该数组包含有从自己学号开始连续 5 个同学的信息（如果学号在你后面的同学不足 5 个则往前排序），创建完成后查看自己的信息。

Matlab实验二

实验二DFS和DFT 实验任务 1、阅读并输入实验原理中介绍的例题程序，观察输出的图形曲线，理解每一条语句的含义。 2、已知一个信号序列的主值为x(n)=[0,1,2,3,2,1,0]，显示两个周期的信号序列波形，要求： （1）用DFS求信号的幅度频谱和相位频谱，用图形表示； （2）求IDFS的图形，并与原信号进行比较。 N = 7; xn = [0,1,2,3,2,1,0]; xn = [xn xn]; n = 0:2*N-1; k = 0:2*N-1; Xk = xn*exp(-j*2*pi/N).^(n'*k); x = (Xk*exp(j*2*pi/N).^(n'*k))/(2*2*N); subplot(2,2,1),stem(n,xn); title('x(n)');axis([-2,2*N,1.1*min(xn),1.1*max(xn)]); subplot(2,2,2),stem(n,abs(x)); title('IDFS|X(k)|');axis([-2,2*N,1.1*min(xn),1.1*max(xn)]); subplot(2,2,3),stem(k,abs(Xk)); title('|X(k)|');axis([-2,2*N,1.1*min(abs(Xk)),1.1*max(abs(Xk))]); subplot(2,2,4),stem(k,angle(Xk)); title('arg|X(k)|');axis([-2,2*N,1.1*min(angle(Xk)),1.1*max(angle(Xk)) ]);

3、已知有限长序列x(n)=[7,6,5,4,3,2]，求x(n)的DFT 和IDFT ，要求： （1） 画出DFT 对应的()X k 和()arg X k ????的图形。 （2）画出原信号与傅里叶逆变换IDFT[X(k)]图形进行比较 xn = [7 6 5 4 3 2]; N = length(xn); n = 0:N-1; k = 0:N-1; Xk = xn*exp(-j*2*pi/N).^(n'*k); x = (Xk*exp(j*2*pi/N).^(n'*k))/N; figure,subplot(2,2,1),stem(n,xn); title('x(n)'); subplot(2,2,2),stem(n,abs(x)); title('IDFT|X(k)|'); subplot(2,2,3),stem(k,abs(Xk)); title('|X(k)|'); subplot(2,2,4),stem(k,angle(Xk)); title('arg|X(k)|'); axis([0,N,1.1*min(angle(Xk)),1.1*max(angle(Xk))]) 4、一周期序列的主值x(n)=[7,6,5,4,3,2],求x(n)周期周期重复次数为3次时的DFS 。要求： （1）画出原主值和信号周期序列信号； （2）画出序列傅里叶变换对应的()X k %和()arg X k ????%的图形。 xn = [7 6 5 4 3 2];