Matlab例题.doc

MATLAB原理及应用实验报告题目：MATLAB符号运算 专业：通信工程 姓名： 学号：实验程序1、建立和（1）求矩阵A和B的积，矩阵A左除B，以及矩阵A的2次方程序： A=1 2 -1;2 5 -1;4 10 3; B=2 -1 3;3 1 -5;4 -1 1; Y1=A*BY1 = 4 2 -8 15 4 -20 50 3 -35 Y2=ABY2 = 2.8000 -8.8000 31.6000 -0.6000 3.6000 -13.2000 -0.4000 -0.6000 2.2000 Y3=A2Y3 = 1 2 -6 8 19 -10 36 88 -5（2）求数组A和B的积，数组A左除B，及数组A的2次方程序： Z1=A.*BZ1 = 2 -2 -3 6 5 5 16 -10 3 Z2=A.BZ2 = 2.0000 -0.5000 -3.0000 1.5000 0.2000 5.0000 1.0000 -0.1000 0.3333 Z3=A.2Z3 = 1 4 1 4 25 1 16 100 92、建立一个字符串向量，然后对该向量做如下处理：(1) 取第15个字符组成的子字符串。程序： T=fdklsfncdADHCHSW %建立一个字符串T =fdklsfncdADHCHSW t=T(1:5) %取前五个字符组成一个新的字符串t =fdkls(2) 将字符串倒过来重新排列。程序： len=length(T) %先求出字符串的长度，然后对字符串取反len = 16 B=T(16:-1:1)B =WSHCHDAdcnfslkdf(3) 将字符串中的小写字母变成相应的大写字母，其余字符不变。程序： find(T=a&T T(1:1:9)=T(1:1:9)-32T =FDKLSFNCDADHCHSW(4) 统计字符串中小写字母的个数。程序： B=length(find(T=a&T3、使用simple和simplify两个指令分别化简，比较两个结果有什么不同程序：用simple化简 syms f x %定义f x两个变量 f=(1/x3+6/x2+12/x+8)(1/3) f = (1/x3+6/x2+12/x+8)(1/3) simple(f) %用simple函数化简 simplify: (2*x+1)3/x3)(1/3) radsimp: %函数对含根式表达式进行化简 (2*x+1)/x combine(trig): (1+6*x+12*x2+8*x3)/x3)(1/3) factor: (2*x+1)3/x3)(1/3) expand: (1/x3+6/x2+12/x+8)(1/3) combine: (1/x3+6/x2+12/x+8)(1/3) convert(exp): (1/x3+6/x2+12/x+8)(1/3) convert(sincos): (1/x3+6/x2+12/x+8)(1/3) convert(tan): (1/x3+6/x2+12/x+8)(1/3) collect(x): (1/x3+6/x2+12/x+8)(1/3) mwcos2sin: (1/x3+6/x2+12/x+8)(1/3) ans = (2*x+1)/x 用simplify函数化简 simplify(f) %用simplify函数化简 ans = (2*x+1)3/x3)(1/3) 4、求级数和1+x+x2+xk+的和。（使用symsum函数）程序： syms k symsum(1/(k2),1,inf) %变量k的范围时1到无穷 ans = 1/6*pi2 syms x symsum(xk,k,0,inf) %k的范围时0到无穷 ans = -1/(x-1) MATLAB原理及应用实验报告题目：MATLAB程序设计 专业：通信工程 姓名： 学号：实验程序及结果 1、计算的1100中所有素数的和，判断是否为素数是将100以内的每个数都被2整除，不能被整除的就是素数。程序：sum=0; %初始化sum=0flag=0; %初始化flag=0 for m=1:100; %建立1到100的向量 for k=2:ceil(m(1/2); %建立2到sqrt(m)的向量 if rem(m,k)=0; %判断m是否能整除k flag=1; %若都能整除就使flag=1 break; end end if flag=0; %判断是否符合是素数的条件 sum=sum+m; %累加 end flag=0;endfprintf(sum=%d,sum)结果为:sum = 10592、计算1+3+5.+100 的值，当和大于1000时终止计算。sum=0;for m=1:2:100; %建立1 3 5.100的向量 if sum=1000 %如果sum小于1000则可以继续加 sum=sum+m; %累加 else break; %若sum的结果不符合条件就跳出整个循环 endend结果为：sum = 10243、求1！+2！+3！+8！的值程序为: n=1;sum=0; for m=1:8; %循环7次使得得到各次阶乘 n=n*m; sum=sum+n; %累加 end结果为:sum = 462334、若一个数等于它的各个真因子之和，则称该数为完数，如6=1+2+3，所以6是完数。求1,500之间的全部完数。程序为:count=0; for m=1:500; %建立1到500的向量 sum=0; for k=1:ceil(m/2); %要建立1到m/2的向量 if rem(m,k)=0; %看是否整除 sum=sum+k; end end if sum=m; count=count+1; end end结果为:count = 4m = 1m = 6m = 28m = 4965、写程序，判断一年是否为闰年，符合下面两条件之一：A、能被4整除，不能被100整除B、能被4整除，不能被400整除程序为:count=0; for y=1988:2011; if(rem(y,4)=0&rem(y,100)=0)| (rem(y,4)=0&rem(y,400)=0); count=count+1;endend结果为: count =5MATLAB原理及应用实验报告题目：MATLAB图形操作 专业：通信工程 姓名： 学号：实验程序及结程序为1：t=0:0.1:4*pi; y=sin(t); y1=cos(t); plot(t,y,r,t,y1,-b) text(1.7*pi,-0.3,sint) text(1.6*pi,0.8,cost) 2：x=0:0.1:15;y1=sin(1*x);subplot(2,2,1);plot(x,y1);y2=cos(1*x);subplot(2,2,2);plot(x,y2);y3=(1*exp(x);subplot(2,2,3);plot(x,y3);y4=x.2;subplot(2,2,4);plot(x,y4);3 程序为： t=0:0.01:4*pi;y0=exp(-t/3);y1=-exp(-t/3);y=exp(-t/3).*sin(3*t);plot(t,y0,-b,t,y,r,t,y1,-b);title(y=exp(t/3).*sin(3*t)及它的包络线,fontsize,10,fontweigh,bold,fontangle,italic);text(1,0.5,学号0605074201,fontsize,10,fontwe