matlab程序設(shè)計(數(shù)學班)

1、Matlab程序設(shè)計例一：例二：例三：顯示函數(shù)y=e-xsin(2pix)的圖形例四：設(shè)計一個計算和的程序用下列方法得出的結(jié)果有何不同？從左到右各項相加 從右到左各項相加從左到右各個正項和負項分別相加；從右到左各個正項和負項分別相加例四：求自然對數(shù)e的近似值當余項小于精度要求即可以。例五、例六：求正整數(shù)a,b之間的所有完全數(shù)。完全數(shù)例七、歌德巴赫猜想 DJS-130是中小規(guī)模集成電路的計算機，運算速度每秒50萬次，機器字長16位，內(nèi)存容量只有32K字節(jié)，并且是磁芯內(nèi)存，也沒有真正意義上的操作系統(tǒng)。程序設(shè)計第二講例一：仔細研究下面例子，具有廣泛的指導與模仿意義，要求能逐條解釋清楚。本例中使用了特

2、定變量nargin表示輸入變量數(shù)目，nargout表示輸出變量數(shù)目，它們是matlab中的永久變量，通常根據(jù)nargin和nargout的數(shù)目不同而調(diào)用不同的程序段，從而體現(xiàn)它的智能作用。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)例三：把上面程序改造成函數(shù)（輸出參數(shù)可以變）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; 例四：把上面程序改造成函數(shù)（輸入/出參數(shù)可以變）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;