实验三-MATLAB-程序设计

《实验三-MATLAB-程序设计》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、实验三MATLAB程序设计题目:1,编写M函数实现求一个数是否为素数，再编写一主程序(脚本文件)，要求通过键盘输入一个整数，然后调用判断素数函数，从而确定它是否素数。2,编写M函数统计一数值屮零的个数，然后编写脚本文件，实现统计从1-2007中零的总个数。3,编写程序计算xe[-3,3],字长0.01：-x2-4*x-3)/2-3

2、：1.%脚本文件x二input('请输入一个整数x:');ifmyprime(x)disp(，您输入的整数x是一个素数。’)elsedispC您输入的数x不是一个素数。’)end运行结果：请输入一个整数x:9您输入的数x不是一个素数。2.%脚本文件s二0;fori二1:2007s=s+my_total_0(i);ends%myprime函数functionyLDiypriiiicd)y=l;fori=2:fix(sqrt(x))ifmod(x,i)=0y=0;endend%my_total_0函数functiony二mytotalO(x)y=0;n=fix(loglO

3、(x))+1;fori=l:ny二y+、inod(x,10);x=fix(x/10);end-0.5242-0.6428-0.58390.08710.11060.9630-o.2309>>a=chap4_4(4)-0.8436-0.4764选帑楚阵形式运行结果:5063.functiony二chap4_3(step,1incstyle)y=0；xl二-3:step:-l-step;x2二T:step:1-step;x3=l:step:3;yl=(-xl.J-4*xl-3)/2;y2=-x2."2+1;y3=(-x3.八2+4*x3-3)/2;x二[xl,x2,x3];

4、yy二[yl,y2,y3];plot(x,yy,lincstylc)运行结果：>>y=chap4_3(0.01,'ro,)4.functiona=chap4_4(n)ch二menu('选择矩阵形式','[TO,10]均匀分布','[-5,5]均匀分布','[-1,1]均匀分布','正态分布')；switchchcase1,a=10-20*rand(n);case2,a=5-10^rand(n);case3,a=l~2*rnncl(n);case4,a=randn(n);end运行结果:»a=chap4_4(l)a二-9.0026»a=chap4_4(2)a二2.688

5、60.1402-1.0684-3.9130>>a=chap4_4(3)-0.4326-1.14650.3273-1.66561.19090.17460.12531.1892-0.18670.2877-0.03760.7258functiony=chap4_5(1)n=length(I);fori二l:nifT(i)<=10y(i)二0.l*I(i);elseifl(i)<=20y(i)=0.05*(I(i)-10)+1;-0.58832.1832-0.13640.1139elseifI(i)<=40y(i)=0.02*(I(i)-20)+1.5;elsey(i)二0.

6、01*(I(i)-40)+1.9;endencly=[l:n;I;y];运行结果：»y=chap4_5(25)y=1.000025.00001.60006.%求出各项fori=l:na(i)二x(i+l)/x(i);end%求出前n项和y=sum(a);functiony二chap4_6(n)x⑴二1;x(2)二2;%求出费波纳契序列fori=3:n+1x(i)=x(i-2)+x(i-l);end运行结果：»y=chap4_6(15)y=24.57017.functiony二chap4_7(m,n)x=l:n;fori二l:nN=n+l-i;t=mod(in,N);i

7、ft二二0t二N;endy(i)二x(t);x=[x(t+l:N),x(l:t-1)];end运行结果：»y=chap4_7(7,20)Columns1through111751541682071419Columns12through20131110121961823