| 查看: 1532 | 回复: 15 | ||||
309065816木虫 (正式写手)
|
[求助]
matlab小白，发帖求助，编好的程序出了问题，附编程
|
|
从文献和书上找到合适的程序，但是抄到MATLAB里面却有问题。都是相同的问题，h=(b-a)/n， b没有含义。附上题目，两个原程序和我抄的程序，均是梯形法，应该没抄错吧？ 程序1： function [q,step] = CombineTraprl(f,a,b,eps) %被积函数： f; %积分区间左端点：a; %积分区间右端点：b; %eps 精度； %积分结果：q; %step 积分的子区间数； if(nargin==3); eps=1.0e-4; end n=1; h=(b-a)/2; q1=0; q2=(subs(sym(f),findsym(sym(f)),a)+subs(sym(f),findsym(sym(f)),b))/h; while abs(q2-q1)>eps; n=n+1; h=(b-a)/n; q1=q2; q2=0; for i=0:n-1; x=a+h*i; x1=x+h; q2=q2+(h/2)*(subs(sym(f),findsym(sym(f)),x)+... subs(sym(f),findsym(sym(f)),x1)); end end q=q2; step=n; ??? Input argument "b" is undefined. Error in ==> CombineTraprl at 12 h=(b-a)/2; >> 程序2： function y = trapm(n,a,b) format long h = (b - a)/n; sym x fx fx = sqrt(4 - x^2); f0 = subs(fx,'x',a); fn = subs(fx,'x',b); sum = f0; for i = 1:n-1 xi = a + i* h; fi = subs(fx,'x',xi); sum = sum + 2* fi; end sum = sum + fn; y = h* sum/2; ??? Input argument "b" is undefined. Error in ==> trapm at 3 h = (b - a)/n; >>  题目.jpg  程序1.jpg  程序2.jpg [ Last edited by 309065816 on 2013-7-10 at 20:54 ] |
回复此楼» 猜你喜欢
有多少人是今天查系统知道结果的？
已经有17人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
基金不中，共勉
已经有5人回复
-
看板上这么多中的，有点像50人群里49个人都是骗子的那种感觉……
已经有3人回复
-
国自然面上复盘~欢迎讨论
已经有5人回复
-
我不理解！
已经有23人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有10人回复
-
系统查不到
已经有10人回复
-
怎么查啊
已经有5人回复
-
范进中举一文的中心思想
已经有10人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 小白求助自适应控制中投影算子的MATLAB实现方法
已经有5人回复
309065816
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2494.2
- 散金: 180
- 红花: 3
- 帖子: 516
- 在线: 142.3小时
- 虫号: 1227451
- 注册: 2011-03-09
- 性别: MM
- 专业: 电化学
2楼2013-07-10 20:55:21
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
3楼2013-07-10 21:52:50
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
4楼2013-07-10 21:58:25
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
309065816(jjdg代发): 金币+5, 辛苦了 2013-07-11 04:11:54
309065816: 金币+10, ★★★★★最佳答案 2013-07-11 09:15:43
309065816(jjdg代发): 金币+5, 辛苦了 2013-07-11 04:11:54
309065816: 金币+10, ★★★★★最佳答案 2013-07-11 09:15:43
|
原文的程序没有错，出现问题，在于其一程序抄写不一致，你抄写“syms……”的时候，少了一个字母"s"；其二程序调用方式不对，不是直接运行这些m文件而应该是按调用格式来。 比如，原文中梯形法的程序编写得很巧妙，只需要3个参数：积分的左右端点数值和子区间数目，就可以得到计算结果。把trapm.m存到MATLAB的工作目录下，在主窗口输入: trapm(10000,0,3*pi/2) 即得结果，如附图。  附图1.jpg |
5楼2013-07-10 22:20:21
309065816
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2494.2
- 散金: 180
- 红花: 3
- 帖子: 516
- 在线: 142.3小时
- 虫号: 1227451
- 注册: 2011-03-09
- 性别: MM
- 专业: 电化学
6楼2013-07-11 09:17:02
309065816
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2494.2
- 散金: 180
- 红花: 3
- 帖子: 516
- 在线: 142.3小时
- 虫号: 1227451
- 注册: 2011-03-09
- 性别: MM
- 专业: 电化学
7楼2013-07-11 09:22:05
月只蓝
主管区长 (职业作家)
-
专家经验: +1059 - 应助: 1712 (讲师)
- 贵宾: 8.888
- 金币: 68171.7
- 散金: 1938
- 红花: 443
- 沙发: 4
- 帖子: 4373
- 在线: 3292.4小时
- 虫号: 1122189
- 注册: 2010-10-14
- 专业: 宇宙学
- 管辖: 计算模拟区
8楼2013-07-11 09:23:03
309065816
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2494.2
- 散金: 180
- 红花: 3
- 帖子: 516
- 在线: 142.3小时
- 虫号: 1227451
- 注册: 2011-03-09
- 性别: MM
- 专业: 电化学
9楼2013-07-11 09:23:43
309065816
木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 2494.2
- 散金: 180
- 红花: 3
- 帖子: 516
- 在线: 142.3小时
- 虫号: 1227451
- 注册: 2011-03-09
- 性别: MM
- 专业: 电化学
10楼2013-07-11 09:29:17
投票:
