实验三M函数和M文件的编写与应用以作三道题详解.docxVIP

专业 ： 成绩： 姓名： 学号： 实验三: M函数和M文件的编写与应用 实验目的: ⒈ 熟悉M函数和M文件的编写与应用 ⒉ 实验内容: （1） 编写M函数实现：求一个数是否为素数，再编写一主程序(脚本文件)，要求通过键盘输入一个整数，然后判断其是否为素数。 （2） 编写求解方程ax2+bx+c=0的根的函数，这里应根据b2-4ac的不同取值分别处理，并输入几组典型值加以检验。 （3） 有n个人围成一圈，按序列编号。从第1个人开始报数，数到m时该人退出，并且下一个从1重新开始报数，求出出圈人的顺序(nm，例如n=20，m=7)。 （4） 编写M函数计算（ ，步长0.01） 。 (5)设fx=1(x-2)2+0.1+1(x-3)3+0.01写一个MATLAB函数程序，使得调用此函数时，x可以用矩阵代入，得出f(x)为同阶矩阵。 实验结果： 1. myprime函数 2. %myprimeoˉêy 3. 4. function y=myprime(x) 5. y=1; 6. for i=2:fix(sqrt(x)) 7. if mod(x,i)==0 8. y=0; 9. end 10. end 11. 脚本函数： 实验结果： 编译函数直接开放取余，看是否整出。直接输入一个数会直接判断是否是素数。 2. a=input(??ê?è?aμ??μ￡o); b=input(??ê?è?bμ??μ￡o); c=input(??ê?è?cμ??μ￡o); disp([·?3ì￡o,num2str(a),*x^2+,num2str(b),*x+,num2str(c)]) if a==0 x=-b/c disp([?ù?a￡ox=,num2str(x)]) elseif b^2-4*a*c0 x1=(-b+sqrt(b^2-4*a*c))/(2*a); x2=(-b-sqrt(b^2-4*a*c))/(2*a); disp([?ù?a￡ox1=,num2str(x1),;x2=,num2str(x2)]) elseif b^2-4*a*c==0 x=-b/(2*a); disp([?ù?a￡ox1=x2=,num2str(x)]) else disp(?T?a) end 直接由公式先判断b^2-4ac然后根据不同的情况去取算法或者输出值。 3. function y=baoshu(m,n) x=1:n; y=zeros(1,n-m+1); t=0; p=size(x,2); s=fix(p/m); while s0 j=1; while j=s y(t+j)=x(m*j); j=j+1; end x([m:m:m*s])=[]; if 3*sp x=x([(m-1)*s+1:p-s,1:(m-1)*s]); end t=t+s; p=size(x,2); s=fix(p/m); end 循环嵌套，取出。 4. function [y] = fun3() x1=-3:0.01:-1; y1=(-x1.^2-4.*x1-3)./2 x2=-1:0.01:1; y2=-x2.^2+1 x3=1:0.01:3; y3=(-x3.^2+4.*x3-3)./2 y=[y1(:,[1:200]),y2(:,[1:200]),y3] end y1 = Columns 1 through 15 0 0.0099 0.0198 0.0295 0.0392 0.0488 0.0582 0.0675 0.0768 0.0859 0.0950 0.1040 0.1128 0.1216 0.1302 Columns 16 through 30 0.1387 0.1472 0.1555 0.1638 0.1719 0.1800 0.1879 0.1958 0.2035 0.2112 0.2188 0.2262 0.2336 0.2408 0.2479 Columns 31 through 45 0.2550 0.2620 0.2688 0.2755