控制流与M文件.docVIP

使用循环语句计算1+2+…+500. 在matlab的M文件编辑窗口输入以下程序代码： t=0; %t的初始值 for i=1:500 %i为循环变量 t=t+i; %叠加 end t %输出t的值 运行结果为：t=125250 使用循环语句计算矩阵所有元素的和。 在matlab的M文件编辑窗口输入以下程序代码： format rat %以有理数形式输出 A=[1,42,43,46;0.5,0.96,7,8;9,-9,-5.8,5;2,2,4,4]; %矩阵A sum=0; %sum的初始值 [a,b]=size(A); %a,b分别为矩阵A的行和列 for i=1:a for j=1:b sum=sum+A(i,j); end end sum %输出sum的值 运行结果为：sum=7983/50 使用循环语句计算。（matlab的计算阶乘命令为：factorial(n) 例如：计算5!,可在matlab命令行中（或在M文件中）直接输入：factorial(5). 在matlab的M文件编辑窗口输入以下程序代码： y=0; %y的初始值 for n=1:100 %循环变量 y=y+(-1)^(n-1)*factorial(n+1)/n^(n+1); %叠加 end y %输出y的值 运行结果为：y=1142/781 使用循环语句和if语句编程：将向量a=[2 3 -1 9 3 9 8 4 6 0]由大到小进行排列组成一个新的向量。 在matlab的M文件编辑窗口输入以下程序代码： syms temp %定义变量 a=[2,3,-1,9,3,9,8,4,6,0] %向量a for i=1:length(a)-1; for j=i+1:length(a) if a(i)a(j); %判断条件 temp=a(i); a(i)=a(j); a(j)=temp; %交换数值 end end end a %输出排列好的新的向量 运行结果为：a=9 9 8 6 4 3 3 2 0 -1 已知，建立一个函数文件，并计算. 在matlab的M文件编辑窗口输入以下程序代码： function y=f(x) if x1 y=x^2+1; else y=-x-1; end 在命令窗口中输入：y1=f(3) 运行结果为：y1=10 在命令窗口中输入：y2=f(-2) 运行结果为：y2=1 编写一个函数文件，函数的自变量（即输入变量）为一元二次方程中的，函数值（即输出变量）为这个方程的两个根所组成的一个行向量，这个向量的第一个元素为，第二个元素为，并计算当时方程的根。（可参考课件及作业写法参考例题3） 在matlab的M文件编辑窗口输入以下程序代码： function y=f(a,b,c) y(1)=(-b+sqrt(b^2-4*a*c))/(2*a);????%二元一次方程的第一个根 y(2)=(-b-sqrt(b^2-4*a*c))/(2*a);????%二元一次方程的第二个根 ? 在命令窗口中输入：y=f(3,7,2) 运行结果为：y =-0.3333???-2.0000 编写一个函数文件，并使用循环语句，使得函数能计算自变量（即向量）的平均值和标准差，并将平均值和标准差构成一个向量作为函数值。其中，平均值作为这个向量的第一个元素，标准差作为这个向量的第二个元素。并计算向量x=[ 3.5 4.1 0 9.5 8.4 9 4 8 9.8 4.7 1.2 5.5 3.4 ]的平均值和标准差。（参考作业写法例题3） 其中平均值公式是： 标准差计算公式是： 在matlab的M文件编辑窗口输入以下程序代码： function y=f(x) aver=0;std=0; for i=1:length(x) ????aver=aver+x(i); end aver=aver/length(x); ? for i=1:length(x) ????std=std+(x(i)-aver)^2; end std=sqrt(std/length(x)); y(1)=aver; y(