- 1、本文档共16页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实验三MATLAB程序设计
实验三MATLAB程序设计
PAGE
实验三MATLAB程序设计
2014秋2012级《MATLAB程序设计》实验报告
班级:软件C121 姓名:冯杨腾 学号:125692
实验三 MATLAB程序设计
一、实验目的
1、掌握建立和执行M文件的方法。
2、掌握选择结构、多分支选择结构程序的编程方法。
3、掌握多种循环结构程序的编程方法。
4、掌握定义函数文件和调用函数文件的方法。
二、实验内容
1、MATLAB M文件的创建与使用。
建立自己工作目录,如D:\MATLAB。在File菜单中,单击Set Path菜单项,在弹出对话框中建立目标文件夹,:单击Add Folder按钮,将自己的工作目录D:\MATLAB添加到搜索目录之中,按Save按钮保存,如图1。
图1 工作目录设置
在M文件编辑器中编制命令行文件。在File菜单中,单击New—M-File,弹出M文件编辑器,如图2:
图2 编辑器
在编辑器中输入如下命令程序:
%求小于2000且为2的整数次幂的正整数。
f(1)=2;
k=1;
while f(k)1000
f(k+1)=f(k)*2;
k=k+1;
end
f,k
调试后,将其以形式保存在D:\MATLAB中。
命令行程序运行。进入命令窗口,键入test1,观察运行结果:
f =
Columns 1 through 9
2 4 8 16 32 64 128 256 512
Column 10
1024
k =
10
按上述步骤,在在编辑器中输入如下函数文件:
function f=tt(n)
%求小于任何正整数且为2的整数次幂的正整数。
%c=n(n)
%n可取任意正整数。
%2004年8月。
f(1)=2;
k=1;
while f(k)ceil(n/2) %ceil函数取最近的正整数
f(k+1)=f(k)*2;
k=k+1;
end
f
调试后,将其以形式保存在D:\MATLAB中。
进入命令窗口,键入test2(9),观察程序运行结果。
test2(9)
f =
2 4 8
ans =
2 4 8
2、编写一个函数M文件,计算下面函数的值,给出标量x的值,调用该函数后,返回y的值。
选择一些数据测试编写的函数。
源程序 :
function [y]=myfun(x)
if x=0
y=sin(x);
elseif x3
y=-x+6;
else
y=x;
end
运行结果:
test(0)
ans =
0
test(2)
ans =
2
test(5)
ans =
1
3、编写一个函数M文件,返回向量x中元素的和、平均值、最大值、最小值、均方根值。其中方均根值(Root Mean Square)的计算公式为:
并用下面数据测试你写的函数:
x=sin(0::6*pi)
x=rand(1,200),得到的x为200个(0,1)之间均匀分布的随机数。
源程序:
function [sum_x,m_x,max_x,min_x,rsm_x] = myfun7(x)
sum_x=0;
a=x;
for k=a;
sum_x=sum_x+k;
n=length(x);
for i=1:n
m_x=sum_x/n;
end
end
min_x=x(1);
for t=1:n
if min_xx(t)
min_x=x(t);
end
end
max_x=x(1);
for t=1:n
if max_xx(t)
max_x=x(t);
end
end
sum_x1=0;
for t=1:n
sum_x1=sum_x1+x(t).^2;
rms_x=sqrt(sum_x1/n);
end
sum_x
m_x
max_x
min_x
rms_x
运行结果:
(1)
x=sin(0::6*pi);
myfun7(x)
sum_x =
m_x =
max_x =
min_x =
rms_x
文档评论(0)