matlab程序设计(数学班).docVIP

Matlab程序设计 例一： 例二： 例三：显示函数y=e-xsin(2pix)的图形 例四：设计一个计算和的程序 用下列方法得出的结果有何不同？ 从左到右各项相加 从右到左各项相加 从左到右各个正项和负项分别相加；从右到左各个正项和负项分别相加 例四：求自然对数e的近似值 当余项小于精度要求即可以。 例五、 例六：求正整数a,b之间的所有完全数。完全数 例七、歌德巴赫猜想 DJS--130是中小规模集成电路的计算机，运算速度每秒50万次，机器字长16位，内存容量只有32K字节，并且是磁芯内存，也没有真正意义上的操作系统。 程序设计第二讲 例一：仔细研究下面例子，具有广泛的指导与模仿意义，要求能逐条解释清楚。 本例中使用了特定变量nargin表示输入变量数目，nargout表示输出变量数目，它们是matlab中的永久变量，通常根据nargin和nargout的数目不同而调用不同的程序段，从而体现它的智能作用。 function [out1,out2]=sincos(x) if nargin==0 x=0:0.01:2*pi; end; y=sin(x).*cos(x); if nargout==2 out1=x; out2=y; else out1=y; end 例二：细读下列不动点迭代程序，仔细体会输出格式 clear all; x(1)=1.5;x(2)=((x(1))^2+1)^(1/3); i=1;k=[1,2]; while abs(x(i+1)-x(i))=0.5e-10; i=i+1;k=[k,i]; x(i+1)=((x(i))^2+1)^(1/3); end; sprintf(%20.17f\n,x); sprintf(i=%10d\n,i); sprintf(i=%5d , x=%20.15f\n,[k;x]) %sprintf(i=%5d , x=%20.15f\n,k,x) %sprintf(i=%5d , x=%20.15f\n,[k,x]) 例三：把上面程序改造成函数（输出参数可以变） function [xstar,cs]=budongdiang2(x0,wc) x(1)=x0;x(2)=((x(1))^2+1)^(1/3); i=1;k=[1,2]; while abs(x(i+1)-x(i))=wc; i=i+1;k=[k,i]; x(i+1)=((x(i))^2+1)^(1/3); end; if nargout==1 xstar=x(i); end; if nargout==2 xstar=x(i); cs=i; end; 例四：把上面程序改造成函数（输入/出参数可以变） function [xstar,cs]=budongdiang2(x0,wc) if nargin==1 wc=1e-15; end; if nargin==0 x0=0;wc=1e-15; end; x(1)=x0;x(2)=((x(1))^2+1)^(1/3); i=1;k=[1,2]; while abs(x(i+1)-x(i))=wc; i=i+1;k=[k,i]; x(i+1)=((x(i))^2+1)^(1/3); end; if nargout==1 xstar=x(i); end; if nargout==2 xstar=x(i); cs=i; end; 例五、歌德巴赫猜想 function [flag]=isprime(x) flag=1; for i=2:x-1 if rem(x,i)==0 flag=0; break; end; end function k=gdbhcx(n) k=[]; for i=2:(n/2) if isprime(i)isprime(n-i) k=[k,sprintf(%d=%d+%d\n,n,i,n-i)]; end end k=[k,sprintf(---------------\n)]; function s=cx(m,n) s=[]; for i=m:2:n s=[s,gdbhcx(i)]; end 例：三角形面积及三个角 function [out1,out2,out3,out4]=triangle(a,b,c) if(a+b=c)|(a+c=b)|(b+c=a) sprintf(2??ü113éèy??D?); else p=(a+b+c)/2; s=sqrt(p*(p-a)*(