Matlab在制导控制系统设计中的应用实验- Matlab基础及科学计算.docxVIP

Matlab基础及科学计算 1、用Matlab语句创建矩阵A，，如果给出“A(5,6)=你的学号”命令将得到什么结果？ 答:A方阵将扩充成5×6的矩阵,在第5行第6列处填充2015300005,多出来的其他未赋值处填充0,原位置处为初始值. 源代码: A=[2,3,4,5;4,3,2,1;6,4,3,7;8,9,4,6] A(5,6)=2015300005 2、用magic命令生成一个六阶的魔方矩阵，并将其全部偶数行提取出来。 答: 提取的结果为[8,30,4,32,28,36,6,2,34,26,22,12,10,14,18,24,20,16] 源代码: Magic=reshape(magic(6),[1,36]); result=[]; for i=1:length(Magic) if rem(Magic(i),2)==0 result(length(result)+1)=Magic(i); end end result 3、输入A为3×3的魔方阵，B为3×3的单位阵，由小矩阵组成3×6的大矩阵C和6×3的大矩阵D，将D矩阵的最后1列构成向量E。 答:A= B= C= D= E= 源代码: A=magic(3) B=ones(3,3) C=zeros(3,6); D=zeros(6,3); for i=1:3 C(i,:)=[A(i,:),B(i,:)]; end C D=[A;B] E=D(:,3) 4、输入矩阵A=，使用全下标方式取出元素“3”，使用单下标方式取出元素“8”，取出后2行子矩阵块，使用逻辑矩阵方式取出。 答:详细操作见代码. 源代码: A=[1,2,3;4,5,6;7,8,9] three=A(1,3) three_=A(7) last_two_line=A(2:3,:) A.*(rem(A,2)~=0).*(A~=5) 5、用数值方法求出，试不采用循环的形式给出和式的数值解。 答:S=2251799813685247 源代码: syms k symsum(2^k,k,0,50) 6、试用循环结构找出“学号后三位”以下的所有质数。 答:[2,3,5] 源代码: for i=1:005 if(isprime(i)) i end end 7、Fibonacci数列定义为，为一个数值增长很快的数列，试编写程序寻找Fibonacci数列中第一个大于10000的元素 答: 10946 源代码: i=1; while fibonacci(i)10000 i=i+1; end fibonacci(i) 8、使用while循环计算满足的最小值及相应的积。 答:mmin=3 乘积为:15 源代码: clear clc i=1; c=1; while c*(2*i-1)005 c=c*(2*i-1); i=i+1; end i c*(2*i-1) 9、下面给出一个迭代模型 写出求解该模型的M函数，如果取迭代初值为，那么请进行30000次迭代求出一组和向量，然后在所有的和坐标处点亮一个点（注意不要连线），最后绘制出所需的图形。（注：绘制出的图形又称为，它将迭代出来的随机点吸引到一起，最后得出貌似连贯的引力线图） 源代码: T9(30000) function T9(n) x=zeros(size(1,n)); y=zeros(size(1,n)); for i=1:n-1 x(i+1)=1+y(i)-1.4*x(i)^2; y(i+1)=0.3*x(i); end plot(x,y,.) title( Henon引力线图) 10、绘制出下面的图形，其中，由于该函数在t=0附近震荡剧烈，因此，尝试在在0附近进行调整步长。绘制并对比两种时间序列下的结果。第一个是采用0.02的时间间隔，另外一个在附近将时间间隔设置为0.001。 源代码: t=[-1:0.02:1]; subplot(1,2,1) plot(t,sin(1./t)) title(step=0.02) t=-0.2:0.001:0.2; subplot(1,2,2) plot(t,sin(1./t)) title(step=0.001) 11、绘出函数的图像。 源代码: a=[-10:0.01:10]; plot(a(a=1),2.*a(a=1).^2+1,a(a1a-1),0,a(a=-1),-a(a=-1).^3) 12、试绘制出的极坐标曲线，选择变量的范围为和，观察极坐标曲线的情况 源代码: subplot(1,2,1) theta=[0:0.01:4*pi]; rho=sin(8*theta/3)./(005-cos(3*theta/2).^2); p