您所在位置： &
&nbsp&&nbsp&nbsp&&nbsp
MATLAB编程(第二版)-菜鸟入门教材.pdf 339页
本文档一共被下载：
次 ,您可全文免费在线阅读后下载本文档。
下载提示
1.本站不保证该用户上传的文档完整性，不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值，立即自动返金币，充值渠道很便利
MATLAB编程(第二版)-菜鸟入门教材.pdf
你可能关注的文档：
··········
··········
Stephen J. Chapman 著
邢树军 郑碧波译
国外高校电子信息类优秀教材（英语影印版）
MATLAB 编程
（第二版）
MATLAB Programming for Engineers
（Second Edition）
Stephen J. Chapman
邢树军 郑碧波 译
本文为 Stephen
《MATLAB 编程（第二版）》英文影印版的中文译本，这本
书对初学者很好的入门教材。之前从网上下载到邢树军兄弟译的 pdf 版本，邢的这个版本用
起来有以下几方面不便：
1.文件较大。
2.打开之后总是提示播放某首歌，很烦。
3.源代码和图形采用了大量的图片，不方便学习者直接复制到 MATLAB 中进行调试学
习。即使是文本的源代码，也经常出错，主要是每一行语句末尾的“；”为中文编码的分号。
4.个别字翻译错误（sorry，我的这个版本也没有完全消除），但我有发现的已经更正。
5.第十章翻译不完全。本来全书只有十章，在第十章还没有翻译结束时插进了第十一章。
6.附录 A 和 B 没有（翻译）提供。其中附录 B 为练习答案，其实这对初学者来说，是
极其重要的。
于是就想着改进该版本的。本版本在邢译本的基础上做如下修改：
1.转换成 doc 格式，重新排版（其实也没怎么排版）。
2.更换封面，去掉歌曲。
3.所有源代码测试正确。
4.用了较多新的彩色截图代替原来的。
5.提供了附录 A 和附录 B 。
6.删除了原来 pdf 版本的第 11 章（原著没有）。
本版本有如下缺点
1.由于本人截图用MATLAB 是 7.0 和 6.5 版本，在 winxp 中运行。一些截图与原英文版
本不一样，原因一是与原英文版本 MATLAB
6.0 不同，原因二是因于产生图形的数据无法
找到相同的，好在基础知识方面的学习差别不大。
2.第十章前几节仍为树军兄弟翻译，后面部分由我翻译，但由于下载不到 6.0 版本，只
好用 6.5 版本，但使用 6.5 版本会出现 GUI 编程代码不一样的情况，从忠于原著角度考虑感
觉不爽，后来想想，既然网上都找不到 6.0 版本，那么使用 6.5 版本岂不更方便读者？！
3.附录 B 与前面的专有名词翻译可能会不一致。例如 cell
array 一词，邢把它译为“单
元阵列”，而我把它称为“胞数组”。我是一个初学者，不知权威的书是如何翻译的。不管如
何，我已经做了相应说明，希望不会引起读者太大的不便。
本人英语也忘记得七七八八了。在翻译过程中肯定有不少不恰当的地方，如果你发现，
请通知我，谢谢！
本书的绝大部分是邢树军翻译的，感谢邢树军兄弟的无私奉献，我仅仅是做些小小修补。
其实翻译之后才发现，翻译过程是一项很繁琐很枯燥的事情！
有这么一个故事：
在非洲草原上，有一种不起眼的动物叫吸血蝙蝠。它身体极小，却是野马的天敌。这种
蝙蝠靠吸动物的血生存，它在攻击野马时，常
正在加载中，请稍后...只需一步，快速开始
扫一扫，访问微社区
请完成以下验证码
查看: 5298|回复: 5|关注: 0
菜鸟求助：用fsolve解超越方程（怎么添加循环？）
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
菜鸟求助：fsolve解超越方程M文件：
function output = myfun(x,k)
& &output(1) =1-2*cos(pi*x(1)/180)+2*cos(pi*x(2)/180)-k/100
&&tput(2) =1-2*cos(pi*5*x(1)/180)+2*cos(pi*5*x(2)/180);
再用fsolve把x1,x2,x3解出来，并且把x1，x2，x3与k的函数图输出来。
现在的问题是，不知道把k从1到100循环怎么添加进去，并且把图输出来。
菜鸟提问，请大家帮忙详细解答啊，谢谢！
<h1 style="color:#6 麦片财富积分
关注者： 3
for k=1:100
x=fsolve(f,x0)
<h1 style="color:# 麦片财富积分
for k=0:1:100
& && &fsolve(@myfun,[30,30],[],k);
& && &plot(x(1),k),hold on
这是我写的，貌似不对，帮我看看....谢谢
for k=1:100
x=fsolve(f,x0)
end&&好像也不不对啊。
<h1 style="color:#0 麦片财富积分
关注者： 49
这个可以利用嵌套函数结构来解决：
function PlotX_k
function output = myfun(x)
& &output(1) =1-2*cos(pi*x(1)/180)+2*cos(pi*x(2)/180)-k/100;
&&output(2) =1-2*cos(pi*5*x(1)/180)+2*cos(pi*5*x(2)/180);
end
X = zeros(95,2);
for k = 1:95& &
& & X(k,:) = fsolve(@myfun,[0.01,12]);
end
& & plot(1:95,X(:,1));
& & xlabel('k');ylabel('x1')
end
复制代码此外，由于楼主的目标函数形式比较简单，也可以利用匿名函数，这样的话，就不用单独建立m文件了。我的书里关于嵌套函数和匿名函数有较多篇幅的介绍，嵌套函数和匿名函数都是自MATLAB7.0以后推出的新的结构，Mathworks公司一直大力提倡使用。
fun = @(k)@(x) [1-2*cos(pi*x(1)/180)+2*cos(pi*x(2)/180)-k/100;1-2*cos(pi*5*x(1)/180)+2*cos(pi*5*x(2)/180);]
X = zeros(95,2);
for k = 1:95
X(k,:) = fsolve(fun(k),[30,30]);
end
plot(1:100,X(:,1));
xlabel('k');ylabel('x1')
plot(1:95,X(:,1));
xlabel('k');ylabel('x1')
复制代码还需要说明一点的是，当k接近100的时候，方程变得病态起来，譬如k = 100的时候，由原方程组第一个方程可得
cos(x(2)*pi/180) = cos(x(1)*pi/180);
而第二个方程又要满足：cos(5*x(2)*pi/180) = cos(5*x(1)*pi/180)-1/2;
我们可以用MATLAB符号运算验证下当cos(a) = cos(b)时，cos(5*a) = cos(5*b),所以方程是不成立的。
&& syms a
&& expand(cos(5*a))
ans =
cos(a)^5 - 10*cos(a)^3*sin(a)^2 + 5*cos(a)*sin(a)^4
复制代码
[ 本帖最后由 吴鹏 于
22:35 编辑 ]
<h1 style="color:# 麦片财富积分
非常感谢你！！！
<h1 style="color:# 麦片财富积分
这个可以利用嵌套函数结构来解决：此外，由于楼主的目标函数形式比较简单，也可以利用匿名函数，这样的话， ...
您好，您这里的K只可以取整数值，请问怎么实现非整数的情况，如K=1:0.01:3
站长推荐 /3
Simulink工具定制实现高效模型验证
MATLAB中文论坛是全球最大的 MATLAB & Simulink 中文社区。用户免费注册会员后，即可下载代码，讨论问题，请教资深用户及结识书籍作者。立即注册加入我们吧！
MATLAB官方社交平台
MATLAB中文论坛微社区菜鸟求助，一些比较复杂的微分方程用matlab的dsolve也解不了，有没有其他方法呢【微信】
全部答案（共1个回答）
冬糖极蜜、山野杂花,梅花，然后做成的糖，希望对你有用，麻烦给与好评，谢谢
这个是成立合作社的程序，你说的具体多少亩我不太清楚，（b2cf.cn）但是你可以去这里面看看，我不知道有没有，我一直都关注这个的，也很有用，希望能帮到你，还望采...
一旦开始办护照了,立马同时准备签证需要的材料,可能还来得及呦!
如果是函数的形式，可以用diff求出对应的导数，按照dy/dx=dy/dt/(dx/dt) ...的形式求出来
答: 我也移植后第37天，一个单卵单胎，一个单卵双胎，必需需要减胎吗
答: 对于那些有志于穷尽数学奥秘的学生，他总是循循善诱地予以启发和教育，而对于那些急功近利、在学习上不肯刻苦钻研的人，则毫不客气地予以批评
答: 友情帮顶,祝楼主早日找到自己想要的答案.
祝你身体健康,笑口常开!!!
答: 中国人的数学理应比外国人好！ 这是我的个人观点，这在于中国人对数字的发音是单音，因此，对数字的记忆较为简单，提高了学习数学的效率！
而科学的发展，往往受制于社会...
海鸟的种类约350种，其中大洋性海鸟约150种。比较著名的海鸟有信天翁、海燕、海鸥、鹈鹕、鸬鹚、鲣鸟、军舰鸟等。海鸟终日生活在海洋上，饥餐鱼虾，渴饮海水。海鸟食量大，一只海鸥一天要吃6000只磷虾，一只鹈鹕一天能吃(2～2.5)kg鱼。在秘鲁海域，上千万只海鸟每年要消耗?鱼400×104t，它们对渔业有一定的危害，但鸟粪是极好的天然肥料。中国南海著名的金丝燕，用唾液等作成的巢被称为燕窝，是上等的营养补品。
tann转成假名就是たん，拼音的话，相当于tang吧……
bakka转成假名是ばっか
kkou转成假名是っこう
benn转成假名是べん
kyo转成假名是きょ
系统学过五十音的话，看罗马音就能直接读了。用拼音来学是不合适的，有些发音没法用拼音标的。
P.S.罗马音里，“nn”就是波音“ん”，“kka”是前面带促音小写“っ”的“か”，同理“tta”就等于“った”。
做鲫鱼汤很重要的一点是注意火候的把握。
步骤如下：
买新鲜现杀的鲫鱼两条，个头要适中。洗的时候要把鱼鳞全部弄干净，鱼肚里也要洗净，免得汤有腥味；
洗好后，在鱼身上涂抹适当食盐，腌放十分钟；
准备好香葱三根，洗净，打结备用；
切好姜片若干（根据鱼的大小和量）；
均匀涂抹姜汁于锅内（防止鱼皮粘锅），倒入色拉油，点火；
油不宜太热，将火旋小，轻轻放鱼入锅，同时放入姜片，把火调大；
煎至鱼皮微露金黄色，将鱼轻轻翻身，直至也微呈金黄色；
煎的过程中，注意转动锅，使鱼均匀煎透；
把火调小，加冷水至淹没鱼为止，放入备好的葱结，开大火，煮沸；
把鱼翻身，再煮五分钟，放入适量的盐，继续煮，直至汤呈现奶白色；
加味精，煮两分钟。
同时准备好吃鱼的料：蘸鱼的陈醋少许倒入碗中，放少许盐，糖，味精，搅拌均匀。
将鱼单独盛在大碗里，鲫鱼汤盛在汤碗里；鱼蘸着料吃，汤即喝。
^_^，美味的鲫鱼汤呈现在你的眼前了，还有香喷喷的鱼肉……
有2个拉电子基团啊，怎么稳定，联想一下“五马分尸”~~
考虑是由于天气比较干燥和身体上火导致的，建议不要吃香辣和煎炸的食物，多喝水，多吃点水果，不能吃牛肉和海鱼。可以服用（穿心莲片，维生素b2和b6）。也可以服用一些中药，如清热解毒的。
确实没有偿还能力的，应当与贷款机构进行协商，宽展还款期间或者分期归还； 如果贷款机构起诉到法院胜诉之后，在履行期未履行法院判决，会申请法院强制执行； 法院在受理强制执行时，会依法查询贷款人名下的房产、车辆、证券和存款；贷款人名下没有可供执行的财产而又拒绝履行法院的生效判决，则有逾期还款等负面信息记录在个人的信用报告中并被限制高消费及出入境，甚至有可能会被司法拘留。
第一步：教育引导
不同年龄阶段的孩子“吮指癖”的原因不尽相同，但于力认为，如果没有什么异常的症状，应该以教育引导为首要方式，并注意经常帮孩子洗手，以防细菌入侵引起胃肠道感染。
第二步：转移注意力
比起严厉指责、打骂，转移注意力是一种明智的做法。比如，多让孩子进行动手游戏，让他双手都不得闲，或者用其他的玩具吸引他，还可以多带孩子出去游玩，让他在五彩缤纷的世界里获得知识，增长见识，逐渐忘记原来的坏习惯。对于小婴儿，还可以做个小布手套，或者用纱布缠住手指，直接防止他吃手。但是，不主张给孩子手指上“涂味”，比如黄连水、辣椒水等，以免影响孩子的胃口，黄连有清热解毒的功效，吃多了还可导致腹泻、呕吐。
合肥政务区网络广告推广网络推广哪家公司比较好 一套能在互联网上跑业务的系统，被网络营销专家赞为目前最 有效的网络推广方式！
1、搜索引擎营销：分两种SEO和PPC，即搜索引擎优化，是通过对网站结构、高质量的网站主题内容、丰富而有价值的相关性外部链接进行优化而使网站为用户及搜索引擎更加友好，以获得在搜索引擎上的优势排名为网站引入流量。
良工拥有十多位资深制冷维修工程师，十二年生产与制造经验，技术力量雄厚，配有先进的测试仪器，建有系列低温测试设备，备有充足的零部件，包括大量品牌的压缩机，冷凝器，蒸发器，水泵，膨胀阀等备品库，能为客户提供迅捷，优质的工业冷水机及模温机维修和保养。
楼主，龙德教育就挺好的，你可以去试试，我们家孩子一直在龙德教育补习的，我觉得还不错。
成人可以学爵士舞。不过对柔软度的拒绝比较大。　　不论跳什么舞，如果要跳得美，身体的柔软度必须要好，否则无法充分发挥出理应的线条美感，爵士舞也不值得注意。在展开暖身的弯曲动作必须注意，不适合在身体肌肉未几乎和暖前用弹振形式来做弯曲，否则更容易弄巧反拙，骨折肌肉。用静态方式弯曲较安全，不过也较必须耐性。柔软度的锻炼动作之幅度更不该超过疼痛的地步，肌肉有向上的感觉即可，动作(角度)保持的时间可由10馀秒至30-40秒平均，时间愈长对肌肉及关节附近的联结的组织之负荷也愈高。
正在加载...
Copyright &
Corporation, All Rights Reserved
确定举报此问题
举报原因(必选)：
广告或垃圾信息
激进时政或意识形态话题
不雅词句或人身攻击
侵犯他人隐私
其它违法和不良信息
报告，这不是个问题
报告原因(必选)：
这不是个问题
这个问题分类似乎错了
这个不是我熟悉的地区
相关问答：1234566666666菜鸟求助 编写matlab s函数，运行时出现Undefined function or variable &#39;p&#39;.这样的错误，为何？_百度知道
菜鸟求助 编写matlab s函数，运行时出现Undefined function or variable &#39;p&#39;.这样的错误，为何？
我的代码是：function[sys,x0,str,ts]=fm0code(t,x,u,flag,p)switchflag,%判断仿真状态case0,%初始化[sys,xO,str,ts]=mdlInitializeSizes(p);case1,%导数计算sys=mdlDerivatives(t,x,...
我的代码是：function[sys,x0,str,ts]=fm0code(t,x,u,flag,p)switch flag,
%判断仿真状态 case 0,
%初始化[sys,xO,str,ts]=mdlInitializeSizes(p);case 1,
%导数计算sys=mdlDerivatives(t,x,u);case 2,
%更新离散状sys=mdlUpdate(t,x,u);case 3,
%计算输出sys=mdlOutputs(t,x,u);case 4,
%计算下一个采样时刻sys=mdlGetTimeOfNextVarHit(t,x,u);case 9,sys=mdlTerminate(t,x,u);
%仿真结束otherwise
error([&#39;Unhandled flag=&#39;,num2str(flag)]);%错误处理endfunction[sys,x0,str,ts]=mdlInitializeSizes(p) %定义基本属性sizes =
%返回空的SIZE结构sizes.NumContStates=O;
%连续状态数目为0sizes.NumDiscStates=2;
%离散状态数目为2sizes.NumOutputs=1;
%输出信号数目为1sizes.Numlnputs=1;
%输入信号数目为1sizes.DirFeedthrough=1;
%直接反馈sizes.NumSampleTimes=1;
%采样周期数目为1sys = simsizes(sizes);
%初始化离散状态str = [];ts =[p 0];
%记录采样时间和模块的延时时间function sys=mdlDerivatives(t,x,u) %计算离散状态变量的导数sys = []function sys=mdlUpdate(t,x,u) %更新离散状变量、采样时间和最大步长x(1)=x(1)+1; if rem(x(1),2)==0
%每lbit周期开始电平反向
if x(2)==l
%输出为0，电平跳变
if x(2)==O
end endendsys=x;function sys=mdlOutputs(t,x,u)
%计算输出sys=x(2);function sys=mdlGetTimeOfNextVarHit(t,x,u)
%计算下一次仿真时间sampleTime=1;sys=t+sampletimefunction sys=mdlTerminate(t,x,u)
%仿真结束sys=[];这个p是用于simulink的
s函数，还需要定义吗？我看了好多例子都没定义我是这样用的仿真的时候出现：Error evaluating parameter &#39;p&#39; in &#39;FM0code1/FM0code&#39;: Undefined function or variable &#39;p&#39;.初学中，希望大师帮忙仿真模型如图
&#xe6b9;答题抽奖
首次认真答题后
即可获得3次抽奖机会，100%中奖。
江湖人称雪姐知道合伙人
江湖人称雪姐
擅长：暂未定制
我现在也被这个FM0编码解码给难住了，不知道p参数怎么就错了。请问您现在知道怎么设置了吗？
mechart知道合伙人
采纳数：62
获赞数：501
擅长：暂未定制
FM0code模块的实例中的参数p没有定义。你该打开模块检查参数
松少听听知道合伙人
擅长：暂未定制
你好，我在别的例子中也出现这样的错误，请问解决了吗？怎么解决的，请教一下，谢谢！
为你推荐：
其他类似问题
个人、企业类
违法有害信息,请在下方选择后提交
色情、暴力
我们会通过消息、邮箱等方式尽快将举报结果通知您。简单验证码识别(matlab)昨天晚上一个朋友给我发了一些验证码的图片，希望能有一个自动识别的程序。
我看了看这些样本，发现都是很规则的印刷体数字，而且还没有角度旋转，所以我就直接使用数字的面积和周长两个特征量来进行检测，发现效果还是蛮不错的。
在实验中，主要问题是‘6’和‘9’两个数字的面积和周长都是完全一样的，所以这时候我又添加了一个重心的特征。
有些蛋疼的是数字‘4’和‘0’竟然面积是一样的。。。所以只好再引入椭圆离心率特征来判别。
图像也没有怎么预处理，就是简单的阈值化，再使用ostu分割。
具体代码如下
部分实验结果
代码里主要函数
功能：regionprops(英文get the properties of region)用来度量图像区域属性的函数.常用来统计被标记的区域的面积分布，显示区域总数。
语法：STATS = regionprops(L,properties)
描 述：测量标注矩阵 L中每一个标注区域的一系列属性。L 中不同的正整数元素对应不同的区域， 例如：L 中等于整数1的元素对应区域1；L 中等于整数2的元素对应区域2;以此类推。
返回值STATS 是一个长度为 max(L(:))的结构数组,结构数组的相应域定义了每一个区域相应属性下的度量。 properties 可以是由逗号分割的字符串列表，包含字符串的单元数组，单个字符串 'all' 或者 'basic'。如果 properties 等于字符串 'all'，则所有下述字串列表中的度量数据都将被计算，如果properties 没有指定或者等于 'basic'，则属性 'Area'、'Centroid' 和'BoundingBox' 将被计算。
'Area'图像各个区域中像素总个数
　　'BoundingBox' 包含相应区域的最小矩形
　　'Centroid' 每个区域的质心（重心）
　　'MajorAxisLength' 与区域具有相同标准二阶中心矩的椭圆的长轴长度（像素意义下）
　　'MinorAxisLength' 与区域具有相同标准二阶中心矩的椭圆的短轴长度（像素意义下）
　　'Eccentricity' 与区域具有相同标准二阶中心矩的椭圆的离心率（可作为特征）
　　'Orientation' 与区域具有相同标准二阶中心矩的椭圆的长轴与x轴的交角（度）
　　'Image' 与某区域具有相同大小的逻辑矩阵
　　'FilledImage' 与某区域具有相同大小的填充逻辑矩阵
　　'FilledArea' 填充区域图像中的on像素个数
　　'ConvexHull' 包含某区域的最小凸多边形
　　'ConvexImage' 画出上述区域最小凸多边形
　　'ConvexArea' 填充区域凸多边形图像中的on像素个数
　　'EulerNumber' 几何拓扑中的一个拓扑不变量——欧拉数
　　'Extrema' 八方向区域极值点
　　'EquivDiameter' 与区域具有相同面积的圆的直径
　　'Solidity' 同时在区域和其最小凸多边形中的像素比例
　　'Extent' 同时在区域和其最小边界矩形中的像素比例
　　'PixelIdxList' 存储区域像素的索引下标
　　'PixelList' 存储上述索引对应的像素坐标
　　'Perimeter' 图像各个区域边界地区的周长
这个问题里面的样本实在是太简单，对于有旋转的，大小不同的数字验证码识别应该对数字提取旋转不变的特征点信息，以及离心率，方向像素直方图等多种信息融合才能得到较好结果。
另外还有一点就是图像的预处理，因为更复杂的验证码包含各种噪声，这就要求字符检测器有更好的鲁棒性。
阅读(...) 评论()