matlab基础实验

实验一MATLAB基础知识一、实验目的1．熟练掌握Matlab的启动与退出；2.熟悉Matalb的命令窗口，常用命令和帮助系统；3.熟悉Matalb的数据类型；4.熟悉Matlab的基本矩阵操作，运算符和字符串处理二、实验设备1．方正电脑2．MATLAB软件三、实验内容1.已知矩阵A=[34-11-910;65074-16;1-47-16-8;2-45-612-8;-36-78-11;8-49130]写出完成下列操作的命令:(1)将矩阵A的第2-5行中第1，3，5列元素赋给矩阵B；(2)删除矩阵A的第7号元素；(3)将矩阵A的每个元素加上30；(4)求矩阵A的大小和维数；(5)将矩阵A的右下角3*2矩阵构成矩阵C；(6)输出[-5,5]范围内的所有元素；程序：A=[34-11-910;65074-16;1-47-16-8;2-45-612-8;-36-78-11;8-49130]%(1)B=A(2:5,[135])%(2)A2=A;A2(7)=[]%(3)A3=A+30%(4)length(A)size(A)%(5)C=A(end-2:end,end-1:end)%(6)A6=A.*(A=-5&A=5)结果:A=34-11-91065074-161-47-16-82-45-612-8-36-78-118-49130B=6041762512-3-7-1A2=Columns1through73612-385Columns8through14-4-46-4-107Columns15through215-7917-1-6Columns22through2881-94612-1Columns29through35310-16-8-810A3=333429312140363530373414312637293622322635244222273623382931382639313330ans=6ans=66C=12-8-1130A6=34-11000500401-40-1002-45000-3000-110-401302.设x=-74o,y=27o,求22sin()tanxyxy的值。程序：x=(-74/180)*pi;y=(27/180)*pi;sin(x^2+y^2)/(sqrt(tan(abs(x+y)))+pi)结果:ans=0.22733.当a取-3.0，-2.9，-2.8，…，2.8，2.9，3.0时，求0.3sin(0.3)aea在各点的函数值。程序：a=[-3.0:0.1:3.0]m=exp(-0.3.*a).*sin(a+0.3)结果：a=Columns1through4-3.0000-2.9000-2.8000-2.7000Columns5through8-2.6000-2.5000-2.4000-2.3000Columns9through12-2.2000-2.1000-2.0000-1.9000Columns13through16-1.8000-1.7000-1.6000-1.5000Columns17through20-1.4000-1.3000-1.2000-1.1000Columns21through24-1.0000-0.9000-0.8000-0.7000Columns25through28-0.6000-0.5000-0.4000-0.3000Columns29through32-0.2000-0.100000.1000Columns33through360.20000.30000.40000.5000Columns37through400.60000.70000.80000.9000Columns41through441.00001.10001.20001.3000Columns45through481.40001.50001.60001.7000Columns49through521.80001.90002.00002.1000Columns53through562.20002.30002.40002.5000Columns57through602.60002.70002.80002.9000Column613.0000m=Columns1through4-1.0512-1.2305-1.3863-1.5184Columns5through8-1.6267-1.7116-1.7734-1.8129Columns9through12-1.8309-1.8285-1.8069-1.7675Columns13through16-1.7117-1.6411-1.5572-1.4617Columns17through20-1.3564-1.2428-1.1228-0.9978Columns21through24-0.8696-0.7397-0.6095-0.4804Columns25through28-0.3538-0.2308-0.11260.0000Columns29through320.10600.20470.29550.3779Columns33through360.45150.51600.57140.6174Columns37through400.65430.68210.70100.7115Columns41through440.71380.70850.69590.6768Columns45through480.65160.62100.58560.5460Columns49through520.50300.45720.40930.3597Columns53through560.30930.25860.20800.1582Columns57through600.10970.06280.0180-0.0245Column61-0.06414.已知3542126734457,2874879015930AB求下列表达式的值：(1)A*B和A.*B(2)A^3和A.^3(3)A/B和B\A(4)[A,B]和[A([1,3]);B^2]程序：A=[3542;34-457;879015];B=[1-267;2874;930];A*BA.*BA^3A.^3A/BA\B[A,B][A([1,3],:);B^2]结果;ans=12943241977-407-105240259112489ans=3-10813468-3605187832700ans=-28917240246-4368137883-2591012766917133325250438673ans=27157464839304-911253436585037290003375ans=-5.88455.3549-0.20286.3554-5.65964.3293-6.63256.20789.0241ans=0.1071-0.0667-8.71560.0172-0.0708-0.1894-0.12401.011751.6870ans=35421-26734-4572874879015930ans=3542879015600183-816842827261568255.设a=1+2i,b=-3+4i,6ice,求abcab的值。程序:a=1+2*i;b=3+4*i;c=exp((pi/6)*i);c+a*b/(a+b)结果:ans=1.6353+1.8462i实验二MATLAB程序设计一、实验目的1．学会编写简单函数m文件，调用及调试函数m文件；2．学会MATLAB的简单程序编写。二、实验设备1．方正电脑2．MATLAB软件三、实验内容1．从键盘输入一个4位整数，按如下规则加密后输出。加密规则：每位数字都加上7，然后用和除以10的余数取代该数字；然后将第1位数与第3位数交换，第2位数与第4位数交换。程序：a=input('请输入一个四位整数:');a1=fix(a/1000);%千a2=rem(fix(a/100),10);%百a3=rem(fix(a/10),10);%十a4=rem(a,10);%个%加密a1=a1+7;a2=a2+7;a3=a3+7;a4=a4+7;a1=rem(a1,10);a2=rem(a2,10);a3=rem(a3,10);a4=rem(a4,10);m=a3*1000+a4*100+a1*10+a2;ifa3==0disp(['输出：','0',num2str(m)])elsedisp(['输出：',num2str(m)])end结果：请输入一个四位整数:1234输出：01892.硅谷公司员工的工资计算方法如下：(1)工作时数超过120h者，超过部分加发15%；(2)工作时数低于60h者，扣发700元；(3)其余按每小时84元计发。试编程按输入的工号和该员工的工时数，计算应发工资。程序：clearall;clc;a=input('请输入工号:');b=input('请输入工时:');if60=b&b=120m=b*84;endifb60m=b*84-700;endifb120m=(b-120)*84*(1+0.15)+120*84;enddisp(['工号为',num2str(a),'的员工应发工资',num2str(m),'元'])结果：请输入工号:2453请输入工时:24工号为2453的员工应发工资1316元3．根据222211116123n，求的近似值。当n分别取100，1000，10000时，结果是多少？（要求：分别用循环结构和向量运算（使用sum函数）来实现）程序(1):clc;;clearall;n=input('输入n:');v=0;fori=1:nv=[v,1/(i^2)];endm=sum(v);disp(sqrt(6*m))%pi的近似值结果：输入n:10003.1406程序(2):clc;clearall;n=input('输入n:');v=0;fori=1:nv=v+1/(i^2);enddisp(sqrt(6*v))%pi的近似值结果：输入n:10003.14064．根据11113521yn，求(1)y3时的最大n值；(2)与(1)的n值对应的y值。考虑程序：clc;clearall;y=0;n=0;while(y=3)n=n+1;y=y+1/(2*n-1);endy=y-1/(2*n-1);n=n-1;disp(['y3时的最大的n=',num2str(n)]);disp(['y=',num2str(y)]);结果：y3时的最大的n=56y=2.99445．考虑一下迭代公式1nnaxbx，其中，ａ，ｂ为常数。(1)编写程序求迭代结果，迭代的终止条件为5110nnxx，迭代初值01.0x，迭代次数不超过500次。(2)如果迭代过程收敛于r，那么r的准确值是242bba，当(a,b)的值取(1,1)、(8,3)、(10,0.1)时，分别对迭代结果和准确值进行比较。程序：clc;clearall;a=input('请输入正数a:');b=input('请输入正数b:');n=0;%迭代次数x=1.0;%初值xl=1.0;while(n=499)x=a/(b+x);n=n+1;h=abs(x-xl);ifh=10^(-5)o=1;break;endxl=x;enddisp(['x=',num2str(x)]);disp(['n=',num2str(n)]);m=(-b+sqrt(b^2+4*a))/2结果：请输入正数a:123请输入正数b:134x=0.91171n=3m=0.9117实验三图形绘制一、实验目的1．加深认识二维、三维绘图函数plot,plot3,mesh,surf；2．要求熟练掌握以上绘图函数的用法、简单图形标注、简单颜色设定等。二、实验设备1．方正电脑2．MATLAB软件三、实验内容1．绘制下列曲线(1)33!xyx程序：x=-100: