应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » matlab小白，发帖求助，编好的程序出了问题，附编程
161/2返回列表12下一页
查看: 1531  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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;

>>
matlab小白，发帖求助，编好的程序出了问题，附编程
题目.jpg


matlab小白，发帖求助，编好的程序出了问题，附编程-1
程序1.jpg


matlab小白，发帖求助，编好的程序出了问题，附编程-2
程序2.jpg

[ Last edited by 309065816 on 2013-7-10 at 20:54 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-07-10 20:48:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)
在线坐等解答。。。
回复此楼
2楼2013-07-10 20:55:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★
感谢参与，应助指数 +1
jjdg: 金币+1, 感谢参与 2013-07-11 04:11:09
309065816: 金币+2 2013-07-11 09:15:27
你的a，b，f都没有被定义。
a和b是积分区间的左右端点数值。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
3楼2013-07-10 21:52:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★
jjdg: 金币+1, 感谢解释 2013-07-11 04:11:26
309065816: 金币+3 2013-07-11 09:15:33
引用回帖:
2楼: Originally posted by 309065816 at 2013-07-10 20:55:21
在线坐等解答。。。

这两个程序，都不是直接运行的。
而是把m文件分别保存在MATLAB的工作目录下，然后按照格式调用。
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
4楼2013-07-10 21:58:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
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)
即得结果，如附图。
matlab小白，发帖求助，编好的程序出了问题，附编程-3
附图1.jpg
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
5楼2013-07-10 22:20:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)
引用回帖:
5楼: Originally posted by 月只蓝 at 2013-07-10 22:20:21
原文的程序没有错，出现问题，在于其一程序抄写不一致，你抄写“syms……”的时候，少了一个字母"s"；其二程序调用方式不对，不是直接运行这些m文件而应该是按调用格式来。
比如，原文中梯形法的程序编写 ...

太感谢你了。之前也搜过类似的帖子，发现很多都有你的回复，大神呀。想再问一下呀，我想把结果的精确到小数点后4位怎么添加程序呢？
赞一下 回复此楼
6楼2013-07-11 09:17:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)
引用回帖:
5楼: Originally posted by 月只蓝 at 2013-07-10 22:20:21
原文的程序没有错，出现问题，在于其一程序抄写不一致，你抄写“syms……”的时候，少了一个字母"s"；其二程序调用方式不对，不是直接运行这些m文件而应该是按调用格式来。
比如，原文中梯形法的程序编写 ...

[q,s]=CombineTraprl('sqrt(4-x^2)',0,2,1.0e-4)
这是第一个程序里面最后的命令，而且按照编的程序可以看出，默认精度就是小数点后4位的，为什么出来的结果是这样的：
q =

   3.138365760415841


s =

    51
赞一下 回复此楼
7楼2013-07-11 09:22:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

月只蓝

主管区长 (职业作家)
引用回帖:
6楼: Originally posted by 309065816 at 2013-07-11 09:17:02
太感谢你了。之前也搜过类似的帖子，发现很多都有你的回复，大神呀。想再问一下呀，我想把结果的精确到小数点后4位怎么添加程序呢？...

将vpa与digits联用
如：圆周率pi取9位有效数字：
>> digits 9
>> vpa(pi)

ans =

                                  3.14159265
赞一下 回复此楼
MATLAB、MS小问题、普通问题请发帖求助！时间精力有限，恕不接受无偿私信求助。
8楼2013-07-11 09:23:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)
引用回帖:
5楼: Originally posted by 月只蓝 at 2013-07-10 22:20:21
原文的程序没有错，出现问题，在于其一程序抄写不一致，你抄写“syms……”的时候，少了一个字母"s"；其二程序调用方式不对，不是直接运行这些m文件而应该是按调用格式来。
比如，原文中梯形法的程序编写 ...

第二个程序是要怎么写才能控制精度呢？
>> trapm(10000,0,2)

ans =

   3.141591477611321
赞一下 回复此楼
9楼2013-07-11 09:23:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

309065816

木虫 (正式写手)
引用回帖:
8楼: Originally posted by 月只蓝 at 2013-07-11 09:23:03
将vpa与digits联用
如：圆周率pi取9位有效数字：
>> digits 9
>> vpa(pi)

ans =

                                  3.14159265...

嗯，我要计算的积分里面没有圆周率，上限跟下限是2跟0。那要怎么添加程序呢？还是像你上面写的那样吗？
赞一下 回复此楼
10楼2013-07-11 09:29:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 309065816 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +15 爱看书的可乐 2026-08-26 17/850 2026-08-27 00:25 by xxniao123
[基金申请] 国自然面上复盘~欢迎讨论 +4 晴天加油 2026-08-26 5/250 2026-08-26 21:27 by 木风9012
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 今天务委会开完了，明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[职场人生] 学生撞见辅导员送外卖，第二天全班都沉默了 +3 鲸鱼融金_浙江_� 2026-08-22 3/150 2026-08-26 09:01 by zzuzxg
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭