MATLAB-实验代码

1、实验 11、ones语句：Y = ones(n) %生成n×n 全1 阵Y = ones(m,n) %生成m×n 全1 阵Y = ones(m n) %生成m×n 全1 阵Y = ones(d1,d2,d3) %生成d1×d2×d3×全1 阵或数组Y = ones(d1 d2 d3) %生成d1×d2×d3×全1 阵或数组Y = ones(size(A) %生成与矩阵A 相同大小的全1 阵2、find语句：k = find(x) %按行检索X 中非零元素的点，若没有非零元素，将返回空矩阵。i,j = fi

2、nd(X) %检索X 中非零元素的行标i 和列标j。i,j,v = find(X) %检索X 中非零元素的行标i 和列标j 以及对应的元素值v。实验 31、编写一M函数，a和x作为函数参数输入，函数里面分别用if结构实现函数表示function output=function1(x,a)result=0;if x<=-a result=-1;elseif x>-a&x<a result=x/a; else x>=a result=1; end output=result;2、编写一M函数，迭代计算，给出可能的收敛值，其中x的初值作为函数的参数输入。functio

3、n output=function2(x)y=0;while 1 y=3/(x+2); if abs(y-x)<0.000001 break; else x=y; endendoutput=y;end3、编写一M函数，实现近似计算指数，其中x为函数参数输入，当n+1步与n步的结果误差小于0.00001时停止，分别用for和while 结构实现。用for结构实现：function output=function3(x)e=0;s=0;i=0;while 1 e=s+xi/factorial(i); if abs(e-s)<0.000001 break; else s=e; i=i+1

4、; endendoutput=e;用while 结构实现：function output=function4(x)s=1;e=0;for i=1: 1000000 e=s+xi/factorial(i); if abs(e-s)<0.00001 break; else s=e; endendoutput=e;实验 41 绘制函数的曲线，其中曲线为绿虚线，并进行标注x=-2:0.1:1;y=x.2;plot(x,y,'-g');hold onx=1:0.1:2;y=exp(-(x-1).2);plot(x,y,'-g')text(-1,1,'曲线y1

5、=x2'); text(1.5,0.5,'曲线y1=e(-(x-1)2)');2. 将用于绘制曲线的数据分别保存成MAT、二进制的文本文件中t=0:pi/10:2*pi;x=sin(t) sin(t);y=cos(t) cos(t);z=sin(t).*cos(t);save mydatafile x y zclearload mydatafiled=zt=0:pi/10:2*pi;x=sin(t);y=cos(t);z= sin(t).* cos(t);fid=fopen('text.txt','w')count=fwrite(fid,

6、z,'float32')closestatus=fclose(fid)3. 重启Matlab，从上述保存的文件中依次读取变量z的前10个数据fid=fopen('mydatafile.mat','r')data1=fread(fid,10)fid=fopen('text.txt','r');data2=fread(fid,10)实验 51. 将多项式A的系数向量形式1 3 6 3 1转换为完整形式，并求其根。同时在0-5内随机产生150组自变量，计算他们的对应取值；A=1 3 6 3 1;s,len=poly2st

7、r(A,'x') 求根：r=roots(A) 产生0-5内150组自变量：q=5*rand(1,150)计算对应取值：w=polyval(A,q)2. 对于上述150组数据，采用多项式进行拟合，并对分别采用最邻近、双线性和三次样条插值方法进行插值；x=1:4;y=x.4 + 3.* x.3 + 6.* x.2 + 3.* x + 1 ;xi=1:0.1:10;yi_nearest=interp1(x,y, xi,'nearset');yi_linear=interp1(x,y, xi);yi_spline=interp1(x,y, xi,'spline&

8、#39;);figure;hold on;subplot(1,3,1);plot(x,y,'ro',xi, yi_nearest,'b-'); title('最邻近插值');subplot(1,3,2);plot(x,y,'ro',xi, yi_linear,'b-');title('线性插值'); subplot(1,3,3);plot(x,y,'ro',xi, yi_spline,'b-');title('三次样条插值'); 3、 计算syms x

9、 yt=int(y*exp(x),y,2*y)int(t,1,2)实验 61、连续信号的产生与可视化，直流信号、正弦交流信号、单位阶跃信号、单位冲击信号、符号信号、斜坡信号、单位衰减信号、复指数信号等实现以及可视化直流信号：x1=-5:0.01:0;y1=1;plot(x1,y1); hold onx2=0:0.01:5;y2=1;plot(x2,y2)正弦交流信号：x=0:0.01:2*pi;y=sin(x);plot(x,y)单位阶跃信号：t=-4:0.01:4;f=(t>0);stairs(t,f);axis(-4,4,-1.1,1.1);单位冲击信号：t=-4:0.01:4;n=

10、length(t);f=zeros(1,n);f(1,(-t0+4)/0.01+1)=1; plot(t,f);axis(-4,4,-1.1,1.1);符号信号：t=-4:0.01:4;f=sign(t); plot(t,f);axis(-4,4,-1.1,1.1);斜坡信号：t=0:0.01:4;f=t;plot(t,f);axis(0,4,0,4);单位衰减信号：t=0:0.01:10;y=10*exp(-10*t).*sin(10*t);plot(t,y);复指数信号：t=-10:0.1:10;y=exp(10+j*t);plot(t,y);2、将信号与信号进行加、减、乘运算，并且将结果

11、可视化。加法运算：t=-10:0.1:10;y=exp(-3*t)+0.2*sin(4*pi*t);plot(t,y);减法运算：t=-10:0.1:10;y=exp(-3*t)-0.2*sin(4*pi*t);plot(t,y);乘法运算：t=-10:0.1:10;y=(exp(-3*t).*(0.2*sin(4*pi*t);plot(t,y);3已知信号，试通过反褶、移位、尺度变换由的波形得到的波形。定义符号函数f(t)=sin(t)/t：syms t;f=sym('sin(t)/t'); 对f进行移位：f1=subs(f,t,t+3); 对f1进行尺度变换： f2=subs(f1,t,2*t); 对f2进行反褶：f3=subs(f2,t,-t); 绘制函数波形： subplot(2,2,1);ezplot(f,-8,8);grid on;subplot(2,2,2);ezplot(f1,-8,8);grid on;subplot(2,2,3);ezplot(f2,-8,8);grid on;subplot(2,2,4);ezplot(f3,-8,8);grid on;4、分别产生两个方波信号，并且求这两个方波的卷积。产生两个方波：y1=ones(1,10),zeros(1,20);y2=ones(1,30),zeros(1,20);两个方波卷积：y=c