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

Matlab在电子信息类中的应用实验专业 ： 成绩： 姓名： 学号： 实验三: M函数和M文件的编写与应用实验目的: 1 熟悉M函数和M文件的编写与应用2 实验内容:（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. %myprime3. 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);b=input(b);c=input(c);disp(,num2str(a),*x2+,num2str(b),*x+,num2str(c)if a=0 x=-b/c disp(x=,num2str(x)elseif b2-4*a*c0 x1=(-b+sqrt(b2-4*a*c)/(2*a); x2=(-b-sqrt(b2-4*a*c)/(2*a); disp(x1=,num2str(x1),;x2=,num2str(x2)elseif b2-4*a*c=0 x=-b/(2*a); disp(x1=x2=,num2str(x)else disp()end直接由公式先判断b2-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 s0j=1;while j=sy(t+j)=x(m*j);j=j+1;endx(m:m:m*s)=;if 3*spx=x(m-1)*s+1:p-s,1:(m-1)*s); endt=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),y3endy1 =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.1302Columns 16 through 300.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 450.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 Columns 46 through 60 0.3487 0.3542 0.3596 0.3648 0.3700 0.3750 0.3800 0.3848 0.3896 0.3942 0.3987 0.4032 0.4076 0.4118 0.4160 Columns 61 through 75 0.4200 0.4240 0.4278 0.4316 0.4352 0.4387 0.4422 0.4455 0.4488 0.4520 0.4550 0.4579 0.4608 0.4636 0.4662Columns 76 through 900.4688 0.4712 0.4735 0.4758 0.4779 0.4800 0.4819 0.4838 0.4855 0.4872 0.4888 0.4902 0.4916 0.4928 0.4939Columns 91 through 1050.4950 0.4960 0.4968 0.4975 0.4982 0.4987 0.4992 0.4996 0.4998 0.5000 0.5000 0.4999 0.4998 0.4995 0.4992Columns 106 through 120 0.4988 0.4982 0.4975 0.4968 0.4960 0.4950 0.4940 0.4928 0.4915 0.4902 0.4888 0.4872 0.4855 0.4838 0.4820Columns 121 through 135 0.4800 0.4780 0.4758 0.4735 0.4712 0.4688 0.4662 0.4635 0.4608 0.4580 0.4550 0.4520 0.4488 0.4455 0.4422Columns 136 through 150 0.4387 0.4352 0.4315 0.4278 0.4240 0.4200 0.4160 0.4118 0.4076 0.4032 0.3987 0.3942 0.3896 0.3848 0.3800 Columns 151 through 1650.3750 0.3700 0.3648 0.3596 0.3542 0.3487 0.3432 0.3376 0.3318 0.3260 0.3200 0.3140 0.3078 0.3016 0.2952Columns 166 through 180 0.2888 0.2822 0.2755 0.2688 0.2620 0.2550 0.2479 0.2408 0.2336 0.2262 0.2188 0.2112 0.2035 0.1958 0.1879Columns 181 through 195 0.1800 0.1719 0.1638 0.1555 0.1472 0.1387 0.1302 0.1216 0.1128 0.1040 0.0950 0.0859 0.0768 0.0676 0.0582Columns 196 through 2010.0488 0.0392 0.0295 0.0198 0.0099 0y2 =Columns 1 through 15 0 0.0199 0.0396 0.0591 0.0784 0.0975 0.1164 0.1351 0.1536 0.1719 0.1900 0.2079 0.2256 0.2431 0.2604 Columns 16 through 300.2775 0.2944 0.3111 0.3276 0.3439 0.3600 0.3759 0.3916 0.4071 0.4224 0.4375 0.4524 0.4671 0.4816 0.4959Columns 31 through 450.5100 0.5239 0.5376 0.5511 0.5644 0.5775 0.5904 0.6031 0.6156 0.6279 0.6400 0.6519 0.6636 0.6751 0.6864Columns 46 through 600.6975 0.7084 0.7191 0.7296 0.7399 0.7500 0.7599 0.7696 0.7791 0.7884 0.7975 0.8064 0.8151 0.8236 0.8319Columns 61 through 750.8400 0.8479 0.8556 0.8631 0.8704 0.8775 0.8844 0.8911 0.8976 0.9039 0.9100 0.9159 0.9216 0.9271 0.9324 Columns 76 through 900.9375 0.9424 0.9471 0.9516 0.9559 0.9600 0.9639 0.9676 0.9711 0.9744 0.9775 0.9804 0.9831 0.9856 0.9879Columns 91 through 1050.9900 0.9919 0.9936 0.9951 0.9964 0.9975 0.9984 0.9991 0.9996 0.9999 1.0000 0.9999 0.9996 0.9991 0.9984 Columns 106 through 120 0.9975 0.9964 0.9951 0.9936 0.9919 0.9900 0.9879 0.9856 0.9831 0.9804 0.9775 0.9744 0.9711 0.9676 0.9639Columns 121 through 1350.9600 0.9559 0.9516 0.9471 0.9424 0.9375 0.9324 0.9271 0.9216 0.9159 0.9100 0.9039 0.8976 0.8911 0.8844Columns 136 through 150 0.8775 0.8704 0.8631 0.8556 0.8479 0.8400 0.8319 0.8236 0.8151 0.8064 0.7975 0.7884 0.7791 0.7696 0.7599Columns 151 through 165 0.7500 0.7399 0.7296 0.7191 0.7084 0.6975 0.6864 0.6751 0.6636 0.6519 0.6400 0.6279 0.6156 0.6031 0.5904 Columns 166 through 180 0.5775 0.5644 0.5511 0.5376 0.5239 0.5100 0.4959 0.4816 0.4671 0.4524 0.4375 0.4224 0.4071 0.3916 0.3759 Columns 181 through 1950.3600 0.3439 0.3276 0.3111 0.2944 0.2775 0.2604 0.2431 0.2256 0.2079 0.1900 0.1719 0.1536 0.1351 0.1164 Columns 196 through 201 0.0975 0.0784 0.0591 0.0396 0.0199 0y3 = Columns 1 through 15 0 0.0099 0.0198 0.0295 0.0392 0.0488 0.0582 0.0676 0.0768 0.0859 0.0950 0.1040 0.1128 0.1216 0.1302Columns 16 through 300.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 450.2550 0.2620 0.2688 0.2755 0.2822 0.2888 0.2952 0.3016 0.3078 0.3140 0.3200 0.3260 0.3318 0.3376 0.3432Columns 46 through 600.3487 0.3542 0.3596 0.3648 0.3700 0.3750 0.3800 0.3848 0.3896 0.3942 0.3987 0.4032 0.4076 0.4118 0.4160Columns 61 through 75 0.4200 0.4240 0.4278 0.4315 0.4352 0.4387 0.4422 0.4455 0.4488 0.4520 0.4550 0.4580 0.4608 0.4635 0.4662Columns 76 through 900.4688 0.4712 0.4735 0.4758 0.4780 0.4800 0.4820 0.4838 0.4855 0.4872 0.4888 0.4902 0.4915 0.4928 0.4940Columns 91 through 1050.4950 0.4960 0.4968 0.4975 0.4982 0.4988 0.4992 0.4995 0.4998 0.4999 0.5000 0.5000 0.4998 0.4996 0.4992Columns 106 through 120 0.4987 0.4982 0.4975 0.4968 0.4960 0.4950 0.4939 0.4928 0.4916 0.4902 0.4888 0.4872 0.4855 0.4838 0.4819 Columns 121 through 135 0.4800 0.4779 0.4758 0.4735 0.4712 0.4688 0.4662 0.4636 0.4608 0.4579 0.4550 0.4520 0.4488 0.4455 0.4422 Columns 136 through 150 0.4387 0.4352 0.4316 0.4278 0.4240 0.4200 0.4160 0.4118 0.4076 0.4032 0.3987 0.3942 0.3896 0.3848 0.3800Columns 151 through 1650.3750 0.3700 0.3648 0.3596 0.3542 0.3487 0.3432 0.3376 0.3318 0.3260 0.3200 0.3140 0.3078 0.3016 0.2952Columns 166 thr