matlab程序设计(数学班)

1、Matlab程序设计例一：例二：例三：显示函数y=e-xsin(2pix)的图形例四：设计一个计算和的程序用下列方法得出的结果有何不同？从左到右各项相加 从右到左各项相加从左到右各个正项和负项分别相加；从右到左各个正项和负项分别相加例四：求自然对数e的近似值当余项小于精度要求即可以。例五、例六：求正整数a,b之间的所有完全数。完全数例七、歌德巴赫猜想 DJS-130是中小规模集成电路的计算机，运算速度每秒50万次，机器字长16位，内存容量只有32K字节，并且是磁芯内存，也没有真正意义上的操作系统。程序设计第二讲例一：仔细研究下面例子，具有广泛的指导与模仿意义，要求能逐条解释清楚。本例中使用了特

2、定变量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)

3、-x(i)=0.5e-10; i=i+1;k=k,i; x(i+1)=(x(i)2+1)(1/3);end;sprintf(%20.17fn,x);sprintf(i=%10dn,i);sprintf(i=%5d , x=%20.15fn,k;x)%sprintf(i=%5d , x=%20.15fn,k,x)%sprintf(i=%5d , x=%20.15fn,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

4、)-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

5、=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;endfunction k=gdbhcx(n)k=;for i=2:(n/2) if isprime(i)&isprime(n-i) k=k,sprintf(%d=%d+%dn,n,i,n-i); endendk=k,sprintf(-n);

6、 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();else p=(a+b+c)/2; s=sqrt(p*(p-a)*(p-b)*(p-c); ja=acosd(b*b+c*c-a*a)/(2*b*c); jb=acosd(a*a+c*c-b*b)/(2*a*c); jc=acosd(b*b+a*a-c*c)/(2*b*a);end;if nargout=0 send;if nargout=1 out1=s;end;if nargout=2 out1=s; out2=ja;