24小时热门版块排行榜

>论坛更新日志 (3778)
>文献求助 (379)
>导师招生 (302)
>虫友互识 (289)
>硕博家园 (153)
>博后之家 (147)
>休闲灌水 (137)
>考博 (122)
>论文投稿 (109)
>基金申请 (83)
>招聘信息布告栏 (50)
>教师之家 (48)
>考研 (45)
>公派出国 (39)
>绿色求助(高悬赏) (37)
>找工作 (33)

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

309065816

木虫 (正式写手)

应助: 2 (幼儿园)
金币: 2494.2
散金: 180
红花: 3
帖子: 516
在线: 142.3小时
虫号: 1227451
注册: 2011-03-09
性别: MM
专业: 电化学

[求助] 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 ]

回复此楼

» 猜你喜欢

到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有7人回复
申请2026年博士已经有5人回复
天津工业大学郑柳春团队欢迎化学化工、高分子化学或有机合成方向的博士生和硕士生加入已经有5人回复
寻求一种能扛住强氧化性腐蚀性的容器密封件已经有6人回复
2025冷门绝学什么时候出结果已经有7人回复
请问有评职称，把科研教学业绩算分排序的高校吗已经有6人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有7人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有4人回复
请问下大家为什么这个铃木偶联几乎不反应呢已经有5人回复
康复大学泰山学者周祺惠团队招收博士研究生已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

小白求助自适应控制中投影算子的MATLAB实现方法已经有5人回复

1楼2013-07-10 20:48:53

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

专家经验: +1059
应助: 1712 (讲师)
贵宾: 8.888
金币: 68120.7
散金: 1938
红花: 443
沙发: 4
帖子: 4373
在线: 3291.4小时
虫号: 1122189
注册: 2010-10-14
专业: 宇宙学
管辖: 计算模拟区

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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)
即得结果，如附图。