MATLAB上机练习西电版1
matlab上机练习答案
实验一MA TLAB的基本命令与基本函数１已知矩阵a =11 12 13 1421 22 23 2431 32 33 3441 42 43 44求(1) A(:,1) (2) A(2,:)(3) A(:,2:3) (4) A(2:3,2:3)(5) A(:,1:2:3) (6) A(2:3)(7) A(:) (8) A(:,:)(9) ones(2,2) (10) eye(2)(11) [A,[ones(2,2);eye(2)]](12) diag(A) (13) diag(A,1)(14) diag(A,-1) (15) diag(A,2)2(1)输入如下矩阵Ａ０π/3Ａ＝π/6 π/2(2) 求矩阵B1,B1中每一元素为对应矩阵Ａ中每一元素的正弦函数(3) 求矩阵B2, B2中每一元素为对应矩阵Ａ中每一元素的余弦函数(4) 求B12+B22(5) 求矩阵Ａ的特征值与特征矢量：称特征矢量为Ｍ，而特征值矩阵为Ｌ(6) 求Msin(L)M-13已知水的黏度随温度的变化公式为μ=μ0/(1+at+bt2)其中μ0=1.785×10-3，a=0.03368,b=0.000221,求水在0，20，40，80℃时的黏度。
程序如下：miu0=1.785e-3;a=0.03368;b=0.000221;t=0:20:80miu=miu0./ (1+a*t+b*t.^2)（2）一个长管，其内表面半径为a，温度为Ta ；外表面半径为Tb；则其径向和切向应力可分别表示为：⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛+----=⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛-⎪⎭⎫ ⎝⎛⎪⎪⎭⎫ ⎝⎛----=r b a b r b a b a a b v T T E r b a b r b a b a a b v T T E b a t b a r ln ln 11)/ln()1(2)(ln ln 1)/ln()1(2)(2222222222ασασ式中r 为管子的径向坐标，E 为管子材料的弹性模量，ɑ为热膨胀系数。
西安电子科技大学数字信号处理上机作业
数字信号处理MATLAB上机作业M 2.21.题目The square wave and the sawtooth wave are two periodic sequences as sketched in figure ing the function stem. The input data specified by the user are: desired length L of the sequence, peak value A, and the period N. For the square wave sequence an additional user-specified parameter is the duty cycle, which is the percent of the period for which the signal is positive. Using this program generate the first 100 samples of each of the above sequences with a sampling rate of 20 kHz ,a peak value of 7, a period of 13 ,and a duty cycle of 60% for the square wave.2.程序% 用户定义各项参数参数A = input('The peak value =');L = input('Length of sequence =');N = input('The period of sequence =');FT = input('The desired sampling frequency =');DC = input('The square wave duty cycle = ');% 产生所需要的信号t = 0:L-1;T = 1/FT;x = A*sawtooth(2*pi*t/N);y = A*square(2*pi*(t/N),DC);% Plotsubplot(2,1,1)stem(t,x);ylabel('幅度');xlabel('n');subplot(2,1,2)stem(t,y);ylabel('幅度');xlabel('n');3.结果4.结果分析M 2.41.题目(a)Write a matlab program to generate a sinusoidal sequence x[n]= Acos(ω0 n+Ф) and plot thesequence using the stem function. The input data specified by the user are the desired length L, amplitude A, the angular frequency ω0 , and the phase Фwhere 0<ω0 <pi and 0<=Ф<=2pi. Using this program generate the sinusoidal sequences shown in figure 2.15. (b)Generate sinusoidal sequences with the angular frequencies given in Problem 2.22.Determine the period of each sequence from the plot and verify the result theoretically. 2.程序%用户定义的参数L = input('Desired length = ');A = input('Amplitude = ');omega = input('Angular frequency = ');phi = input('Phase = ');%信号产生n = 0:L-1;x = A*cos(omega*n + phi);stem(n,x);xlabel('n');ylabel('幅度');title(['\omega_{o} = ',num2str(omega)]);3.结果(a)ω0=0ω0=0.1πω0=0.8πω0=1.2π(b)ω0=0.14πω0=0.24πω0=0.34πω0=0.68πω0=0.75π4.结果分析M 2.51.题目Generate the sequences of problem 2.21(b) to 2.21(e) using matlab.2.程序(b)n = 0 : 99;x=sin(0.6*pi*n+0.6*pi);stem(n,x);xlabel('n');ylabel('幅度');(c)n = 0 : 99;x=2*cos(1.1*pi*n-0.5*pi)+2*sin(0.7*pi*n);stem(n,x);xlabel('n');ylabel('幅度');(d)n = 0 : 99;x=3*sin(1.3*pi*n-4*cos(0.3*pi*n+0.45*pi));stem(n,x);xlabel('n');ylabel('幅度');(e)n = 0 : 99;x=5*sin(1.2*pi*n+0.65*pi)+4*sin(0.8*pi*n)-cos(0.8*pi*n);stem(n,x);xlabel('n');ylabel('幅度');(f)n = 0 : 99;x=mod(n,6);stem(n,x);xlabel('n');ylabel('幅度');3.结果(b)(c)(d)(e)(f)4.结果分析M 2.61.题目Write a matlab program to plot a continuous-time sinusoidal signal and its sampled version and verify figure 2.19. You need to use the hold function to keep both plots.2.程序%用户定义的参数fo = input('Frequency of sinusoid in Hz = ');FT = input('Samplig frequency in Hz = ');%产生信号t = 0:0.001:1;g1 = cos(2*pi*fo*t);plot(t,g1,'-')xlabel('时间t');ylabel('幅度')holdn = 0:1:FT;gs = cos(2*pi*fo*n/FT);plot(n/FT,gs,'o');hold off3.结果4.结果分析M 3.11.题目Using program 3_1 determine and plot the real and imaginary parts and the magnitude and phase spectra of the following DTFT for various values of r and θ:G(e jω)=1, 0<r<1.1−2r(cosθ)e−jω+r2e−2jω2.程序%program 3_1%discrete-time fourier transform computatition%k=input('Number of frequency points = ');num=input('Numerator coefficients= ');den=input('Denominator coefficients= ');%computer the frequency responsew=0:pi/k:pi;h=freqz(num,den,w);%plot the frequency responsesubplot(221)plot(w/pi,real(h));gridtitle('real part')xlabel('\omega/\pi');ylabel('Amplitude') subplot(222)plot(w/pi,imag(h));gridtitle('imaginary part')xlabel('\omega/\pi');ylabel('Amplitude') subplot(223)plot(w/pi,abs(h));gridtitle('magnitude spectrum')xlabel('\omega/\pi');ylabel('magnitude') subplot(224)plot(w/pi,angle(h));gridtitle('phase spectrum')xlabel('\omega/\pi');ylabel('phase,radians')3.结果(a)r=0.8 θ=π/6(b)r=0.6 θ=π/34.结果分析M 3.41.题目Using matlab verify the following general properties of the DTFT as listed in Table 3.2:(a) Linearity, (b) time-shifting, (c) frequency-shifting, (d) differentiation-in-frequency, (e) convolution, (f) modulation, and (g) Parseval’s relation. Since all data in matlab have to be finite-length vectors, the sequences to be used to verify the properties are thus restricted to be of finite length.2.程序%先定义两个信号N = input('The length of the sequence = ');k = 0:N-1;%g为正弦信号g = 2*sin(2*pi*k/(N/2));%h为余弦信号h = 3*cos(2*pi*k/(N/2));[G,w] = freqz(g,1);[H,w] = freqz(h,1);%*************************************************************************%% 线性性质alpha = 0.5;beta = 0.25;y = alpha*g+beta*h;[Y,w] = freqz(y,1);figure(1);subplot(211),plot(w/pi,abs(Y));xlabel('\omega/\pi');ylabel('|Y(e^j^\omega)|');title('线性叠加后的频率特性');grid;% 画出Y 的频率特性subplot(212),plot(w/pi,alpha*abs(G)+beta*abs(H));xlabel('\omega/\pi');ylabel('\alpha|G(e^j^\omega)|+\beta|H(e^j^\omega)|');title('线性叠加前的频率特性');grid;% 画出alpha*G+beta*H 的频率特性%*************************************************************************% % 时移性质n0 = 10;%时移10个的单位y2 = [zeros([1,n0]) g];[Y2,w] = freqz(y2,1);G0 = exp(-j*w*n0).*G;figure(2);subplot(211),plot(w/pi,abs(G0));xlabel('\omega/\pi');ylabel('|G0(e^j^\omega)|');title('G0的频率特性');grid;% 画出G0的频率特性subplot(212),plot(w/pi,abs(Y2));xlabel('\omega/\pi');ylabel('|Y2(e^j^\omega)|');title('Y2的频率特性');grid;% 画出Y2 的频率特性%*************************************************************************% % 频移特性w0 = pi/2; % 频移pi/2r=256; %the value of w0 in terms of number of samplesk = 0:N-1;y3 = g.*exp(j*w0*k);[Y3,w] = freqz(y3,1);% 对采样的512个点分别进行减少pi/2,从而生成G(exp(w-w0))k = 0:511;w = -w0+pi*k/512;G1 = freqz(g,1,w);figure(3);subplot(211),plot(w/pi,abs(Y3));xlabel('\omega/\pi');ylabel('|Y3(e^j^\omega)|');title('Y3的频率特性');grid;% 画出Y3的频率特性subplot(212),plot(w/pi,abs(G1));xlabel('\omega/\pi');ylabel('|G1(e^j^\omega)|');title('G1的频率特性');grid;% 画出G1 的频率特性%*************************************************************************% % 频域微分k = 0:N-1;y4 = k.*g;[Y4,w] = freqz(y4,1);%在频域进行微分y0 = ((-1).^k).*g;G2 = [G(2:512)' sum(y0)]';delG = (G2-G)*512/pi;figure(4);subplot(211),plot(w/pi,abs(Y4));xlabel('\omega/\pi');ylabel('|Y4(e^j^\omega)|');title('Y4的频率特性');grid;% 画出Y4的频率特性subplot(212),plot(w/pi,abs(delG));xlabel('\omega/\pi');ylabel('|delG(e^j^\omega)|');title('delG的频率特性');grid;% 画出delG的频率特性%*************************************************************************% % 相乘性质y5 = conv(g,h);%时域卷积[Y5,w] = freqz(y5,1);figure(5);subplot(211),plot(w/pi,abs(Y5));xlabel('\omega/\pi');ylabel('|Y5(e^j^\omega)|');title('Y5的频率特性');grid;% 画出Y5的频率特性subplot(212),plot(w/pi,abs(G.*H));%频域乘积xlabel('\omega/\pi');ylabel('|G.*H(e^j^\omega)|');title('G.*H的频率特性');grid;% 画出G.*H的频率特性%*************************************************************************% % 帕斯瓦尔定理y6 = g.*h;%对于freqz函数，在0到2pi直接取样[Y6,w] = freqz(y6,1,512,'whole');[G0,w] = freqz(g,1,512,'whole');[H0,w] = freqz(h,1,512,'whole');% Evaluate the sample value at w = pi/2% and verify with Y6 at pi/2H1 = [fliplr(H0(1:129)') fliplr(H0(130:512)')]';val = 1/(512)*sum(G0.*H1);% Compare val with Y6(129) i.e sample at pi/2 % Can extend this to other points similarly% Parsevals theoremval1 = sum(g.*conj(h));val2 = sum(G0.*conj(H0))/512;% Comapre val1 with val23.结果(a)(b)(c)(d)(e)4.结果分析M 3.81.题目Using matlab compute the N-point DFTs of the length-N sequences of Problem 3.12 for N=3, 5, 7, and 10. Compare your results with that obtained by evaluating the DTFTs computed in Problem 3.12 at ω= 2pik/N, k=0, 1,……N-1.2.程序%用户定义N的长度N = input('The value of N = ');k = -N:N;y1 = ones([1,2*N+1]);w = 0:2*pi/255:2*pi;Y1 = freqz(y1, 1, w);%对y1做傅里叶变换Y1dft = fft(y1);k = 0:1:2*N;plot(w/pi,abs(Y1),k*2/(2*N+1),abs(Y1dft),'o');grid;xlabel('归一化频率');ylabel('幅度');(a)clf;N = input('The value of N = ');k = -N:N;y1 = ones([1,2*N+1]);w = 0:2*pi/255:2*pi;Y1 = freqz(y1, 1, w);Y1dft = fft(y1);k = 0:1:2*N;plot(w/pi,abs(Y1),k*2/(2*N+1),abs(Y1dft),'o');xlabel('Normalized frequency');ylabel('Amplitude');(b)%用户定义N的长度N = input('The value of N = ');k = -N:N;y1 = ones([1,2*N+1]);y2 = y1 - abs(k)/N;w = 0:2*pi/255:2*pi;Y2 = freqz(y2, 1, w);%对y1做傅里叶变换Y2dft = fft(y2);k = 0:1:2*N;plot(w/pi,abs(Y2),k*2/(2*N+1),abs(Y2dft),'o');grid;xlabel('归一化频率');ylabel('幅度');(c)%用户定义N的长度N = input('The value of N = ');k = -N:N;y3 =cos(pi*k/(2*N));w = 0:2*pi/255:2*pi;Y3 = freqz(y3, 1, w);%对y1做傅里叶变换Y3dft = fft(y3);k = 0:1:2*N;plot(w/pi,abs(Y3),k*2/(2*N+1),abs(Y3dft),'o');grid;xlabel('归一化频率');ylabel('幅度');3.结果(a)N=3N=5 N=7N=10 (b)N=3N=5 N=7N=10 (c)N=3N=5 N=7N=104.结果分析M 3.191.题目Using Program 3_10 determine the z-transform as a ratio of two polynomials in z-1 from each of the partial-fraction expansions listed below:(a)X1(z)=−2+104+z−1−82+z−1,|z|>0.5,(b)X2(z)=3.5−21−0.5z−1−3+z−11−0.25z−2,|z|>0.5,(c)X3(z)=5(3+2z−1)2−43+2z−1+31+0.81z−2,|z|>0.9,(d)X4(z)=4+105+2z−1+z−16+5z−1+z−2,|z|>0.5.2.程序% Program 3_10% Partical-Fraction Expansion to rational z-Transform %r = input('Type in the residues = ');p = input('Type in the poles = ');k = input('Type in the constants = ');[num, den] = residuez(r,p,k);disp('Numberator polynominal coefficients');disp(num) disp('Denominator polynomial coefficients'); disp(den)4.结果分析M 4.61.题目Plot the magnitude and phase responses of the causal IIR digital transfer functionH(z)=0.0534(1+z−1)(1−1.0166z−1+z−2) (1−0.683z−1)(1−1.4461z−1+0.7957z−2).What type of filter does this transfer function represent? Determine the difference equation representation of the above transfer function.2.程序b=[0.0534 -0.00088644 -0.00088644 0.0534];a=[1 -2.1291 1.7833863 -0.5434631];figure(1)freqz(b,a);figure(2)[H,w]=freqz(b,a);plot(w/pi,abs(H)),grid;xlabel('Normalized Frequency (\times\pi rad/sample)'),ylabel('Magnitude');幅度化成真值之后：4.结果分析H(z)=0.0534−0.00088644z−1−0.00088644z−2+0.0534z−31−2.1291z−1+1.7833863z−2−0.5434631z−3M 4.71.题目Plot the magnitude and phase responses of the causal IIR digital transfer functionH(z)=(1−z−1)4(1−1.499z−1+0.8482z−2)(1−1.5548z−1+0.6493z−2).2.程序b=[1 -4 6 -4 1];a=[1 -3.0538 3.8227 -2.2837 0.5472]; figure(1)freqz(b,a);figure(2)[H,w]=freqz(b,a);plot(w/pi,abs(H)),grid;xlabel('Normalized Frequency (\times\pi rad/sample)'), ylabel('Magnitude');3.结果4.结果分析。
MATLAB 上机 习题及答案
15、今有多项式P1(x)=x4-2x+1,P2(x)=x2+4x-0.5，要求先求得P(x)=P1(x)+P2(x)，然后计算xi=0.2*i各点上的P(xi)（i=0,1,2,…,5）值。
p1=[1.0 0.0 0.0 -2.0 1.0];>> p2=[0.0 0.0 1.0 4.0 -0.5];>> p1x=poly2sym(p1);p2x=poly2sym(p2);>> p=p1x+p2xp =x^4+2*x+1/2+x^2>> x=0:5;>> x.^4+2*x+1/2+x.^2ans =0.5000 4.5000 24.5000 96.5000 280.5000 660.50001、试个MATLAB的工作空间中建立以下2个矩阵：A＝[1 2]1234B⎡⎤=⎢⎥⎣⎦，求出矩阵A和B的乘积，并将结果赋给变量C。
>> A=[1 2]A =1 2>> B=[1 23 4]B =1 23 4>> C=A*BC =7 102、利用MATLAB提供的帮助信息，了解inv命令的调用格式，并作简要说明。
help invINV Matrix inverse.INV(X) is the inverse of the square matrix X.A warning message is printed if X is badly scaled ornearly singular.See also SLASH, PINV, COND, CONDEST, LSQNONNEG, LSCOV. Overloaded methodshelp gf/inv.mhelp zpk/inv.mhelp tf/inv.mhelp ss/inv.mhelp lti/inv.mhelp frd/inv.mhelp sym/inv.mhelp idmodel/inv.m3、使用help命令查询函数plot的功能以及调用方法，然后利用plot命令绘制函数y=sin(x)的图形，其中0xπ≤≤。
matlab上机习详解_试题答案
P 第一次实验答案1． 设要求以0.01秒为间隔，求出y 的151个点，并求出其导数的值和曲线。
clc clearx=0:0.01:1.5;y=sqrt(3)/2*exp(-4*x).*sin(4*sqrt(3)*x+pi/3)y1=diff(y) subplot(2,1,1) plot(x,y)subplot(2,1,2) plot(x(1:150),y1)2绘制极坐标系下曲线（a,b,n 自定数据）clc clear a=10; b=pi/2; n=5;theta=0:pi/100:2*pi; rho=a*cos(b+n*theta); polar(theta,rho)3. 列出求下列空间曲面交线的程序clc clearx=[-5:0.5:5];[X,Y]=meshgrid(x); z1=X.^2-2*Y.^2;z2=X.*2-Y.*3; xlabel('x') ylabel('y') zlabel('z') surf(X,Y,z1) hold onsurf(X,Y,z2)k=find(abs(z1-z2)<0.5); x1=X(k) y1=Y(k)z3=x1.^2-2*y1.^2 hold onplot3(x1,y1,z3,'*')4、设 ⎥⎦⎤⎢⎣⎡++=)1(sin 35.0cos 2x x x y 把x=0~2π间分为101点，画出以x 为横坐标，y为纵坐标的曲线，要求有图形标注。
clc clearx=-2*pi:0.1: 2*pi;y=cos(x).*(0.5+sin(x)*3./(1+x.^2)); plot(x,y,'b*-'); title('绘图'); xlabel('x 坐标'); ylabel('y 坐标'); legend('原函数')gtext('y=cos(x)(0.5+3*sin(x)/(1+x^2))')5、求下列联立方程的解 81025695832475412743-=+-+-=-+-=++-=--+w z y x w z x w z y x w z y xclc cleara=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10];b=[4,-3,9,-8]; c=b/a; x=c(1,1) y=c(1,2) z=c(1,3) w=c(1,4)6. 假设一曲线数据点为x = 0:2:4*pi;y = sin(x).*exp(-x/5);试将x 的间距调成 0.1，采用不同插值方法进行插值，并通过子图的形式将不同插值结果和原始数据点绘制在同一图形窗口。
Matlab上机练习题及答案
Matlab 上机练习题及答案---------------------------------------------------------------------1、矩阵Y=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡3472123100451150425，给出元素1的全下标和单下标，并用函数练习全下标和单下标的转换，求出元素100的存储位置。
取出子矩阵⎥⎦⎤⎢⎣⎡21301，并求该矩阵的维数。
解：命令为：Y=[5,2,4;0,15,1;45,100,23;21,47,3]Y(2,3)Y(10)sub2ind([43],2,3)[i,j]=ind2sub([43],10)find(Y==100)sub2ind([43],3,2)B=Y(2:2:4,3:-2:1)或B=Y([24],[31])[m n]=size(Y)---------------------------------------------------------------------2、已知矩阵A=[10-1；241；-205]，B=[0-10；213；112]求2A+B 、A 2-3B 、A*B 、B*A 、A.*B ，A/B 、A\B 解：命令为：A=[10-1；241；-205]B=[0-10；213；112]E=2*A+B F=A^2-3*B G=A*B H=B*A I=A.*B J=A/B K=A\B---------------------------------------------------------------------3、利用函数产生3*4阶单位矩阵和全部元素都为8的4*4阶矩阵，并计算两者的乘积。
解：命令为：A=eye(3,4)B=8*ones(4)C=A*B---------------------------------------------------------------------4、创建矩阵a=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡------7023021.5003.120498601,取出其前两列构成的矩阵b ，取出前两行构成矩阵c ，转置矩阵b构成矩阵d ，计算a*b 、c<d ,c&d,c|d ,~c|~d 解：命令为：a=[-1,0,-6,8;-9,4,0,12.3;0,0,5.1,-2;0,-23,0,-7]b=a(:,[12])c=a([12],:)d=b’e=a*b f=c<d g=c&d h=c|d i=~c|~d---------------------------------------------------------------------5、求!201∑=n n 解：命令文件为sum=0;s=1;for n=1:20s=n*s;sum=sum+s;end sum---------------------------------------------------------------------6、求a aa aaa aa a S n ++++=得值，其中a 是一个数字，由键盘输入，表达式中位数最多项a 的个数，也由键盘输入。
MATLAB上机实验练习题及答案
MATLAB上机实验练习题及答案09级MATLAB上机实验练习题1、给出一个系数矩阵A[2 3 4;5 4 1;1 3 2]，U=[1 2 3]，求出线性方程组的一个精确解。
2、给出两组数据x=[0 0.3 0.8 1.1 1.6 2.3]’y=[0.82 0.72 0.63 0.60 0.55 0.50]’，我们可以简单的认为这组数据在一条衰减的指数函数曲线上，y=C1+C2e-t通过曲线拟合求出这条衰减曲线的表达式，并且在图形窗口画出这条曲线，已知的点用*表示。
3、解线性方程4、通过测量得到一组数据：5、已知一组测量值6、从某一个过程中通过测量得到：分别采用多项式和指数函数进行曲线拟合。
7、将一个窗口分成四个子窗口，分别用四种方法做出多峰函数的表面图（原始数据法，临近插值法，双线性插值法，二重三次方插值法）8、在同一窗口使用函数作图的方法绘出正弦、余弦、双曲正弦、双曲余弦。
分别使用不同的颜色，线形和标识符。
9、下面的矩阵X表示三种产品五年内的销售额，用函数pie显示每种产品在五年内的销售额占总销售额的比例，并分离第三种产品的切片。
X= 19.3 22.1 51.634.2 70.3 82.4 61.4 82.9 90.8 50.5 54.9 59.1 29.4 36.3 47.010、对应时间矢量t ，测得一组矢量y采用一个带有线性参数的指数函数进行拟合，y=a 0+a 1e -t +a 2te -t ，利用回归方法求出拟合函数，并画出拟合曲线，已知点用圆点表示。
11、请创建如图所示的结构数组（9分）12、创建如图所示的元胞数组。
（9分）13、某钢材厂从1990年到2010年的产量如下表所示，请利用三次样条插值的方法计算1999年该钢材厂的产量，并画出曲线，已知数据用‘*’表示。
要求写出达到题目要求的MATLAB 操作过程，不要求计算结果。
14、在一次化学动力学实验中，在某温度下乙醇溶液中，两种化合物反应的产物浓度与反应时间关系的原始数据如下，请对这组数据进行三次多项式拟合，并画出拟合曲线，已知数据如下。
西电MATLAB大作业
MATLAB大作业班级：021131 姓名：XXX 学号：02113XXX实验一傅里叶变换一试验要求信号在传输过程中，由于受信道或环境等客观因素的影响，在接收端得到的是噪声环境下的信号。
利用FFT函数对这一信号进行傅里叶分析，确定原始有用信号的频率。
已知：采样率：fs=100MHz，有用信号频率：f0=10MHz，样本数：1024点，信噪比：S/N=5dB。
二设计思路1、根据已知条件产生一个有用信号2、根据已知的信噪比产生一个高斯噪声3、将产生的有用信号和高斯噪声叠加，将此信号视为被噪声污染的信号4、调用FFT函数对被噪声污染的信号进行快速傅里叶变换5、对比图象并结合已知条件，得出原始有用信号的频率设计思路如下图：三程序代码clear allclose allclcfs=100e6; %采样率fs=100Mhzf0=10e6; %有用信号频率f0=10MhzN=1024; %样本数N=1024n=0:N-1;t=[0:N-1]/fs; %时间序列input=2*sin(2*pi*f0*t)+sin(2*pi*f0*t)-1; %输入信号subplot(2,2,1),plot(t,input); %有用信号随时间变化的输振幅xlabel('时间');ylabel('振幅');title('输入信号');grid on;am=2; %有用信号的最大幅度p=10^(5/20); %信噪比：S/N=5dBb=am/10^(p/20); %b为噪声信号最大幅度,p=20*log10(am/b)noise=b*randn(1,N); %高斯噪声信号subplot(2,2,2),plot(t,noise); %绘出随时间变化的噪声信号的振幅xlabel('时间');ylabel('振幅');title('噪声信号');grid on;sig=input+noise; %被噪音污染的信号subplot(2,2,3),plot(t,sig); %绘出随时间变化的被噪声信号污染后的信号的振幅xlabel('时间');ylabel('振幅');title('被噪声信号污染后的信号');grid on;y=fft(sig,N); %对被噪声信号污染后的信号进行快速傅里叶变换mag=abs(y); %求得傅里叶变换后的振幅f=n*fs/N; %频率序列subplot(2,2,4),plot(f,mag); %绘出随频率变化的振幅xlabel('频率/Hz');ylabel('振幅');title('傅里叶变换后的信号');grid on;四实验结果及结果分析1、实验结果2、结果分析由“傅里叶变换后的信号”图知有用信号的频率是1×107Hz=10MHz=f0,符合实验设定的事实，实验成功。
MATLAB上机习题一
MATLAB上机习题一请按以下步骤完成上机实验：1）在FTP上下载“MATLAB上机习题一.doc”文件，所有习题列在该文件内；2）在MATLAB中完成所有习题，并将屏幕截图粘贴到相应习题后面；3）如果习题是问答题，请将答案写在题目后；4）如果有的习题要求提供脚本文件，请将脚本文件内容拷贝到相应习题后；5）将文件保存并重命名为“自己的学号-姓名”，例如“20110771-张三.doc”；6）上传该文件到FTP的相关目录。
1. 运行MATLAB软件，拖放、关闭界面上的子窗口，并恢复到原始试图。
注：不需截图，但要说明如何操作的2. 采用鼠标及命令两种方式将桌面添加到MATLAB搜索路径列表的起始及最后位置。
3. 采用鼠标及命令方式将当前工作目录设置为桌面指向的文件夹。
4. 通过使用帮助确定内置变量ispc 的含义。
只需用简单语言说明。
6. 观察MATLAB中关键字、字符串、注释的字体显示颜色。
只需用简单语言说明。
5. 创建double类型的变量，并进行计算1）a=87，b=190，计算a+b 、a-b 、a*b ；2）创建uint8类型的变量，数值与（1）中相同，进行相同的计算，观察计算结果与预想的是否一致，并说明为什么。
6. 计算如下表达式：1）()sin 602）3e3）3cos 4π⎛⎫ ⎪⎝⎭4）27562323336+⨯-⨯⨯7. 设u=2，v=3，计算：1）4log uv v2）()22u e v v u+- 38. 计算如下表达式：1）()()3542i i -+2）()sin 28i -9. 判断下面语句的运算结果，并与MATLAB 计算结果做比较：1）420<2）420≤3）420==4）4~20=5）''''b B <10. 设a=39，b=58，c=3，d=7，判断下面表达式的值与MATLAB 计算结果做比较： 1）a b >2）a b <3）&&a b b c >>4）a d ==5）|a b c >6）~~d11. 判断下列变量名哪些是合法的MATLAB 变量：1）fred7）fred! 2）book_18）book-1 3）2ndplace9）Second_Place 4）#110）No_1 5）vel_511）vel.5 6）tan12）while 提示，可以使用isvarname 对上述变量名进行检验。
MATLAB上机练习一参考解答
上机练习一参考解答一、实验目的1、 熟悉Matlab 编程2、 体会数学上恒等，算法上不一定恒等二、实验内容1． Using the Taylor polynomial of degree nine and three-digit rounding arithmetic to find an approximationto 5-e by each of the following methods.(A) ∑=--≈905!)5(n n n e , (B) ∑=-≈=9055!5/11n nn e e An approximate value of 5-e correct to three digits is 31074.6-⨯. Which formula, (A) or (B), gives the most accuracy, and why?1） 算法基础利用x e 的Taylor 公式00!!n nk x n n x x e n n ∞===≈∑∑，x -∞<<+∞ （1）及001/1/1/!!n nk x x n n x x e e n n ∞-====≈∑∑，x -∞<<+∞， （2）其中k 是根据精度要求给定的一个参数。
在本题中将k 取为9， x 取为－5或5即可由公式（1）或（2）得到5-e 的近似计算方法（A ）或（B ）。
2） 程序下述程序用公式（A ）及（B ）分别在Matlab 许可精度下及限定在字长为3的算术运算情况下给出5-e的近似计算结果，其中results_1, results_2为用方法（A ）在上述两种情况下的计算结果，err_1, err_2为相应的绝对误差；类似的，results_3, results_4为用方法（B ）在上述两种情况下的计算结果，err_3, err_4为相应的绝对误差；具体程序如下：% Numerical Experiment 1.1 % by Xu Minghua, May 17, 2008 clc; %Initialize the data x=-5; k=9; m=3; %three-digit rounding arithmetic %------------------------------------ % Compute exp(x) by using Method (A) % with the computer precision results_1=1; power_x=1; for i=1:k factor_x=x/i; power_x=power_x*factor_x; results_1=results_1+power_x; end results_1 err_1=abs(exp(x)-results_1)%------------------------------------% Compute exp(x) by using Method (A) % with the 3-digits precisionresults_2=1;power_x=1;for i=1:kfactor_x=digit(x/i,m);power_x=digit(power_x*factor_x,m); results_2=digit(results_2+power_x,m); endresults_2err_2=abs(exp(x)-results_2)%------------------------------------% Compute exp(x) by using Method (B) % with the computer precisiont=-x;results_3=1;power_x=1;for i=1:kfactor_x=t/i;power_x=power_x*factor_x;results_3=results_3+power_x; endresults_3=1/results_3err_3=abs(exp(x)-results_3)%------------------------------------% Compute exp(x) by using Method (B) % with the 3-digits precisiont=-x; results_4=1;power_x=1;for i=1:kfactor_x=digit(t/i,m);power_x=digit(power_x*factor_x,m);results_4=digit(results_4+power_x,m); endresults_4=digit(1/results_4,m)err_4=abs(exp(x)-results_4)%------------------------------------上述主程序用到一个子程序digit.m, digit(x,m)的作用是将x四舍五入成m位数。
MATLAB上机实验1答案(2021年整理精品文档)
MATLAB上机实验1答案编辑整理：尊敬的读者朋友们：这里是精品文档编辑中心，本文档内容是由我和我的同事精心编辑整理后发布的，发布之前我们对文中内容进行仔细校对，但是难免会有疏漏的地方，但是任然希望（MATLAB上机实验1答案）的内容能够给您的工作和学习带来便利。
同时也真诚的希望收到您的建议和反馈，这将是我们进步的源泉，前进的动力。
本文可编辑可修改，如果觉得对您有帮助请收藏以便随时查阅，最后祝您生活愉快业绩进步，以下为MATLAB上机实验1答案的全部内容。
实验1 Matlab 初步一、问题已知矩阵A 、B 、b 如下：⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡-------------=031948118763812654286174116470561091143A ⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎣⎡------=503642237253619129113281510551201187851697236421B []1187531=b应用Matlab 软件进行矩阵输入及各种基本运算.二、实验目的学会使用Matlab 软件构作已知矩阵对应的行（列)向量组、子矩阵及扩展矩阵，实施矩阵的初等变换及线性无关向量组的正交规范化,确定线性相关相关向量组的一个极大线性无关向量组，且将其余向量用极大线性无关向量组线性表示，并能编辑M 文件来完成所有的实验目的。
三、预备知识1、 线性代数中的矩阵及其初等变换、向量组的线性相关性等知识。
2、 Matlab 软件的相关命令提示如下；（1） 选择A 的第i 行做一个行向量：ai=A(i,:）;（2）选择A的第j行做一个列向量：ai=A(j,:）;（3）选择A的某几行、某几列上的交叉元素做A的子矩阵：A(［行号]，[列号］）；（4）n阶单位阵：eye(n);n阶零矩阵：zeros(n）；（5）做一个n维以0或1为元素的索引向量L，然后取A（:,L）,L中值为1的对应的列将被取到。
MATLAB上机习题一 解答
9. 判断下面语句的运算结果，并与 MATLAB 计算结果做比较： 1） 4 < 20 2） 4 ≤ 20 3） 4 == 20 4） 4 ~ = 20 5） ' b ' < ' B ' 解答： 1）true 2）true 3）false 4）true 5）false
注意，字符比较实际上比较的是字符的 ASCII 码的大小。
提示，可以使用 isvarname 对上述变量名进行检验。
解答：1、2、5、6、9、10 是合法的，其它不能作为变量名。 注意，第 6 问的函写脚本程序，命名为 barbell.m，完成如下计算： 如图 1 所示的杠铃，
每个球的半径是 10cm，两个求直接链接杆的长度是 15cm，杆的直径为 1cm，计 算杠铃的体积和表面积；使用 disp 函数显示计算结果。 （提示：disp 函数的用法 参考教材的例 2-31， 面积计算时相交部分球冠可近似为圆， 体积计算时可忽略重 叠的球冠体积） 解答： barbell.m（拷贝文件内容到如下方框中） ： %注意，将圆柱与球相交部分面积近似为圆形 clear all clc R=10; %球半径,单位cm l=15; %杆长度,单位cm d=1; %杆直径,单位为cm
n2a + p nRT (V − nb ) = V2
该方程新增了 a 和 b 两个变量，用来表示纯净气体的属性。 1）编写脚本程序，命名为 state_equation.m，利用如下数据： 压强，p 摩尔数，n 体积，V a b 理想气体常数，R 220 mbar 2 mol 1L 5.536 L2bar/mol2 0.03049 L/mol 8.314 J/(mol.K)
1. 运行 MATLAB 软件，拖放、关闭界面上的子窗口，并恢复到原始视图。 解答：略。注意采用下图方式恢复视图
Matlab上机题库及详细答案
case num2cell(25:49) %价格大于等于2500但小于5000
Байду номын сангаасrate=10/100;
otherwise %价格大于等于5000
rate=14/100;
end
price=price*(1-rate) %输出商品实际销售价格
例19已知，当n=100时，求的值。
解：程序如下：
y=0;n=100;for i=1:n；y=y+1/(2*i-1);End
例20：一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数
解：for m=100:999
m1=fix(m/100);m2=rem(fix(m/10),10);m3=rem(m,10);
if m==m1*m1*m1+m2*m2*m2+m3*m3*m3；disp(m)；end
D=diag(1:5);D*A %用D左乘A，对A的每行乘以一个指定常数
例15：输入x,y的值，并将它们的值互换后输出。
程序如下：
x=input('Input x please.')；y=input('Input y please.');z=x;x=y;y=z;
disp(x);disp(y);
例16：:求一元二次方程ax2 +bx+c=0的根。
switch fix(price/100)
case {0,1} %价格小于200
rate=0;
case {2,3,4} %价格大于等于200但小于500
rate=3/100;
case num2cell(5:9) %价格大于等于500但小于1000
(完整word版)2018秋西安电子科技大学数字图像处理上机实验题答案及matlab程序
西安电子科技大学数字图像处理上机实验题李XX2018-12-16数字图像处理上机实验题1. 产生右图所示图像 f1(m,n)，其中图像大小为256×256，中间亮条为128×32，暗处=0，亮处=100。
对其进行FFT：(matlab程序附于文末附录)①同屏显示原图f1(m,n)和FFT(f1)的幅度谱图；②若令f2(m,n)=(-1)m+n f1(m,n)，重复以上过程，比较二者幅度谱的异同，简述理由；结论：相同点：图像f1,f2的幅度谱的实质没有改变，幅度等都没有发生变化。
不同点：f2的频谱是对f1频谱的移位，它是f1的频谱从原点(0,0)移动到中心点(64,64)得到的频谱。
原因：f2图像是对f1图像做频谱中心化变换得来的，这一过程不改变图像的幅度等特性。
③若将f2(m,n)顺时针旋转90 度得到f3(m,n)，试显示FFT(f3)的幅度谱，并与FFT(f2)的幅度谱进行比较；结论：由图可知，FFT(f3)的幅度谱是由FFT(f2)的幅度谱顺时针旋转90度得到的。
④若将f1(m,n) 顺时针旋转90 度得到f4(m,n)，令f5(m,n)=f1(m,n)+f4(m,n)，试显示FFT(f5)的幅度谱，并指出其与 FFT(f1)和FFT(f4)的关系；结论：FFT(f5)的幅度谱是FFT(f1)和FFT(f4)幅度谱叠加的结果。
⑤若令f6(m,n)=f2(m,n)+f3(m,n)，试显示FFT(f6)的幅度谱，并指出其与FFT(f2)和FFT(f3)的关系，比较FFT(f6)和FFT(f5)的幅度谱。
结论：FFT(f6)的幅度谱是FFT(f2)和FFT(f3)幅度谱叠加的结果。
FFT(f6)是FFT(f5)频谱中心化的结果。
2. 产生教材104 页题图4.18（右图）所示的二值图像（白为1，黑为0），编程实现习题4.18 所要求的处理（3×3 的平均滤波和中值滤波）功能（图像四周边界不考虑，处理结果按四舍五入仍取0 或1），显示处理前后的图像，比较其异同。
最新matlab 上机练习1-2 (2)
第一次上机:基础准备及入门2.1 指令窗运行入门2.1.1 最简单的计算器使用法【例2.5.1-1】求23)]47(212[÷-⨯+的算术运算结果。
（1）用键盘在MATLAB 指令窗中输入以下内容 >> (12+2*(7-4))/3^2（2）在上述表达式输入完成后，按【Enter 】键，该就指令被执行。
（3）在指令执行后，MATLAB 指令窗中将显示以下结果。
ans = 2【例2.5.1-2】简单矩阵⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=987654321A 的输入步骤。
（1）在键盘上输入下列内容A = [1,2,3; 4,5,6; 7,8,9]（2）按【Enter 】键，指令被执行。
（3）在指令执行后，MATLAB 指令窗中将显示以下结果：A =1 2 3 4 5 6 7 8 9【例2.5.1-3】矩阵的分行输入A=[1,2,3 4,5,6 7,8,9]（以下是显示结果）A =1 2 3 4 5 6 7 8 9【例2.5.1-4】指令的续行输入（以下格式在除Notebook 外的MATLAB 环境中可运行）S = 1 – 1/2 + 1/3 –1/4 + 1/5 – 1/6 ... + 1/7 – 1/8 S =0.63452.1.2数值、变量和表达式2.1.2.1复数和复数矩阵【例2.5.2.5-1】复数ieziziz63212,21,43π=+=+=表达，及计算321zzzz=。
（1）经典教科书的直角坐标表示法z1= 3 + 4iz1 =3.0000 +4.0000i（2）采用运算符构成的直角坐标表示法和极坐标表示法z2 = 1 + 2 * i %运算符构成的直角坐标表示法z3=2*exp(i*pi/6) %运算符构成的极坐标表示法z=z1*z2/z3z2 =1.0000 +2.0000iz3 =1.7321 + 1.0000iz =0.3349 + 5.5801i【例2.5.2.5-2】复数矩阵的生成及运算A=[1,3;2,4]-[5,8;6,9]*iB=[1+5i,2+6i;3+8*i,4+9*i]C=A*BA =1.0000 - 5.0000i 3.0000 - 8.0000i2.0000 - 6.0000i 4.0000 - 9.0000iB =1.0000 + 5.0000i2.0000 + 6.0000i3.0000 + 8.0000i4.0000 + 9.0000iC =1.0e+002 *0.9900 1.1600 - 0.0900i1.1600 + 0.0900i 1.3700【例2.5.2.5-3】求上例复数矩阵C的实部、虚部、模和相角。
(完整版)matlab上机题目一
上机一1、（1（2） 查询所有人的成绩。
2、已知：A=[27,24,0,-1;23,35,1,14;e 4,-8,13,20;10,12,19,21;11,18,ln2,0]， 求：（1）求A 中第4、5、2行与第2、4、1列组成的C 阵，并将C 中第2行元素全改成8，得到D 阵；（2） 创建与A 同维的G 阵，使C 阵左右翻转后位于G 中第2行至第4行、第2列至第4列位置，且G 中第1行第1列元素为1，其余元素为0。
（3） 求A 阵除以3的余数阵H 。
3、已知：A=[5,17,-2,3;-1,0,21,11;14,9,-10,8]，B=[2,5,8,-1;12,-2,7,8;11,13,16,4] 求：（1）将A 除以B 阵的对应元素，并将结果朝负无穷大方向取整。
（2）求A./B ， B.\A ，A.*B 和A*B ’的值。
（3）求2*ones(3,4)-eye(3,4)|A-B>B&B-8的值.4、要求在闭区间[0，2π ]上产生50个等距采样的一维数组。
试用两种不同的指令实现。
5、对矩阵，⎥⎦⎤⎢⎣⎡=4321A 问如下MATLAB 指令：)(),(),5.0(^),5.0(.^A sqrtm A sqrt A A 所得的结果相同吗？他们中的哪个结果为复数矩阵？ 6、说出MATLAB 指令A （3，1，2，： ）=1：4所产生数组的维数、大小和长度；然后对A 进行降维处理；最后指出所有非零元素的“全下标”位置。
7、现有五个数组： (2×3)的实数数组A ； (2×3)的复数数组B ； (2×3)的字符串数组C ； (2×3)的异质（元胞）数组D ； (2×3)的结构（构架）数组E 。
问，它们各占多少字节？8、写出使以下这段文字成为字符串的MATLAB 程序。
注意保持这段文字的格式。
在英式用法中，引号通常是单引号，如‘Fire!’。
2014年下学期MATLAB语言上机练习内容(1)
要求：（1）每个同学复印一份，在做相关实验前提前做完相关题目，一定要写好相关编程。
（2) 上机时验证执行每个练习题的程序，发现问题并修改，获得正确结果。
(3) 把上机执行后的程序理解 并完整的写在课堂作业本上，以备教学检查与自己 考试复习用。
(4) 上机编程操作情况将作为期末课程考核的重要依据。
2014年下学期MATLAB 语言上机实验内容使用班级：1312201、202、203姓名： 班级： 学号：实验1: matlab 基础与符号计算1、指令>>A= [2 2 ; 2 2];>>A .^2= 。
>>A^2= 。
2、对3、数组>>x=10:-2:2; y=1:5;>>(x<5)&(y>2)，ans= 。
>> y==2|x~=2， ans= 。
3、（1）对于0>x ，求12011122+∞=∑⎪⎭⎫⎝⎛+-+k k x x k 。
（提示：理论结果为x ln ）（2）求∑∞=0k kx 的符号解，并进而用该符号解求∑∞=-0)31(k k ，∑∞=0)1(k kπ，∑∞=03k k的准确值。
（提示：注意subs 的使用）4、（1）通过符号计算求t t y sin )(=的导数dtdy。
（2）然后根据此结果，求-=0t dt dy 和2π=t dt dy。
5、指令>>A= [2 1 ; 1 3];>>A .^2-A= 。
>>A^2+A= 。
6、数组>>x=10:-2:2; y=5:-1:1;>>(x>=4)&(y>=2)，ans= 。
>> y==3|x~=3， ans= 。
7计算二重积分⎰⎰+211222)(x dydx y x 。
8在]2,0[π区间，画出dt ttx y x⎰=sin )(曲线，并计算)5.4(y 。
第一、二次上机练习
第一、二次上机练习：目的：1.Matlab 的安装、卸载：（选作，自己有电脑的可以做） 2.建立自己的工作目录，再将自己的工作目录设置到matlab 的搜索路径下。
3.掌握MATLAB 各种表达式的书写规则 4.运行课堂上讲过的例子，熟悉矩阵、表达式的基本操作和运算。
作业：1. 求下列表达式的值，显示MATLAB 工作空间的使用情况并保存全部变量：（1）)1034245.01(26-⨯+⨯=ww=sqrt(2)*(1+0.34245*10^-6)w =1.4142（2）,)tan(22ac b e abc c b a x ++-+++=ππ 其中a=3.5,b=5,c=9.8。
>> a=3.5;b=5;c=9.8;>> x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/(tan(b+c)+a)x =6.6186（3）])48333.0()41[(22απβππα---=y ，其中32.3=α，9.7-=β>> m=3.32;n=-7.9;>> y=2*pi*m^2*((1-pi/4)*n-(0.8333-pi/4)*m)y =（4）)1ln(2122t t e z t ++=，其中⎥⎦⎤⎢⎣⎡--=65.05312i t >> t=[2,1-3*i;5,-0.65];>> z=exp(2*t)*log(t+sqrt(1+t^2))/2z =1.0e+004 *0.0057 - 0.0007i 0.0049 - 0.0027i1.9884 - 0.3696i 1.7706 - 1.0539i2. 已知⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=7613870451A ，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--=023352138B ，求下列表达式的值： （1）B A 6+、I B A +-2（其中I 为单位阵）；（2）A*B 、A.*B 、B*A 、B.*A ；（3）A/B 、B\A ；（4）[A, B]、[A([1 3],:)；B^2]。
matlab上机题目一答案
上机一1、（1（2） 查询所有人的成绩。
J=struct(…number ‟,{9901,9902,9903,9904,9905},‟name ‟,{…Xiao li ‟, … …Xiao zhou ‟,‟Xiao xu ‟,‟Xiao liao ‟,‟Li na ‟},… ‟score ‟,{90,83,71,82,64},‟grade ‟,{…A ‟,‟B ‟,‟B ‟,‟A ‟,‟A ‟}); A=J(3),B=[J.score], 运行结果：2、已知：A=[27,24,0,-1;23,35,1,14;e 4,-8,13,20;10,12,19,21;11,18,ln2,0]， 求：（1）求A 中第4、5、2行与第2、4、1列组成的C 阵，并将C 中第2行元素全改成8，得到D 阵；（2） 创建与A 同维的G 阵，使C 阵左右翻转后位于G 中第2行至第4行、第2列至第4列位置，且G 中第1行第1列元素为1，其余元素为0。
（3） 求A 阵除以3的余数阵H 。
A=[27,24,0,-1;23,35,1,14;exp(4),-8,13,20;10,12,19,21;11,18,log(2),0], C=A([4,5,2],[2,4,1]),D=C;D(2,:)=8, G=zeros(size(A)); G(2:4,2:4)=fliplr(C); G(1,1)=1,H=rem(A,3)3、已知：A=[5,17,-2,3;-1,0,21,11;14,9,-10,8]，B=[2,5,8,-1;12,-2,7,8;11,13,16,4] 求：（1）将A 除以B 阵的对应元素，并将结果朝负无穷大方向取整。
（2）求A./B ， B.\A ，A.*B 和A*B ‟的值。
（3）求2*ones(3,4)-eye(3,4)|A-B>B&B-8的值.a=A./B,b= floor(A./B),c=A./B ， d=B.\A ，e=A.*B, f=A*B ‟, g=2*ones(3,4)-eye(3,4)|A-B>B&B-8,4、要求在闭区间[0，2π ]上产生50个等距采样的一维数组。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
6. 设 f ( x) e2sin( x 3) cos x e3cos x sin x ，试在[-7,7]上求出函数的零点及极大、极 小值。
解：
7. 8.
2 5 88 7 2 5 8 9 2 3 设 A 5 2 , B 10 11 7 5 ，求矩阵方程 88 6 0 3
4. 编制 m 文件，等待键盘输入，输入密码 123，密码正确，显示输 入密码正确，程序结束；否则提示，重新输入。
解：
5. 求解一元六次方程 3x6 12x5 4x4 6x2 8x 9 3 的根。 解：>> syms x
>> x=solve(3*x^6+12*x^5-4*x^4+6*x^2+8*x-12)
By:王刘鄞 1215 上机题目
1. 生成两个 3*4 的矩阵，将矩阵 A 和矩阵 B 横向拼接，并提取合并 后矩阵的主对角元素。
解：
2. 求下列矩阵的列均值、列和、列方差、所有元素的平均值，
2 3 88 7 3 1 2 5 8 9 B 2 A 5 6 4 , 10 11 7 5 。 87 93 17t=1000000; for i=1:4 for j=i+1:5 x=[d(i,1) d(j,1)]; y=[d(i,2) d(j,2)]; plot(x,y,'--') d1=sqrt((x(1)-x(2))^2+(y(1)-y(2))^2); if d1<dist dist=d1; flag=[i j]; end end end plot(d(flag,1),d(flag,2),'-.r');
解：
3.已知 5 所大学的坐标分别为： （-201.1，-99.7） 、 （0,111.5） 、 （13..8， -69.7） 、 （91.8,29.7） 、 （26.9,-199.1） 。用不同标记标出各学校，并用黄 色点划线连接两两学校；然后根据学校间的两两距离，寻找距离最短 的两所学校，并用蓝色虚线连接。 解：d=[-201.1 -99.7
AX-2A=B-X 的解。
将字符矩阵 A=['Hello!';'world';'西电科大']按行转换为一个行向量，转换为一个列向量。
解：
9. 建立一个字符串向量，然后对该向量做如下处理： (1)取第 2-7 个字符组成的子字符串。 (2)将字符串中的小写字母变成相应的大写字母，其余字符不变。 (3)统计字符串中小写字母的个数。
0 111.5 13.8 -69.7 91.8 29.7 26.9 -199.1]; plot(d(1,1),d(1,2),'*'); hold on; plot(d(2,1),d(2,2),'+'); plot(d(3,1),d(3,2),'o'); plot(d(4,1),d(4,2),'s'); plot(d(5,1),d(5,2),'d'); flag=[];