实验三 MATLAB图像处理基本操作及摄像机标定(DLT) 实验三 Matlab图像处理基本操作及摄像机标定(DLT) (DLT)1、实验目的
通过应用Matlab的图像处理基本函数，学习图像处理中的一些基础操作和处理。
理解摄像机标定(DLT)方法的原理，并利用程序实现摄像机内参数和外参数的估计。。
2、实验内容:
1) 读取一幅图像并显示。
2) 检查内存(数组)中的图像。
3) 实现图像直方图均衡化。
4) 读取图像中像素点的坐标值。
5) 保存图像。
6) 检查新生成文件的信息。
7) 使用阈值操作将图像转换为二值图像。
8) 根据RGB图像创建一幅灰度图像。
9) 调节图像的对比度。
10) 在同一个窗口内显示两幅图像。
11) 掌握matlab命令及函数，获取标定块图像的特征点坐标。 12) 根据摄像机标定(DLT)方法原理，编写Matlab程序，估计摄像机内参数和12) 外参数。
3、实验要求:
1) 选取一幅图像，根据实验内容1)—10)给出结果。
2) 根据给定的标定块图像及实验内容11)，12)进行编程实验。 3) 书写实验报告
4、实验设备
1) 微机。
2) Matlab软件。
5、实验原理
DLT变换:
Abdal-Aziz和Karara于70年代初提出了直接线性变换像机定标的方法，他们从摄影测量学的角度深入的研究了像机图像和环境物体之间的关系，建立了像机成像几何的线性模型，这种线性模型参数的估计完全可以由线性方程的求解来实现。
直接线性变换是将像点和物点的成像几何关系在齐次坐标下写成透视投影矩阵的形式:
X，，wu，，,, Y,,w,,,svP 3，4,,,,Zw ,,1,,，，1，，
为图像坐标系下的点的齐次坐标，为世界坐标系下的空其中，，，，
u,v,1X,Y,Zwww
间点的欧氏坐标， P为3*4的透视投影矩阵，为未知尺度因子。
消去S，可以得到方程组:
pX，pY，pZ，p,puX,puY,puZ,pu,011w12w13w1431w32w33w34
pX，pY，pZ，p,puX,puY,puZ,pu,021w22w23w1431w32w33w34
当已知N个空间点和对应的图像上的点时，可以得到一个含有2*N个方程的方程组: AL,0
其中A为(2N*12)的矩阵， L为透视投影矩阵元素组成的向量:
T ,,p,p,p,p,p,p,p,p,p,p,p,p111213142122232431323334
||AL||像机定标的任务就是寻找合适的L，使得为最小，即 min || AL || L
给出约束: p,134
'T,1T L,,(CC)CB
L‘为L的前11个元素组成的向量， C为A前11列组成的矩阵， B为A前12列组成的向量。
6、程序代码
I=imread ('DSCN0831.JPG');
imshow(I);
whos
figure,
I1=rgb2gray(I);
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
4 -9.2702 -9.2702 -9.2702 -9.386
5 -9.379
6 -9.287
-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
-18.9315 -18.6146 -19.7033 -17.5651 -18.4268 -18.9672
-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 基本操作及运算 一、 实验目的 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中所产生的向量和矩阵计算以下等式，注意本题的操作
第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程序设计基本方法 覃照乘自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基本操作 题目： 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. 不能加热的仪器 仪器图形与名称主要用途使用方法及注意事项 用于收集和贮存少量气体上口为平面磨砂，内侧不磨砂，玻璃片 要涂凡士林油，以免漏气，如果在其中 进行燃烧反应且有固体生成时，应在底 部加少量水或细砂。 分装各种试剂，需要避光保存时用 棕色瓶。广口瓶盛放固体，细口瓶 盛放液体。 瓶口内侧磨砂，且与瓶塞一一对应，切 不可盖错。玻璃塞不可盛放强碱，滴瓶 内不可久置强氧化剂等。 制取某些气体的反应器固体＋液 体 固体为块状，气体溶解性小反应无强热 放出，旋转导气管活塞控制反应进行或 停止。 （二）计量仪器 仪器图形与名称主要用途使用方法及注意事项
实验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优化建模过程。 二、实验设备与工具 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矩阵及其运算 学号：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基本操作 实验报告 课程名称： 院系： 专业班级： 学号： 学生姓名： 指导教师： 开课时间：至学年第学期
一、学生撰写要求 按照实验课程培养方案的要求，每门实验课程中的每一个实验项目完成后，每位参加实验的学生均须在实验教师规定的时间内独立完成一份实验报告，不得抄袭，不得缺交。 学生撰写实验报告时应严格按照本实验报告规定的内容和要求填写。字迹工整，文字简练，数据齐全，图表规范，计算正确，分析充分、具体、定量。 二、教师评阅与装订要求 1.实验报告批改要深入细致，批改过程中要发现和纠正学生实验报告中的问题，给出评语和实验报告成绩，签名并注明批改日期。实验报告批改完成后，应采用适当的形式将学生实验报告中存在的问题及时反馈给学生。 2.实验报告成绩用百分制评定，并给出成绩评定的依据或评分标准（附于实验报告成绩登记表后）。对迟交实验报告的学生要酌情扣分，对缺交和抄袭实验报告的学生应及时批评教育，并对该次实验报告的分数以零分处理。对单独设课的实验课程，如学生抄袭或缺交实验报告达该课程全学期实验报告总次数三分之一以上，不得同意其参加本课程的考核。 3.各实验项目的实验报告成绩登记在实验报告成绩登记表中。本学期实验项目全部完成后，给定实验报告综合成绩。 4.实验报告综合成绩应按课程教学大纲规定比例（一般为10-15%）计入实验课总评成绩；实验总评成绩原则上应包括考勤、实验报告、考核（操作、理论）等多方面成绩； 5.实验教师每学期负责对拟存档的学生实验报告按课程、学生收齐并装订，按如下顺序装订成册：实验报告封面、实验报告成绩登记表、实验报告成绩评定依据、实验报告（按教学进度表规定的实验项目顺序排序）。装订时统一靠左侧按“两钉三等分”原则装订。
实验二 信号的表示及其基本运算 一、实验目的 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实验报告 学校：湖北文理学院 学院：物理与电子工程学院 专业：电子信息工程 学号： 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基本操作及运算 一、实验目的 二、实验的设备及条件 三、实验内容 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 个则往前排序），创建完成后查看自己的信息。
实验二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];
MATLAB实验二运算基础答案 以下是为大家整理的mATLAb实验二运算基础答案的相关范文，本文关键词为mATLAb,实验,运算,基础,答案,实验,mATLAb,运，您可以从右上方搜索框检索更多相关文章，如果您觉得有用，请继续关注我们并推荐给您的好友，您可以在综合文库中查看更多范文。 实验二、mATLAb运算基础 一、实验目的 掌握mATLAb各种表达式的书写规则及常用函数的使用。掌握mATLAb中字符串、元胞数组和结构的常用函数的使用。二、实验内
容及步骤 1、设有矩阵A和b，A=[12345;678910;1112131415;1617 181920;2122232425]，b=[3016;17-69;023-4;970;41311] 1）求它们的乘积c 2）将矩阵c的右下角3x2子矩阵赋给D >>A=[1:1:5;6:1:10;11:1:15;16:1:20;21:1:25];b=[3016;17-69;023-4;9 70;41311];...c=A*b,D=c(3:5,2:3)c=9315077258335237423520397588705 557753890717D=5203977055578907172、完成下列操作 1）求[100，999]之间能被61整除的数及其个数（提示：先 利用冒号表达式，再利用find和length函数。） >>A=[100:999];b=find(rem(A,61)==0),b_total=length(b) b=2384145206267328389450511572 633694755816877b_total=15 2）建立一个字符串向量，删除其中的大写字母（提示：利 用find函数和空矩阵。） >>a=['mATLAbisimportant'],b=abs(a);c=find(b A=12222222222s=2047 7、建立矩阵A并回答有关问题 A{1,1}='1';A{1,2}='brenden';A{2,1}=reshape(1:9,3,3);A{2,2}={12,34,2;54,2 1,3;4,23,67}; 1）使用访问的方式如何将brenden改写成bRenDen？2）分别执
实验一 MATLAB 基本操作 一、实验目的 1. 学习和掌握MA TLAB 的基本操作方法 2. 掌握命令窗口的使用 3. 熟悉MATLAB 的数据表示、基本运算 二、实验内容和要求 1. 实验内容 1) 练习MATLAB7.0或以上版本 2) 练习矩阵运算与数组运算 2. 实验要求 1) 每位学生独立完成，交实验报告 2) 禁止玩游戏！ 三、实验主要软件平台 装有MATLAB7.0或以上的PC 机一台 四、实验方法、步骤及结果测试 1. 实验方法：上机练习。 2. 实验步骤： 1) 开启PC ，进入MA TLAB 。 2) 使用帮助命令，查找sqrt 函数的使用方法 答: help sqrt 3) 矩阵、数组运算 a) 已知 ??????????=987654321A ，???? ??????=963852741B ，求)2()(A B B A -?+ 答： A=[1, 2, 3; 4, 5, 6; 7, 8, 9]; B=[1, 4, 7; 2, 5, 8; 3, 6, 9]; (A+B)*(2*B-A) b) 已知?? ????-=33.1x ，??????=π24y ，求T xy ，y x T c) 已知??????????=987654321A ，???? ??????=300020001B ，求A/B, A\B. d) 已知???? ??????=987654321A ，求：(1) A 中第三列前两个元素；(2) A 中所有第二行元素；(3) A 中四个角上的元素；(4) 交换A 的第1、3列。(5) 交换A 的第1、2行。(6) 删除A 的第3列。
e) 已知[]321=x ，[]654=y ，求：y x *.，y x /.，y x \.，y x .^， 2.^x ，x .^2。 f) 给出x=1,2,…,7时，x x sin 的值。 3）常用的数学函数 a ）随机产生一个3x3的矩阵A ，求：(1) A 每一行的最大、最小值，以及最大、最小值所在的列；(2) A 每一列的最大、最小值，以及最大、最小值所在的行；(3) 整个矩阵的最大、最小值；(4) 每行元素之和；(5) 每列元素之和；(6) 每行元素之积；(7) 每列元素之积。 b) 随机产生两个10个元素的向量x ，y 。（1） 求x 的平均值、标准方差。（2） 求x ，y 的相关系数。（3）对x 排序，并记录排序后元素在原向量中的位置。 4) 字符串操作函数 建立一个字符串向量（如‘ABc123d4e56Fg9’），然后对该向量做如下处理： (1) 取第1～5个字符组成的子字符串。 (2) 将字符串倒过来重新排列。 (3) 将字符串中的小写字母变成相应的大写字母，其余字符不变。 (4) 统计字符串中小写字母的个数。