当前位置:首页 > Matlab课程复习大纲(1)
6.1、写出生成下图所示波形的MATLAB脚本文件M文件。图中虚线为正弦波,要求它的负半波被置零,且在3/2处被削顶。
t=linspace(0,3*pi,500); y=sin(t);
z1=((t
w=(t>pi/3&t<2*pi/3)+(t>7*pi/3&t<8*pi/3); w_n=~w; z2=w*sin(pi/3)+w_n.*z1; plot(t,y,'r-.',t,z2,'k-') xlabel('x'); ylabel('z=f(t)');
title('逐段解析函数'); legend('y=sin(t)','z=f(t)');
6.2运用关系运算求近似极限,画出如图所示完整的插值函数图形。
t=-2*pi:pi/10:2*pi; y=sin(t)./t;
subplot(1,2,1),plot(t,y),axis([-7,7,-0.5,1.2]), xlabel('t'),ylabel('y'),title('残缺图形') tt=t+(t==0)*eps; yy=sin(tt)./tt;
subplot(1,2,2),plot(tt,yy),axis([-7,7,-0.5,1.2]) xlabel('t'),ylabel('yy'),title('正确图形')
7、编程实现y=sin(x),x从0到2,x=0.02,求y的最大值、最小值、均值和标准差 x=0:0.02:2; y=sin(x);
ymax=max(y); ymin=min(y); ymean=mean(y); ystd=std(y)
8、编写一个函数,使其能够产生如下的分段函数:
?0.5x?f(x)??1.5?0.25x?0.5?
function f=fun(x) if x<=2
f=0.5*x; elseif x<=6
f=1.5-0.25*x; else
x?22?x?6 x?6 f=0.5; end
9、对下列程序的每句进行注释,说明每句程序完成的功能。 x=0:pi/50:2*pi;
y1=sin(x); y2=cos(x); figure(1)
plot(x, y1,'k-.', x, y2, 'b-'); grid on
legend('sin(x)','cos(x)'); text(pi,0,'sin(x)'); gtext('cos(x)');
title('sin(x) and cos(x)'); xlabel('x');
ylabel('sin(x) and cos(x)');
10、 设a1=[2, 4, 6, 8]为系统分母系数向量,b1=[3, 6, 9]为系统分子系数向量,求此系统的频率响应并画出频率特性图。 a1=[2,4,6,8]; b1=[3,6,9]; freqs(b1,a1 11、
计算二重积分S??0???2?(ysinx?xcosy)dxdy
S=dblquad('y*sin(x)+x*cos(y)',pi,2*pi,0,pi)
Q=dbquad(inline('y*sin(x)+x*cos(y)'),pi,2*pi,0,pi)
12、 生成一个4 x 4的随机矩阵,对其进行三角分解和正交分解,并求其特征根和特征向量。
A=rand(4,4);
[L,U]=lu(A); [Q,R]=qr(A); [V,D]=eig(A)
13、求下列联立方程的解:
3x?4y?7z?12w?45x?7y?4z?2w??3x?8z?5w?9?6x?5y?2z?10w??8i26314、请分别写出用for和while循环语句计算K??2?1?2?2?....?2的程序,以
i?063
及写出一种避免循环的程序。 K=0; K=0; i=0:33; for i=0:33 K=K+2.^i; end K i=0; while i<34 K=K+2^i; i=i+1; end K y=2.^i; K=sum(y)
15、阅读下述MATLAB脚本程序 A=[3 5 7 9;0 2 4 1;6 4 4 8]; B=ones(3,3); b=A([1 2],2:3); B(4,4)=5; A([1,2],:)=[]; C=[A;B(:,1:3)']; D=reshape(C,2,8);
写出b、A、B、C、D的结果。
b=[5 7;2 4] A=[6 4 4 8]
B=[1 1 1 0;1 1 1 0;1 1 1 0;0 0 0 5] C=[6 4 4 8; 1 1 1 1;1 1 1 1;1 1 1 0] D=[6 1 4 1 4 1 8 0;1 1 1 1 1 1 0 0]
16、读程序
A=[1,2,7,11;22,14,10,32] s=0; for a=A s=s+a; end s
S=21;78
17、写出下列程序运行后变量a,b,c,d的值。 a=eye(4) b=fliplr(a) c=diag(a,1)
d=reshape(a,2,8)
a=[1 0 0 0;0 1 0 0;0 0 1 0;0 0 01] b=[0 0 0 1;0 0 1 0;0 1 0 0;1 0 0 0] c=[0 0 0]
d=[1 0 0 0 0 1 0 0;0 0 1 0 0 0 0 1]
共分享92篇相关文档