2010专第一次上机作业： MatLab练习题-基本运算
MATLAB实验：运算基础,并且附有答案
实验二、MATLAB运算基础一、实验目的掌握MATLAB各种表达式的书写规则及常用函数的使用。
掌握MATLAB中字符串、元胞数组和结构的常用函数的使用。
二、实验内容及步骤1、设有矩阵A和B，A=[1 2 3 4 5;6 7 8 9 10;11 12 13 14 15;16 1718 19 20;21 22 23 24 25]，B=[3 0 16;17 -6 9;0 23 -4;9 7 0;4 1311]1）求它们的乘积C >>C=A*B2）将矩阵C的右下角3x2子矩阵赋给D >>I=[3 4 5];J=[2 3];D=C(I,J)也可以用>>D=C([3 4 5],[2 3])D =520 397705 557890 7172、完成下列操作1）求[100，999]之间能被61整除的数及其个数（提示：先利用冒号表达式，再利用find和length函数。
）>> a=100:999;find(rem(a,61)==0)ans =23 84 145 206 267 328 389 450 511572 633 694 755 816 877>> b=a(ans)b =122 183 244 305 366 427 488 549 610 671 732 793 854 915 976>> length(b)ans =152）建立一个字符串向量，删除其中的大写字母（提示：利用find函数和空矩阵。
）a=’I am maying’;a( find(a>’A’&a<’Z’))=[]3、已知A=[23 10 -78 0；41 -45 65 5；32 5 0 32；6 -54 92 14]，取出其前3行构成矩阵B，其前两列构成矩阵C，其左下角3x2子矩阵构成矩阵D，B与C的乘积构成矩阵E，分别求E<D、E&D、E|D、~E|~D。
Matlab上机练习题及答案
Matlab 上机练习题及答案---------------------------------------------------------------------1、矩阵Y=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡3472123100451150425，给出元素1的全下标和单下标，并用函数练习全下标和单下标的转换，求出元素100的存储位置。
取出子矩阵⎥⎦⎤⎢⎣⎡21301，并求该矩阵的维数。
解：命令为：Y=[5,2,4;0,15,1;45,100,23;21,47,3]Y(2,3)Y(10)sub2ind([43],2,3)[i,j]=ind2sub([43],10)find(Y==100)sub2ind([43],3,2)B=Y(2:2:4,3:-2:1)或B=Y([24],[31])[m n]=size(Y)---------------------------------------------------------------------2、已知矩阵A=[10-1；241；-205]，B=[0-10；213；112]求2A+B 、A 2-3B 、A*B 、B*A 、A.*B ，A/B 、A\B 解：命令为：A=[10-1；241；-205]B=[0-10；213；112]E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B---------------------------------------------------------------------3、利用函数产生3*4阶单位矩阵和全部元素都为8的4*4阶矩阵，并计算两者的乘积。
解：命令为：A=eye(3,4)B=8*ones(4)C=A*B---------------------------------------------------------------------4、创建矩阵a=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------7023021.5003.120498601,取出其前两列构成的矩阵b ，取出前两行构成矩阵c ，转置矩阵b构成矩阵d ，计算a*b 、c<d ,c&d,c|d ,~c|~d 解：命令为：a=[-1,0,-6,8;-9,4,0,12.3;0,0,5.1,-2;0,-23,0,-7]b=a(:,[12])c=a([12],:)d=b’e=a*b f=c<d g=c&d h=c|d i=~c|~d---------------------------------------------------------------------5、求!201∑=n n 解：命令文件为sum=0;s=1;for n=1:20s=n*s;sum=sum+s;end sum---------------------------------------------------------------------6、求a aa aaa aa a S n ++++=得值，其中a 是一个数字，由键盘输入，表达式中位数最多项a 的个数，也由键盘输入。
MATLAB上机练习题(一)
MATLAB 上机练习题（一）
1. 计算下式
7
2009/42008/41949222?+-⨯-⨯=
2. 试计算d
3461.1,30,5b d a b c
a b c +=-+===其中
要求只输出d 的值，其他变量的值不输出。
**************************************************************
3. 计算下式：
181(2)sin()cos()6,3
y
x x y π-+-==（，其中 0.1104
2log 805(3)2009
6(4)1,5;
2e
e x x x -⨯+-=其中
4. 完成以上题目后清除工作空间中的所有变量。
**************************************************************
5. 求下列方程的根：
（1）325x x -=
（2）32340x x x -+=
6. 根据课本任意找出4个矩阵在MA TLAB 中输入和输出。
附1：增加的一些常用函数pi：常数
round(x)：四舍五入至最近整数
fix(x) ：无论正负，舍去小数至最近整数min(x): 向量x的元素的最小值
max(x): 向量x的元素的最大值
diff(x): 向量x的相邻元素的差
length(x): 向量x的元素个数
sin(x)：正弦函数
cos(x)：余弦函数
tan(x)：正切函数
asin(x)：反正弦函数
acos(x)：反余弦函数
atan(x)：反正切函数。
MATLAB运算基础答案
实验 01 讲评、参考答案讲评未交实验报告的同学名单数学： 6 人（11、12 级）信科： 12-04, 12-22, 13-47批改情况：问题 1：不仔细，式子中出错。
问题 2：提交的过程不完整。
问题 3：使用语句尾分号 (;)不当，提交的过程中不该显示的结果显示。
问题 4：截屏窗口没有调整大小。
附参考答案：王平《MATLAB软件》课内实验实验 01 MATLAB 运算基础（第 2 章 MATLAB 数据及其运算）一、实验目的1.熟悉启动和退出 MATLAB 的方法。
2.熟悉 MATLAB 命令窗口的组成。
3.掌握建立矩阵的方法。
4.掌握 MATLAB 各种表达式的书写规则以及常用函数的使用。
二、实验内容1.数学表达式计算先求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量。
计算三角函数z12sin 850（注意：度要转换成弧度，e2如何给出）1 e2示例：点击Command Window 窗口右上角的，将命令窗口提出来成悬浮窗口，适当调整窗口大小。
命令窗口中的执行过程：计算自然对数z 21ln( x21 2i 1 x 2) ，其中 x520.45（提示： clc 命令擦除命令窗口， clear则清除工作空间中的所有变量，使用时注意区别，慎用clear 命令。
应用点乘方）命令窗口中的执行过程：求数学表达式的一组值z 3 e0.3ae 0.3 a sin(a0.3) ln0.3 a, a3.0, 2.9,L , 2.9, 3.022提示：利用冒号表达式生成 a 向量，求各点的函数值时用点乘运算。
命令窗口中的执行过程：求分段函数的一组值t 20t1z4 t 211t2，其中 t=0:: t 22t 1 2t3提示：用逻辑表达式求分段函数值。
命令窗口中的执行过程：对工作空间的操作接着显示 MATLAB 当前工作空间的使用情况并保存全部变量提示：用到命令who, whos, save, clear, load ，请参考教材相关内容。
Matlab上机题目及一段软件
《Matlab 上机题目》1、基本数学运算与初等函数(1) p=1,计算 10*log(2510*2⎪⎭⎫⎝⎛-p ) (2)a=2,b=5;计算225131++=b b a y(3)计算y=)/tan()(21b a b a e ab -(4)计算)/tan(*(a)tan(b)ctan a)sin(b)cos(b a ac *acsin(a/b)(5)计算 (1+2i)*(1-7i)*ipi e**2-2、基本逻辑、关系运算a=2,b=5,c=6 (1) c>b||a<b (2) c>b&&a<b (3) c>b&&a>b (4) c>b||a>b(5) c==a&&c~=b (6) c==a||c~=bd=[ 1 2] e= [1 0] (1) d|e (2) d&e (3) d>e (4) d<e (5) d==e (6) d~=e (7) d>=e (8) d<=e3、数组（1） 生成一个-5到5，步长为0.5的随机行数组a,判断其中数值大于等于0，小于1，且不等于0.5的元素个数，并找出这些元素，保存到变量b 中。
（2） 用randn 函数生成一个5行5列的随机数组c ，用rand 生成一个6行，4列的随机数组d ，并用d 中位于第2行到第4行，第3列到第4列的元素替换c 中第3行到第5行，第3列到第4列的元素。
（3） 将c 中所有奇数行，偶数列的元素取出来组成一个新的数组c1（4） 用eye 生成一个5行5列的单位数组e1，用randn 函数生成一个5行5列的随机数组c ，计算f=2e1+12e c ,并将f 中对角线上的元素取出来，保存到变量g 中。
4、脚本程序编写(1) 文件编写程序找出1到2345中能被17整除的数的个数以及这些数的和。
{帮助：求余数的函数用c=mod(a,b) ,其中a 为被除数，b 为除数, c 为余数; 对数组m 求和用he=sum(m), 返回值he 为数组的和}以下是一个用MATLAB 语言编写的声音片段软件>> fs = 44100; % sample rate dt = 1/fs; T16 = 0.125; t16 = [0:dt:T16]; [temp k] = size(t16);t4 = linspace(0,4*T16,4*k); t8 = linspace(0,2*T16,2*k); [temp i] = size(t4); [temp j] = size(t8);% Modification functionsmod4=(t4.^4).*exp(-27*(t4.^0.5)); mod4=mod4/max(mod4);mod8=(t8.^4).*exp(-50*(t8.^0.75)); mod8=mod8/max(mod8);mod16=(t16.^4).*exp(-125*(t16.^1)); mod16=mod16/max(mod16);f0 = 2*146.8; % reference frequency ScaleTable = [2/3 3/4 5/6 15/16 ...1 9/8 5/4 4/3 3/2 5/3 9/5 15/8 ... 2 9/4 5/2 8/3 3 10/3 15/4 4 ... 1/2 9/16 5/8]; % 1/4 notesdo0f = mod4.*cos(2*pi*ScaleTable(21)*f0*t4); re0f = mod4.*cos(2*pi*ScaleTable(22)*f0*t4); mi0f = mod4.*cos(2*pi*ScaleTable(23)*f0*t4);fa0f = mod4.*cos(2*pi*ScaleTable(1)*f0*t4); so0f = mod4.*cos(2*pi*ScaleTable(2)*f0*t4);ti0f = mod4.*cos(2*pi*ScaleTable(4)*f0*t4); do1f = mod4.*cos(2*pi*ScaleTable(5)*f0*t4); re1f = mod4.*cos(2*pi*ScaleTable(6)*f0*t4); mi1f = mod4.*cos(2*pi*ScaleTable(7)*f0*t4); fa1f = mod4.*cos(2*pi*ScaleTable(8)*f0*t4); so1f = mod4.*cos(2*pi*ScaleTable(9)*f0*t4); la1f = mod4.*cos(2*pi*ScaleTable(10)*f0*t4); tb1f = mod4.*cos(2*pi*ScaleTable(11)*f0*t4); ti1f = mod4.*cos(2*pi*ScaleTable(12)*f0*t4); do2f = mod4.*cos(2*pi*ScaleTable(13)*f0*t4); re2f = mod4.*cos(2*pi*ScaleTable(14)*f0*t4); mi2f = mod4.*cos(2*pi*ScaleTable(15)*f0*t4); fa2f = mod4.*cos(2*pi*ScaleTable(16)*f0*t4); so2f = mod4.*cos(2*pi*ScaleTable(17)*f0*t4); la2f = mod4.*cos(2*pi*ScaleTable(18)*f0*t4); ti2f = mod4.*cos(2*pi*ScaleTable(19)*f0*t4); do3f = mod4.*cos(2*pi*ScaleTable(20)*f0*t4); blkf = zeros(1,i);% 1/8 notesfa0e = mod8.*cos(2*pi*ScaleTable(1)*f0*t8); so0e = mod8.*cos(2*pi*ScaleTable(2)*f0*t8); la0e = mod8.*cos(2*pi*ScaleTable(3)*f0*t8);ti0e = mod8.*cos(2*pi*ScaleTable(4)*f0*t8); do1e = mod8.*cos(2*pi*ScaleTable(5)*f0*t8); re1e = mod8.*cos(2*pi*ScaleTable(6)*f0*t8); mi1e = mod8.*cos(2*pi*ScaleTable(7)*f0*t8); fa1e = mod8.*cos(2*pi*ScaleTable(8)*f0*t8); so1e = mod8.*cos(2*pi*ScaleTable(9)*f0*t8); la1e = mod8.*cos(2*pi*ScaleTable(10)*f0*t8); tb1e = mod8.*cos(2*pi*ScaleTable(11)*f0*t8); ti1e = mod8.*cos(2*pi*ScaleTable(12)*f0*t8); do2e = mod8.*cos(2*pi*ScaleTable(13)*f0*t8); re2e = mod8.*cos(2*pi*ScaleTable(14)*f0*t8); mi2e = mod8.*cos(2*pi*ScaleTable(15)*f0*t8); fa2e = mod8.*cos(2*pi*ScaleTable(16)*f0*t8); so2e = mod8.*cos(2*pi*ScaleTable(17)*f0*t8); la2e = mod8.*cos(2*pi*ScaleTable(18)*f0*t8); ti2e = mod8.*cos(2*pi*ScaleTable(19)*f0*t8); do3e = mod8.*cos(2*pi*ScaleTable(20)*f0*t8); blke = zeros(1,j);% 1/16 notesfa0s = mod16.*cos(2*pi*ScaleTable(1)*f0*t16); so0s = mod16.*cos(2*pi*ScaleTable(2)*f0*t16);ti0s = mod16.*cos(2*pi*ScaleTable(4)*f0*t16);do1s = mod16.*cos(2*pi*ScaleTable(5)*f0*t16);re1s = mod16.*cos(2*pi*ScaleTable(6)*f0*t16);mi1s = mod16.*cos(2*pi*ScaleTable(7)*f0*t16);fa1s = mod16.*cos(2*pi*ScaleTable(8)*f0*t16);so1s = mod16.*cos(2*pi*ScaleTable(9)*f0*t16);la1s = mod16.*cos(2*pi*ScaleTable(10)*f0*t16);tb1s = mod16.*cos(2*pi*ScaleTable(11)*f0*t16);ti1s = mod16.*cos(2*pi*ScaleTable(12)*f0*t16);do2s = mod16.*cos(2*pi*ScaleTable(13)*f0*t16);re2s = mod16.*cos(2*pi*ScaleTable(14)*f0*t16);mi2s = mod16.*cos(2*pi*ScaleTable(15)*f0*t16);fa2s = mod16.*cos(2*pi*ScaleTable(16)*f0*t16);so2s = mod16.*cos(2*pi*ScaleTable(17)*f0*t16);la2s = mod16.*cos(2*pi*ScaleTable(18)*f0*t16);ti2s = mod16.*cos(2*pi*ScaleTable(19)*f0*t16);do3s = mod16.*cos(2*pi*ScaleTable(20)*f0*t16);blks = zeros(1,k);% Blank Blockblkblock = [blkf blkf blkf blkf blkf blkf blkf blkf...blkf blkf blkf blkf blkf blkf blkf blkf];% Base Melodycello = [do1f do1f so0f so0f la0f la0f mi0f mi0f...fa0f fa0f do0f do0f fa0f fa0f so0f so0f];% So-FUCKING-Long Melodyviolin = [mi2f mi2f re2f re2f do2f do2f ti1f ti1f...la1f la1f so1f so1f la1f la1f ti1f ti1f ...%do2f do2f ti1f ti1f la1f la1f so1f so1f...fa1f fa1f mi1f mi1f fa1f fa1f re1f re1f ...%do1f mi1f so1f fa1f mi1f do1f mi1f re1f...do1f la0f do1f so1f fa1f la1f so1f fa1f...%mi1f do1f re1f ti1f do2f mi2f so2f so1f...la1f fa1f so1f mi1f do1f do2f blkf blke ti1e ...%do2e ti1e do2e do1e ti0e so1e re1e mi1e...do1e do2e ti1e la1e ti1e mi2e so2e la2e...fa2e mi2e re2e fa2e mi2e re2e do2e ti1e...la1e so1e fa1e mi1e re1e fa1e mi1e re1e... %%do1e re1e mi1e fa1e so1e re1e so1e fa1e...mi1e la1e so1e fa1e so1e fa1e mi1e re1e...do1e la0e la1e ti1e do2e ti1e la1e so1e...fa1e mi1e re1e la1e so1e la1e so1e fa1e...%mi1f mi2e blke re2f re2f blkf do1f mi2f mi2f...la2f la2f so2f so2f la2f la2f ti2f ti2f...%do3e blke do2e blke ti1f ti1f blkf la1f do2f do2f...do2f do2f do2f do2f do2f fa2f re2f so2f...%so2e mi2s fa2s so2e mi2s fa2s so2s so1s la1s ti1s ...do2s re2s mi2s fa2s mi2e do2s re2s...mi2e mi1s fa1s so1s la1s so1s fa1s so1s mi1s fa1s so1s...fa1e la1s so1s fa1e mi1s re1s mi1s re1s do1s re1s mi1s fa1s so1s la1s...fa2e la1s so1s la1e ti1s do2s so1s la1s ti1s do2s re2s mi2s fa2s so2s...%mi2e do2s re2s mi2e re2s do2s re2s ti1s do2s re2s mi2s re2s do2s ti1s...do2e la1s ti1s do2e do1s re1s mi1s fa1s mi1s re1s mi1s do2s ti1s do2s...la1e do2s ti1s la1e so1s fa1s so1s fa1s mi1s fa1s so1s la1s ti1s do2s...la2e do2s ti1s do2e ti1s la1s ti1s do2s re2s do2s ti1s do1s la1s ti1s...%%do2e blke blkf ti1e blke blkf la1e blke blkf do2e blke blkf...do1e blke blkf do1e blke blkf do1e blke blkf do1e blke blkf...%blkf so1e blke blkf so1e blke blkf mi1e blke blkf so1e blke...blkf fa1e blke blkf mi1e blke blkf fa1e blke blkf re2e blke...%mi2e mi1e fa1e mi1e re1e re2e mi2e re2e do2e mi1e do1e do2e ti1e so0e fa0e so0e...la0e la1e so1e la1e so1e so0e fa0e so0e do1e la1e so1e la1e ti1e ti0e la0e ti0e...%do1e do2e re2e do2e ti1e ti0e do1e ti0e la0e la1e so1e la1e ti1e ti0e mi1e re1e...do1e do2e re2e fa2e mi2e mi1e so1e mi2e do2e fa2e mi2e fa2e re2e so1e fa1e so1e...%mi1e so1e so1e so1e so1e so1e so1e so1e mi1e mi1e mi1e mi1e mi1e mi1e so1e so1e...fa1e fa1e fa1e do2e do2e do2e do2e do2e do2e do2e la1e la1e so1e so1e re2e ti1e...%%so1e mi2e mi2e mi2e re2e re2e re2e re2e do2e do2e do2e do2e so2e so2e so2e so2e...la2e la2e la2e la2e so2e so2e so2e so2e la2e la2e la2e la2e ti2e ti1e ti1e ti1e...%do2e do1s re1s mi1e do1e ti0e ti1s do2s re2e ti1e la1e la0s ti0s do1e la0e ti0e so1s fa1s mi1e re1e...do1e mi1s re1s do1e fa1e mi1e do1s re1s mi1e so1e fa1e la1s so1s fa1e mi1e re1e so1s fa1s mi1e re1e...%mi1e do2s ti1s do2e mi1e so1e so1s la1s ti1e so1e mi1e do2s re2s mi2e do2e mi2e mi2s re2s do2e ti1e...la1e la1s so1s la1e ti1e do2e mi2s re2s do2e mi2e fa2e do2s ti1s la1e la1e so1e re1e so1e so1e...%so1f so1f so1f so1f do1f do1f do1f so1f...fa1f fa1f so1f so1f fa1f do1f do1f do1e ti0e...%do1f do2f ti1f ti1f la1f la1f so1f so1f...do1f do1e re1e mi1f mi1f do2f do2f ti1f ti1f...%%do2f];% celloc1 = [cello cello cello cello cello...cello cello cello cello cello...cello cello cello cello cello...cello cello cello cello cello...cello cello cello blkf];% violin1v1 = [blkblock violin blkblock blkblock];% violin2v2 = [blkblock blkblock violin blkblock];% violin3v3 = [blkblock blkblock blkblock violin];% Get dirtys = c1+v1+v2+v3;s = s/max(s);sound(s,fs); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%>> fs = 50000; % sample ratedt = 1/fs;T16 = 0.125;t16 = [0:dt:T16];[temp k] = size(t16);t4 = linspace(0,4*T16,4*k);t8 = linspace(0,2*T16,2*k);[temp i] = size(t4);[temp j] = size(t8);% Modification functions%mod4=(t4.^4).*exp(-27*(t4.^0.5));mod4=(t4.^5).*exp(-35*(t4.^0.5));mod4=mod4/max(mod4);mod8=(t8.^4).*exp(-50*(t8.^0.75));mod8=mod8/max(mod8);mod16=(t16.^4).*exp(-125*(t16.^1));mod16=mod16/max(mod16);f0 = 2*146.8; % reference frequencyScaleTable = [2/3 3/4 5/6 15/16 ...1 9/8 5/4 4/3 3/2 5/3 9/5 15/8 ...2 9/4 5/2 8/3 3 10/3 15/4 4 ...1/2 9/16 5/8];m=7;ratiom=[0,1/2,1/4,1/6,1/8,1/12,1/18];%基频% 1/4 notesdo0f = mod4.*cos(2*pi*ScaleTable(21)*f0*t4);re0f = mod4.*cos(2*pi*ScaleTable(22)*f0*t4);mi0f = mod4.*cos(2*pi*ScaleTable(23)*f0*t4);fa0f = mod4.*cos(2*pi*ScaleTable(1)*f0*t4);la0f = mod4.*cos(2*pi*ScaleTable(3)*f0*t4);ti0f = mod4.*cos(2*pi*ScaleTable(4)*f0*t4); do1f = mod4.*cos(2*pi*ScaleTable(5)*f0*t4); re1f = mod4.*cos(2*pi*ScaleTable(6)*f0*t4); mi1f = mod4.*cos(2*pi*ScaleTable(7)*f0*t4); fa1f = mod4.*cos(2*pi*ScaleTable(8)*f0*t4); so1f = mod4.*cos(2*pi*ScaleTable(9)*f0*t4); la1f = mod4.*cos(2*pi*ScaleTable(10)*f0*t4); tb1f = mod4.*cos(2*pi*ScaleTable(11)*f0*t4); ti1f = mod4.*cos(2*pi*ScaleTable(12)*f0*t4); do2f = mod4.*cos(2*pi*ScaleTable(13)*f0*t4); re2f = mod4.*cos(2*pi*ScaleTable(14)*f0*t4); mi2f = mod4.*cos(2*pi*ScaleTable(15)*f0*t4); fa2f = mod4.*cos(2*pi*ScaleTable(16)*f0*t4); so2f = mod4.*cos(2*pi*ScaleTable(17)*f0*t4); la2f = mod4.*cos(2*pi*ScaleTable(18)*f0*t4); ti2f = mod4.*cos(2*pi*ScaleTable(19)*f0*t4); do3f = mod4.*cos(2*pi*ScaleTable(20)*f0*t4); blkf = zeros(1,i);% 1/8 notesfa0e = mod8.*cos(2*pi*ScaleTable(1)*f0*t8); so0e = mod8.*cos(2*pi*ScaleTable(2)*f0*t8); la0e = mod8.*cos(2*pi*ScaleTable(3)*f0*t8); ti0e = mod8.*cos(2*pi*ScaleTable(4)*f0*t8); do1e = mod8.*cos(2*pi*ScaleTable(5)*f0*t8); re1e = mod8.*cos(2*pi*ScaleTable(6)*f0*t8); mi1e = mod8.*cos(2*pi*ScaleTable(7)*f0*t8); fa1e = mod8.*cos(2*pi*ScaleTable(8)*f0*t8); so1e = mod8.*cos(2*pi*ScaleTable(9)*f0*t8); la1e = mod8.*cos(2*pi*ScaleTable(10)*f0*t8); tb1e = mod8.*cos(2*pi*ScaleTable(11)*f0*t8); ti1e = mod8.*cos(2*pi*ScaleTable(12)*f0*t8); do2e = mod8.*cos(2*pi*ScaleTable(13)*f0*t8); re2e = mod8.*cos(2*pi*ScaleTable(14)*f0*t8); mi2e = mod8.*cos(2*pi*ScaleTable(15)*f0*t8); fa2e = mod8.*cos(2*pi*ScaleTable(16)*f0*t8); so2e = mod8.*cos(2*pi*ScaleTable(17)*f0*t8); la2e = mod8.*cos(2*pi*ScaleTable(18)*f0*t8); ti2e = mod8.*cos(2*pi*ScaleTable(19)*f0*t8); do3e = mod8.*cos(2*pi*ScaleTable(20)*f0*t8); blke = zeros(1,j);% 1/16 notesfa0s = mod16.*cos(2*pi*ScaleTable(1)*f0*t16);la0s = mod16.*cos(2*pi*ScaleTable(3)*f0*t16);ti0s = mod16.*cos(2*pi*ScaleTable(4)*f0*t16);do1s = mod16.*cos(2*pi*ScaleTable(5)*f0*t16);re1s = mod16.*cos(2*pi*ScaleTable(6)*f0*t16);mi1s = mod16.*cos(2*pi*ScaleTable(7)*f0*t16);fa1s = mod16.*cos(2*pi*ScaleTable(8)*f0*t16);so1s = mod16.*cos(2*pi*ScaleTable(9)*f0*t16);la1s = mod16.*cos(2*pi*ScaleTable(10)*f0*t16);tb1s = mod16.*cos(2*pi*ScaleTable(11)*f0*t16);ti1s = mod16.*cos(2*pi*ScaleTable(12)*f0*t16);do2s = mod16.*cos(2*pi*ScaleTable(13)*f0*t16);re2s = mod16.*cos(2*pi*ScaleTable(14)*f0*t16);mi2s = mod16.*cos(2*pi*ScaleTable(15)*f0*t16);fa2s = mod16.*cos(2*pi*ScaleTable(16)*f0*t16);so2s = mod16.*cos(2*pi*ScaleTable(17)*f0*t16);la2s = mod16.*cos(2*pi*ScaleTable(18)*f0*t16);ti2s = mod16.*cos(2*pi*ScaleTable(19)*f0*t16);do3s = mod16.*cos(2*pi*ScaleTable(20)*f0*t16);blks = zeros(1,k); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%倍频，可任意。
matlab练习题
matlab练习题Matlab是一种高级的技术计算软件，被广泛应用于科学与工程领域。
通过编写Matlab代码，我们可以解决各种数学问题、进行数据分析和可视化、进行仿真等。
本文将提供一些针对初学者的Matlab练习题，帮助读者熟悉Matlab的基本语法和常用操作。
练习一：矩阵运算在Matlab中，矩阵运算是一项重要的操作。
请编写Matlab代码完成以下任务：1. 创建两个3×3的矩阵A和B，元素的值自行选择。
2. 计算矩阵A和B的和，并将结果存储在变量C中。
3. 计算矩阵A和B的乘积，并将结果存储在变量D中。
4. 计算矩阵A的转置，并将结果存储在变量E中。
提示：使用Matlab的矩阵运算符和内置函数可以简化代码编写过程。
练习二：函数编写与图形绘制Matlab是一个功能强大的数学软件，可以用来编写自定义的函数并进行图形绘制。
1. 编写一个Matlab函数，该函数可以计算任意输入整数n的阶乘，并将结果输出。
2. 调用该函数，计算并输出5的阶乘。
3. 使用Matlab绘制函数y=2x+1的图像，其中x的取值范围为-10到10。
提示：使用Matlab的plot函数可以实现图形的绘制，并可使用grid函数添加刻度线。
练习三：数据分析Matlab提供了丰富的数据分析工具和函数，可以用来处理和分析实验数据。
1. 编写一个Matlab脚本，实现对给定实验数据进行统计分析的功能。
2. 实验数据保存在文本文件data.txt中，每行包含一个数据点。
读取数据并将其存储在一个向量中。
3. 计算数据的均值、方差和标准差，并将结果输出。
4. 使用Matlab的直方图函数绘制数据的频率分布直方图。
提示：使用Matlab的readmatrix函数可以读取文本文件中的数据。
调用mean、var和std函数可以计算均值、方差和标准差。
利用histogram函数可以绘制直方图。
练习四：数值解与绘图Matlab具有强大的数值计算和绘图功能，可以用来求解方程和绘制函数图像。
Matlab上机练习一答案
Matlab上机练习⼀答案Matlab 上机练习⼀班级学号姓名按要求完成题⽬，并写下指令和运⾏结果。
（不需要画图）1、求下列联⽴⽅程的解=+-+-=-+=++-=--+41025695842475412743w z y x w z x w z y x w z y x >> a=[3 4 -7 -12;5 -7 4 2;1 0 8 -5;-6 5 -2 10];>> b=[4;4;9;4];>> c=a\bc =5.22264.45701.47181.59942、设 ??++=)1(sin 35.0cos 2x x x y ，把x=0-2π间分为101点，画出以x 为横坐标，y 为纵坐标的曲线。
>> x=linspace(0,2*pi,101);>> y=cos(x)*(0.5+(1+x.^2)\3*sin(x));>> plot(x,y,'r')3、产⽣8×6阶的正态分布随机数矩阵R1, 求其各列的平均值、和、乘积。
并求该矩阵全体数的平均值。
（mean var ）a=randn(8,6)mean(a)var(a)k=mean(a)k1=mean(k)i=ones(8,6)i1=i*k1i2=a-i1i3=i2.*i2g=mean(i3)g2=mean(g)或者u=reshape(a,1,48);p1=mean(u)p2=var(u)4、绘制)(222y x e x z +-=在定义域x=[-2,2],y=[-2,2]内的曲⾯。
（利⽤meshgrid ）x=-2:2;y=x;[X,Y]= meshgrid(x,y);Z=X^2*exp(-(X^2+Y^2));mesh(X,Y,Z)5、求代数⽅程3x 5+4x 4+7x 3+2x 2+9x+12=0的所有根。
（利⽤roots 函数）p=[3 4 7 2 9 12];roots(p)6、把1开五次⽅，并求其全部五个根。
matlab基础练习题及答案
第1章MATLAB基础1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？答：在MATLAB操作桌面上有五个窗口。
在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口。
在独立窗口的Desktop菜单中选择Dock...项就可以将独立的窗口重新放置到桌面上。
1.5 如何启动M文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。
在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。
1.6 存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。
1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。
1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？答：当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。
在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。
1.9 在MATLAB中有几种获得帮助的途径？答：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器。
MATLAB数据及其运算_习题答案
MATLAB数据及其运算_习题答案第2章 MATLAB数据及其运算习题2⼀、选择题1．下列可作为MATLAB合法变量名的是（）。
DA．合计 B．123 C．@h D．xyz_2a 2．下列数值数据表⽰中错误的是（）。
CA．+10 B． C．2e D．2i3．使⽤语句t=0:7⽣成的是（）个元素的向量。
A A．8 B．7 C．6 D．54．执⾏语句A=[1,2,3;4,5,6]后，A(3)的值是（）。
B A．1 B．2 C．3 D．4 5．已知a为3×3矩阵，则a(:,end)是指（）。
D A．所有元素 B．第⼀⾏元素C．第三⾏元素 D．第三列元素6．已知a为3×3矩阵，则运⾏a (1)=[]后（）。
A A．a变成⾏向量 B．a变为2⾏2列C．a变为3⾏2列 D．a变为2⾏3列7．在命令⾏窗⼝输⼊下列命令后，x的值是（）。
B>> clear>> x=i*jA．不确定 B．-1 C．1 D．i*j 8．fix(354/100)+mod(354,10)*10的值是（）。
DA．34 B．354 C．453 D．439．下列语句中错误的是（）。
BA．x==y==3 B．x=y=3C．x=y==3 D．y=3,x=y10．find(1:2:20>15)的结果是（）。
CA．19 20 B．17 19C．9 10 D．8 911．输⼊字符串时，要⽤（）将字符括起来。
CA．[ ] B．{ } C．' ' D．" " 12．已知s='显⽰"hello"'，则s的元素个数是（）。
A A．9 B．11 C．7 D．1813．eval('sqrt(4)+2')的值是（）。
BA．sqrt(4)+2 B．4 C．2 D．2， 214．有3×4的结构矩阵student，每个结构有name（姓名）、scores（分数）两个成员，其中scores是以1×5矩阵表⽰的5门课的成绩，那么要删除第4个学⽣的第2门课成绩，应采⽤的正确命令是（）。
matlab上机考试题及答案
matlab上机考试题及答案1. 题目：编写一个MATLAB函数，计算并返回一个向量中所有元素的平方和。
答案：函数定义如下：```matlabfunction sumOfSquares = calculateSumOfSquares(vector)sumOfSquares = sum(vector.^2);end```2. 题目：使用MATLAB的内置函数，找出一个矩阵中的最大元素及其位置。
答案：可以使用`max`函数来找出矩阵中的最大元素，同时使用`find`函数来获取其位置。
示例代码如下：```matlabA = [1, 2, 3; 4, 5, 6; 7, 8, 9];[maxValue, linearIndex] = max(A(:));[row, col] = ind2sub(size(A), linearIndex);```3. 题目：给定一个向量，使用MATLAB编写代码，实现向量元素的逆序排列。
答案：可以使用`flip`函数来实现向量的逆序排列。
示例代码如下：```matlabvector = [1, 2, 3, 4, 5];reversedVector = flip(vector);```4. 题目：编写一个MATLAB脚本，计算并绘制一个正弦波的图像。
答案：可以使用`sin`函数生成正弦波数据，并使用`plot`函数绘制图像。
示例代码如下：```matlabx = linspace(0, 2*pi, 100);y = sin(x);plot(x, y);xlabel('x');ylabel('sin(x)');title('Sine Wave');```5. 题目：给定一个3x3的矩阵，使用MATLAB编写代码，计算其行列式。
答案：可以使用`det`函数来计算矩阵的行列式。
示例代码如下：```matlabmatrix = [1, 2, 3; 4, 5, 6; 7, 8, 9];determinant = det(matrix);```结束语：以上是MATLAB上机考试的题目及答案，希望能够帮助大家更好地掌握MATLAB的编程技巧和函数使用。
matlab试题及答案开卷
matlab试题及答案开卷1. MATLAB基础操作在MATLAB中，如何创建一个名为“myMatrix”的3x3矩阵，其元素从1到9按顺序排列？答案：在MATLAB中，可以使用以下命令创建名为“myMatrix”的3x3矩阵：```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```或者使用内置函数`magic`：```matlabmyMatrix = magic(3);```2. 矩阵运算给定两个矩阵A和B，其中A = [2 3; 4 5]，B = [6 7; 8 9]，计算A+B的结果。
答案：在MATLAB中，可以使用加号`+`来计算两个矩阵的和：```matlabA = [2 3; 4 5];B = [6 7; 8 9];C = A + B;```计算结果C为：```matlabC = [8 10; 12 14];```3. 向量操作创建一个从0到1的等差数列，步长为0.1。
答案：在MATLAB中，可以使用`linspace`函数或`:`操作符来创建等差数列：```matlabv = 0:0.1:1;```或者使用`linspace`函数：```matlabv = linspace(0, 1, 11);```两种方法都将得到一个包含11个元素的向量，从0开始，到1结束，步长为0.1。
4. 条件语句编写一个MATLAB脚本，判断一个给定的数n是否为素数。
答案：可以通过以下MATLAB脚本来判断一个数n是否为素数：```matlabfunction isPrime = isPrimeNumber(n)if n <= 1isPrime = false;elsefor i = 2:sqrt(n)if mod(n, i) == 0isPrime = false;return;endendisPrime = true;endend```调用此函数时，传入一个数值n，函数将返回一个布尔值，表示n 是否为素数。
matlab程序设计上机练习(2010研)C
1MATLAB 程序设计上机练习C1、按顺序进行如下的操作（1）用Matlab 命令产生一个4阶魔术方阵A ，并计算'A 与1-A ；（2）求的A 特征值；（3）计算A 的各列的总和与平均值；（4）若[1 2 3 6]'b =，求方程组A x b =的解。
2、 绘制下列曲线的图形：2cos ,[0,2]y x x π=∈。
在命令窗口输入：ezplot(cos(x^2),[0,2*pi]);123456-1-0.50.51xcos(x 2)3、假如有一组实测数据：[0.1 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.2];[42 43.5 45 45.5 45 47.5 49 53 50 55]; xy==求其回归直线，画回归直线图形并计算最小误差平方和。
1）回归直线求解：>> A=[0.1 1;0.11 1;0.12 1;0.13 1;0.14 1;0.15 1;0.16 1;0.17 1;0.18 1;0.2 1] >> b=[42;43.5;45;45.5;45;47.5;49;53;50;55]>> a=A\ba =124.675329.3474>> syms x>> y=124.6753*x+29.3474;2）回归直线2>> ezplot('y',[0.1,0.2]) 0.10.110.120.130.140.150.160.170.180.190.2yy3）残差平方和（最小误差平方和）;求残差的M文件：function S=fangcha(x,y)b=zeros(1,10);for i=1:103b(i)=124.6753*x(i)+29.3474;endbc=b-yS=0;for i=1:10S=S+c(i)^2;end命令窗口;>> a=[0.1 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.2]; >> b=[42 43.5 45 45.5 45 47.5 49 53 50 55];>> S=fangcha(a,b)计算结果;b =Columns 1 through 94541.8149 43.0617 44.3084 45.5552 46.8019 48.0487 49.2954 50.5422 51.7890 Column 10 54.2825 c =Columns 1 through 9-0.1851 -0.4383 -0.6916 0.0552 1.8019 0.5487 0.2954 -2.4578 1.7890 Column 10 -0.7175 S = 14.0990 4、求下列各结果：（1）用Matlab 展开：242(1)(1)x x x +-+；（2）用Matlab 求极限：lim (1)nn xL n→∞=-。
(完整版)matlab上机练习题答案
6 9 3 2 4 11.计算a 与b 的数组乘积2 7 5 4 6 8>> a=[6 9 3;2 7 5];>> b=[2 4 1;4 6 8];>> a.*bans =12 36 38 42 404 9 2 372.对于AX B，如果A 7 6 4 ，B 26，求解X o3 5 7 287 16 185计算多项式除法（3乂3+13*+6乂+8）/0+4）>> d=deconv([3 13 6 8],[1 4])6求欠定方程组298的最小范数解5>> a=[2 4 7 4;9 3 5 6];>> b=[8 5]';>> x=pinv(a)*b-0.2151>> A=[4 9 2;7 6 4;3 5 7]; >> B=[37 26 28]'；>> X=A\BX =-0.51184.04271.33180.44590.79490.27077用符号函数法求解方程at2+b*t+c=0 >> r=solve('a*t A2+b*t+c=0','t')1 2 5 8 7 43.a ，b ，观察a与b之间的3 64 3 6 2六种关系运算的结果[1/2/a*(-b+(bA2-4*a*c)A(1/2))] [1/2/a*(-b-(bA2-4*a*c)A(1/2))]>> a=[1 2 3;4 5 6]; >> b=[8 - 4;3 6 2]; >> a>bans =0 11 0 >> a>=bans =0 11 0 >> a<bans =1 00 1 >> a<=bans =1 00 1 >> a==bans =0 00 0 >> a~=bans =1 11 1 4计算多项式乘法111111(X2+2X+2)(X2+5X+4)8求矩阵A a11 a12的行列式值、逆和特征根a 21 a22>> syms a11 a12 a21 a22;>> A=[a11,a12;a21,a22]>> AD=det(A) % 行列式>> AI=inv(A) % 逆>> AE=eig(A) % 特征值A =[a11, a12][a21, a22]AD =a11*a22-a12*a21AI =[-a22/(-a11*a22+a12*a21), a12/(-a11*a22+a12*a21)][a21/(-a11*a22+a12*a21), -a11/(-a11*a22+a12*a21)]AE =>> c=conv([1 2 2],[1 5 4]) [1/2*a11+1/2*a22+1/2*(a11A2-2*a11*a22+a22A2+4*a12*a21)A(1/2)][1/2*a11+1/2*a22-1/2*(a11A2-2*a11*a22+a22A2+4*a12*a21)A(1/2)] 9 因式分解：X45X35X25X 6>> syms X;>> f=xA4-5*xA3+5*xA2+5*x-6;>> factor(f)ans =(X-1)*(X-2)*(X-3)*(X+1)x ，用符号微分求df/dx。
matlab程序设计上机练习(2010研)B
MATLAB 程序设计上机练习题（B 卷）
1．按顺序进行如下的操作：
（1）用Matlab 命令产生一个4阶魔术方阵A ，并计算'A 与1A -；
（2）求的A 特征值；
（3）计算A 的各列的总和与平均值；
（4）若[1 2 3 4]'b =，求方程组Ax b =的解。
2．在同一个图形上画出两个函数的图形：
22(1)sin(),[0,];(2)sin (),[0,]y x x y x x ππ=∈=∈
3．求下列各结果：
（1）用Matlabq 计算二阶导数：2cos2y x =；
（2）用Matlab 求极限：10(1)lim x
x x e L x →+-=。
4．用Matlab 命令求方程323210x x x --+=的根。
5．求微分方程''2'5cos()x y y y e x -+=的通解。
6．求微分方程组2530t dx x y e dt dy x y dt ⎧++=⎪⎪⎨⎪--=⎪⎩在初始条件001,0t t x y ====下的特解，并画出图
形。
7．某人做一种材料的伸缩实验,x 为温度(℃),y 为长度(mm),实验数据见下表 x y
18 79.6
20 81
25 82.3
30 84
35 86.8
40 89
求x与y的表达式.要求:1.画出散点图;2.写出x与y的关系式并画出回归直线方程.
8. 编程题目：水仙花数是指一个n 位数( n≥3 )，它的每个位上的数字的n 次幂之和等于它本身。
（例如：1^3 + 5^3 + 3^3 = 153）,试找出三位、四位、五位、六位、七位水仙花数。
matlab基础考试题及答案
matlab基础考试题及答案1. MATLAB中，如何创建一个名为"myMatrix"的3x3矩阵，其元素为1到9的连续整数？答案：可以使用以下命令创建矩阵：```matlabmyMatrix = [1 2 3; 4 5 6; 7 8 9];```或者使用`magic`函数：```matlabmyMatrix = magic(3);```2. 在MATLAB中，如何计算一个向量的元素平方？答案：可以使用`.^`运算符来计算向量元素的平方：```matlabv = [1 2 3 4];v_squared = v.^2;```3. MATLAB中，如何将一个向量逆序？答案：可以使用`flip`函数或者`end`索引来逆序一个向量：```matlabv = [1 2 3 4];v_reversed = flip(v);```或者```matlabv_reversed = v(end:-1:1);```4. 在MATLAB中，如何计算矩阵的行列式？答案：可以使用`det`函数来计算矩阵的行列式：```matlabA = [1 2; 3 4];detA = det(A);```5. MATLAB中，如何使用循环结构来计算1到100的自然数之和？答案：可以使用`for`循环来计算总和：```matlabsum = 0;for i = 1:100sum = sum + i;end```或者使用`sum`函数直接计算：```matlabsum = sum(1:100);```6. 在MATLAB中，如何创建一个复数？答案：可以通过指定实部和虚部来创建一个复数：```matlabcomplexNumber = 3 + 4i;```7. MATLAB中，如何绘制一个正弦波？答案：可以使用`plot`函数和`sin`函数来绘制正弦波：```matlabx = linspace(0, 2*pi, 100);y = sin(x);plot(x, y);```8. 如何在MATLAB中计算两个矩阵的点积？答案：可以使用`dot`函数或`*`运算符来计算两个矩阵的点积：```matlabA = [1 2; 3 4];B = [5 6; 7 8];dotProduct = dot(A, B);```或者```matlabdotProduct = A * B;```9. MATLAB中，如何使用条件语句来检查一个数是否为偶数？答案：可以使用`if`语句和模运算符来检查一个数是否为偶数：```matlabnum = 10;if mod(num, 2) == 0disp('The number is even.');elsedisp('The number is odd.');end```10. 在MATLAB中，如何将一个字符串转换为小写？答案：可以使用`lower`函数来将字符串转换为小写：```matlabstr = 'Hello World';lowerStr = lower(str);```结束语：以上是MATLAB基础考试题及答案，希望这些题目能帮助你更好地理解和掌握MATLAB的基础知识。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
一、书上命令练习：
二、MatLab 练习1
注意: 1,请按要求作题,写出程序并把结果附上,并保存文件.
2,仔细填写 “学生信息”.
3,把此word 文档的文件名更名为你的 “学号 姓名(*)”-*表示第几次作业.在下课前5分钟提交.
1, (5*2+1.3-0.8)*10^2/25
答案：
2，y =23log )10sin(4)43.0(2+⋅+-e
π 答案：
3，z = 10sin(π/3) ...
cos(π/3)
答案：
4，x = sin(223/3)， y = x^2，z = y*10 ；
求x+2y-5z
答案：
5，i e z i z i z ⋅==+=π23215,2,72,计算3
221z z z z z += 答案：
6, 建立起始值=3,增量值=5.5,终止值=44的一维数组x
(填写程序语句)
答案：
7, 建立等差一维数组x ：首项为0,末项为π,项数为15 (填写程序语句) 答案：
8, 计算行列式的值1
876243
23=A
答案：
9, 矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=187624323A ,矩阵⎥⎥⎥⎦
⎤⎢⎢⎢⎣⎡=333222111B ;分别求出A x B 及A 与B 中对应元素之间的乘积的程序语句。
答案
10,写出计算 Sin(36o )的程序语句.
答案：
三、MatLab 练习2-矩阵及数组运算
1.输入矩阵。
2.输入一个与A同阶的随机矩阵B。
要求元素为整数。
3.输入一个数组b=（1 3 5 7 9 11 ）
4.输入与A同阶的单位矩阵E
5.输入与B同阶的零矩阵Z
6.构建一个大矩阵
7.将 46-70这25个整数填入一个五行五列的矩阵数表G中，使其各行、各
列以及主对角线和次对角线的和等相等。
8.分别求d1=|A|、d2=|B|、d3=|D|、d4=|G|
9.对以上四个矩阵A、B、D、G中的可逆矩阵求逆。
10.(1) 计算d5=A+B
(2) 计算d6=A*B
11.取A的1，3，5行做一个子阵a1.
12.取A的2，4，6行，4，5，6列上的交叉元素做一个子阵a2.
13.由a2产生列向量a3.
14.取a3绝对值大于4的元素构成列向量a4,
15.改变A的元素A（4，5）为0 值。
16.改变A的第6列为b’
17.建立三角函数表阵S。
第一列是自变量X的取值，从0-2 共取20个值。
第二列是对应y=cos x的值。
18.分别取a2的1、2行作成两个向量A1、A2
19.作向量的点积C1=A1·A2
20.作向量的叉积C2=A1×A2。