三Matlab编程入门实例讲解和讨论

三Matlab编程入门实例讲解和讨论 Matlab编程实例讲解和讨论编程实例讲解和讨论主讲老师闫培光电子科学与技术学院主讲老师闫培光电子科学与技术学院yanpgszu.?几个函数的用法find、tic和toc，等（掌握），等（掌握）结题实例和讨论（重点）主要内容?结题实例和讨论（重点）?编写高效编程代码（掌握） 一、几个函数的用法（掌握）find:找到符合条件的元素所在的位置当作用在矩阵找到符合条件的元素所在的位置当作用在矩阵A上，它将找到A中满足条件的元素的角（是）标（是单下标形式）x=nan,nan,1,2,inf,5,nan,inf;i=find(isnan(x)|isinf(x);%Find badelements inxx(i)=;%and deletethemi=find(?isnan(x)&?isinf(x);%Find elementsthat arenot NaNand notinfinitex=x(i);%Keep thoseelementsx(isnan(x)|isinf(x)=;%Delete badelementstic和toc用来得到程序运行所花费时间 二、实例讲解和讨论例编写m函数，当输入一个变量x时能自动根据据x值来计算函数值第第11节节If的使用据据x值来计算函数值x(xA (3)&A (1)+A (3)A (2)&A (2)+A (3)A (1)p=(A (1)+A (2)+A (3)/2;s=sqrt(p*(p-A (1)*(p-A (2)*(p-A (3);disp(s);elsedisp(不能构成一个三角形。 )end例题将A矩阵用for循环写出来A=12345 二、实例讲解和讨论第2节for和while的使用12345234563456745678说明For循环更倾向于在已知计算步数的情况下使用，while更多应用在已知条件满足情况下的循环。 相比之下，更多应用在已知条件满足情况下的循环。 相比之下，while循环会更加快速的收敛。 讨论题如何编制程序？找出如何编制程序？找出A矩阵中值等于5的元素所在位置？的元素所在位置？sum=0;t=0;val=input(Enter anumber(end in0):);while(val=0)例题从键盘输入若干个数，当输入0时结束输入，求这些数的平均值和它们之和时结束输入，求这些数的平均值和它们之和sum=sum+val;t=t+1;val=input(Enter anumber(end in0):);endif(t0)summean=sum/tends=1;n=100:500;for i=1:length(n)if rem(n(i),21)=0例题求100，500之间前5个能被21整除的整数。 整除的整数。 请改成用while循环编程求解？循环编程求解？m(s)=n(i);s=s+1;continue,endif s=6breakendenddisp(m)更简洁的方法？n=100;a=ones(1,n);for i=3:n例题用for循环指令来寻找Fibona数组中第一个大于数组中第一个大于10000的元素a (1)=1,a (2)=1,a(i)=a(i-2)+a(i-1)a(i)=a(i-2)+a(i-1);if a(i)10000a(i),break;endenddisp(i)for和while可以用来迭代（重要）%求1，2.,100的和s=0;for i=1:100,s=s+i;enddisp(s);如何来求解数的阶乘？如何使用如何来求解数的阶乘？如何使用while来做循环？知已知y=1+1/3+1/5+1/2n-1，当n=100时，求y的值？程序阅读a=zeros(1,80000);%Preallocationb=zeros(1,80000);)a (1)=1;b (1)=0;for k=2:80000a(k)=0.99803*a(k-1)-0.06279*b(k-1);b(k)=0.06279*a(k-1)+0.99803*b(k-1);end编写高效编程代码（掌握）在程序编制时，如果变量的大小随着计算会扩展，预先给变量设定好尺寸会使得代码更加高效在程序编制时，如果变量的大小随着计算会扩展，预先给变量设定好尺寸会使得代码更加高效clear,tica=zeros(1,10000);%Preallocatecount=0;for k=1:10000v=exp(rand (1)*rand (1);if v0.5%Conditionally addto arraycount=count+1;a(count)=v;endenda=a(1:count);%Trim theresulttoc在实际MATLAB编程中，采用循环语句会降低其执行速度，用其内在函数批量处理大矩阵可以有效提高计算速度编程中，采用循环语句会降低其执行速度，用其内在函数批量处理大矩阵可以有效提高计算速度矢量化运算 （1）知已知y=1+1/3+1/5+1/2n-1，当n=100时，求y的值？n=100;i=1:2:2*n-1;y=sum(1./i);矢量化运算 （2）function d=minDistance(x,y,z)%Find themin distancebetween aset ofpoints andthe originnPoints=length(x);d=zeros(nPoints,1);%Preallocater kfor k=1:nPoints%Compute distancefor everypointd(k)=sqrt(x(k)?2+y(k)?2+z(k)?2);endd=min(d);%Get theminimum distancefunctiond=minDistance(x,y,z)%Find themin distancebetween aset ofpoints andthe