min=f(i); end end sum max min p n z
5、qmss.m sum=0;n=0; for i=2:50
m=i*(i+1)-1; k=sqrt(m); flag=1; for j=2:k
if rem(m,j)==0 flag=0; end end
if flag~=0 n=n+1;
sum=sum+m; end end n sum
===================================================================================
MATLAB教程实验报告
实验项目名称 实验五 函数文件
学生姓名 专业班级 学号 实验成绩 日期
五. 实验目的和要求
1、理解函数文件的概念。
2、掌握定义和调用MATLAB函数的方法。
3、填写实验报告,实验报告文件取名为report5.doc。
4、于邮件附件形式将实验报告文件report5.doc 发到邮箱matlab0909@sina.com,邮件主题为班级学号姓名,如:09数教1班15号张三。
二、实验内容
1、定义一个函数文件,求给定复数的指数、对数、正弦和余弦,并在命令文件中调用该函数文件。 2、一个物理系统可用下列方程组来表示:
?m1cos??msin??1?0??0?m10m20?sin?cos??sin??cos?0??a1??0??a??mg?0???2???1? 0??N1??0??????1??N2??m2g?从键盘输入m1、m2和?的值,求a1、a2、N1和N2的值。其中g取9.8,输入?时以角度为单位。 要求:定义一个求解线性方程组AX=B的函数文件,然后在命令文件中调用该函数文件。
3、一个自然数是素数,且它的数字位置经过任意对换后仍是素数,则称是绝对素数。例如13是绝对素数。试求所有两位绝对素数。 要求:定义一个判断素数的函数文件。 4、设f(x)?11?,编写一个MATLAB函数文件fx.m,使得调用f(x)时,x24(x?2)?0.1(x?3)?0.01可用矩阵代人,得出的f(x)为同阶矩阵。 5、已知y?f(40)
f(30)?f(20)(1)当f(n)?n?10ln(n2?5)时,求y的值。
???n?(n?1)(2)当f(n)?1?2?2?3?3?4??时,求y的值。
三. 实验结果与分析
(包括运行结果截图、结果分析等) 四. 心得体会
(记录实验感受、上机过程中遇到的困难及解决办法、遗留的问题、意见和建议等。) 【附录----源程序】 1、clear all
a=input('请输入一个复数:'); [e,l,s,c]=fushu(a); 2、clear all
m1=input('请输入m1的值:'); m2=input('请输入m2的值:'); m3=input('请输入θ的值:'); J=jiefangcheng(m1,m2,m3); 3、clear all for n=10:99
a=sushu(n); end
4、clear all
y=input('请输入一个数或矩阵:'); disp('输入的数或矩阵x是:') disp(y) L=fx(y);
5、clear all for n=20:10:40 if n==20;
[f1,f2]=f(n); a=f1; b=f2;
elseif n==30;
[f1,f2]=f(n); c=f1; d=f2; else
[f1,f2]=f(n); e=f1; f=f2; end end
y1=e/(a+c); y2=f/(b+d);
disp(['(1) y=',num2str(y1)]) disp(['(2) y=',num2str(y2)])