matlab程序集(1)

《matlab程序集(1)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、1、编程求满足的最小m值。解：m=1;s=0;whiles<=10000s=s+2^m;m=m+1;endm=m-1;disp('m=');disp(m);运行结果：m=132、编写一个函数，计算下面函数的值，给出标量x的值，调用该函数后，返回y的值。function[y]=myfun1(x)选择一些数据测试你编写的函数。functiony=myfun1(x)ifx<=0y=sin(x);elseifx>0&x<=3y=x;elsey=6-x;endreturn;3、编写一个函数，给出一个向量，生成如下范德蒙矩阵。

2、function[v]=myvander(x)例如：>>v=myvander([2345])得v=生成一些数据测试你写的函数。function[v]=myvander(x)n=length(x);v(1:n)=1;fori=n+1:n*na=i/n;r=mod(i,n);b=floor(a);ifr==0r=n;b=b-1;endv(i)=x(r)^b;endv=reshape(v,n,n)';>>myvander([123])ans=111123149>>myvander([1468])ans=111114681

3、1636641642165124.求从1到100的数字的和。a=1;s=0;whilea>0&&a<=100;s=s+a;a=a+1;end;s第三个实验的实验内容5．Xb(n)的时域和幅频特性：n=0:50;x=zeros(1,51);x(1)=1;closeall;subplot(3,1,1);stem(n,x);title('单位冲击信号序列x(n)');k=-25:25;X=x*(exp(-j*pi/25)).^(n'*k);magX=abs(X);subplot(3,1,2);stem(n,magX);t

4、itle('单位冲击信号的幅度谱');angX=angle(X);subplot(3,1,3);stem(n,angX);title('单位冲击信号序列的相位谱');************************************************6．Hb（n）的时域和幅频特性：n=1:50;x=zeros(1,50);x(1)=1;x(2)=2.5;x(3)=2.5;x(4)=1;closeall;subplot(3,1,1);stem(x);title('系统单位脉冲响应信号序列');k=-24:

5、25;X=x*(exp(-j*pi/25)).^((n-1)'*k);magX=abs(X);subplot(3,1,2);stem(magX);title('系统频率响应的相位谱')*************************************************7.xC（n）的时域和幅频特性：Subplot(3,1,1);Stem(x);Title('矩形脉冲序列');K=-24:25;X=x*(exp(-j*pi/25)).^((n-1))*k;magX=abs(X);Subpolt(3,1,

6、2);stem(magX);Title('矩形脉冲序列博立叶变换的幅度谱');angX=angle(X);Subpolt(3,1,3);stem(angX);Title('矩形脉冲序列博立叶变换的相位谱');**************************************************************************************************8．a=[1-10.9];b=[1];n=0:30;x=(1/2).^n;h=filter(b,a,x);stem(

7、n,h,'fill'),gridonxlabel('n'),title('系统单位取样响应h(n)')****************************************************************************************************9.matlap程序；n=0:50;%定义序列的长度是50A=444.128;%设置信号有关的参数a=50*sqrt(2.0)*pi;T=0.001;%采样率w0=50*sqrt(2.0)*pi;%ω符号在MatLab

8、中不能输入，用w代替x=A*exp(-a*n*T).*sin(w0*n*T);%pi是MATLAB定义的π，信号乘可采用“.*”closeall%清除已经绘制的x(n)图形subplot(3,1,1);stem(x);%绘制x(n)的图形title('理想采样信号序列');%设置结果图形的标题*********************************