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

实验三M函数和M文件的编写与应用以作三道题 Matlab在电子信息类中的应用实验 专业：成绩： 姓名：学号： 实验三:M函数和M文件的编写与应用 实验目的： 熟悉M函数和M文件的编写与应用 2?实验内容： （1） 编写M函数实现：求一个数是否为素数，再编写一主程序（脚 本文件），要求通 过键盘输入一个整数，然后判断其是否为素数。 （2） 编写求解方程ax+bx+c二0的根的函数，这里应根据b-4ac的不 同取值分别处理， 并输入几组典型值加以检验。 （3） 有n个人围成一圈，按序列编号。从第1个人开始报数，数到 m时该人退出， 并且下一个从1重新开始报数，求出出圈人的顺序（ngt;m,例如n=20, m二刀。 （4） 编写M函数计算（x?[?3,3],步长0.01） ?（?x2?4x?3）/2?3?x??l??y???x2?l?l?x?lo ?(?x2?4x?3)/21?x?3??22 ⑸设 fx=ll+写一个MATLAB函数程序，使得调用此函数时，x可以用矩阵 代入，得出f(x)为同阶矩阵。实验结果： myprime 函数 %myprimeo 一 ey 3. function y二myprime(x) y=i; fori=2:fix(sqrt(x)) if mod(x,i)==0 y=0; end end 11. Matlab在电子信息类中的应用实验 脚本函数： 实验结果： 编译函数直接开放取余，看是否整出。直接输入一个数会直接判断是 否是素数。 2. a=input(#39;??e?e?a u ?? u ￡o#39;); b=input(#39;??e?e?b u ?? u ￡ o#39;); c=input(#39;??e?e?c u ?? u ￡ o#39;); disp([#39; ? ?3i ￡ o#39;,num2str(a),#39;*xA2+#39;,num2str(b),#39;*x+#39;,num2str( c)]) if a==0 x=-b/c disp([#39;?u?a ￡ ox=#39;,num2str(x)]) elseif bA2-4*a*cgt;O xl=(-b+sqrt(bA2-4*a*c))/(2*a); x2=(-b-sqrt(bA2-4*a*c))/(2*a); disp([#39;?u?a ￡ oxl=#39;?num2str(xl),#39;;x2=#39;,num2str(x2)]) elseif bA2-4*a*c==0 x=-b/(2*a); disp([#39;?u?a ￡ oxl=x2=#39;/num2str(x)]) else disp(#39;?T?a #39;) end Matlab在电子信息类中的应用实验 直接由公式先判断bA2-4ac然后根据不同的情况去取算法或者输出值。 function y二baoshu(m,n) x=l:n; y=zeros(l, m+1); t=0; p=size(x/2); s=fix(p/m); while sgt;O j=l； while jlt;=s y(t+j)=x(m*j); j=j+l； end x([m:m:m*s])=[]; 讦 3*slt;p x=x([(m-l)*s+l:p-s,l*(m-l)*s]); end t=t+s; p=size(x/2); s=fix(p/m); end Matlab在电子信息类中的应用实验 循环嵌套，取出。 function [y]二 fun3() xl=-3:0.01:-l; yl=(-xl.A2-4.*xl-3)./2 x2=-l:0.01:l; y2=-x2.A2+l x3=l:0.01:3; y3=(-x3.A2+4.*x3-3)./2 y=[yl(:/[l:200])/y2(:/[l:200])/y3] end yi = Colum ns 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 0.2822 0.2887 0.2952 0.3016 0.3078 0.3140 0.3200 0.3260 0.3318 0.3376 0.3432 Colu