精华区文章阅读

发信人: hahn (有奇@选择遗忘的自由), 信区: Matlab
标  题: [合集] 秋天的萝卜，再次请教一个关于曲线画图的问题！
发信站: 哈工大紫丁香 (Sat Aug  5 19:16:25 2006), 站内

────────────────────────────────────────
  zhaoxl (小刀)                    于  (Fri May 12 09:55:38 2006)  说道:

你好，秋天的萝卜，我还有一个要画图的问题是关于而重积分的自己又搞不定了，请你帮
我看看。谢谢了！

────────────────────────────────────────
  zjliu (秋天的萝卜)               于  (Fri May 12 10:26:51 2006)  说道:

下面的程序待修改[你可参考之]：

clear;clc;close all;
mu=0.115;
d0=0.115^2;
b0=0.158;
Rt=linspace(-25,10,40);
k=1;
for R=Rt;
    if R<0;
        C_R(k)=1;
    else
        [r,z]=meshgrid(linspace(0,R,200),linspace(0.1,80,500));
        F=r./z.*exp(-[log(z)-mu].^2/2/d0-(r.^2+z.^2)/2/b0);
        F=F.*besseli(0,r.*z/b0)/b0/sqrt(2*pi*d0);
        dr=r(1,2)-r(1,1);
        dz=z(2,1)-z(1,1);
        C_R(k)=1-sum(F(:))*dr*dz;
    end
    k=k+1;
end
plot(Rt,C_R)

【在 zhaoxl (小刀) 的大作中提到: 】
: 你好，秋天的萝卜，我还有一个要画图的问题是关于而重积分的自己又搞不定了，请你帮
: 我看看。谢谢了！

────────────────────────────────────────
  sammi (非洲老头跳高)             于  (Fri May 12 14:55:42 2006)  说道:

萝卜尊强

【在 zjliu (秋天的萝卜) 的大作中提到: 】
: 标  题: Re: 秋天的萝卜，再次请教一个关于曲线画图的问题！
: 发信站: 哈工大紫丁香 (Fri May 12 10:26:51 2006), 转信
:
: 下面的程序待修改[你可参考之]：
:
: clear;clc;close all;
: mu=0.115;
: d0=0.115^2;
: b0=0.158;
: Rt=linspace(-25,10,40);
: k=1;
: for R=Rt;
:     if R<0;
:         C_R(k)=1;
:     else
:         [r,z]=meshgrid(linspace(0,R,200),linspace(0.1,80,500));
:         F=r./z.*exp(-[log(z)-mu].^2/2/d0-(r.^2+z.^2)/2/b0);
:         F=F.*besseli(0,r.*z/b0)/b0/sqrt(2*pi*d0);
:         dr=r(1,2)-r(1,1);
:         dz=z(2,1)-z(1,1);
:         C_R(k)=1-sum(F(:))*dr*dz;
:     end
:     k=k+1;
: end
: plot(Rt,C_R)
:
: 【在 zhaoxl (小刀) 的大作中提到: 】
: : 你好，秋天的萝卜，我还有一个要画图的问题是关于而重积分的自己又搞不定了，请你帮
: : 我看看。谢谢了！
:
:
: --
: 学习是最重要的，抓最主要的事做
:
: ※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: luobo.yculblog.com]
:
:
: ※ 来源:·哈工大紫丁香 bbs.hit.edu.cn·[FROM: 202.118.229.162]

────────────────────────────────────────
  zjliu (秋天的萝卜)               于  (Fri May 12 14:57:31 2006)  说道:

其实这个应该zhaoxl自己去写的，我不想再写下去了
如果再找3重积分的最好模仿这两个吧~~
【在 sammi (非洲老头跳高) 的大作中提到: 】
: 萝卜尊强

────────────────────────────────────────
  feifeifool (学习、学习、再学习)  于  (Fri May 12 15:43:13 2006)  说道:

萝卜，我的意见是像数值积分的问题，简单的求和并不总是正确的，程序只能作为一种参考，
因为这种简单的累加，有可能引起累计误差，而解决这些累计误差的问题应该是数值分析里
面的东西，这个应该参考相应数值分析中的方法，自己编写程序。我没做过这方面的东西，只
是个人意见，仅供参考。
【在 zjliu (秋天的萝卜) 的大作中提到: 】
: 其实这个应该zhaoxl自己去写的，我不想再写下去了
: 如果再找3重积分的最好模仿这两个吧~~

────────────────────────────────────────
  zjliu (秋天的萝卜)               于  (Fri May 12 15:46:03 2006)  说道:

这样负责的表达式作符号解不合适, 这样求和或者用
quad 做数值解还是很实用的方法,只要表达式收敛做
累加和还是可以的

【在 feifeifool (学习、学习、再学习) 的大作中提到: 】
: 萝卜，我的意见是像数值积分的问题，简单的求和并不总是正确的，程序只能作为一种参考，
: 因为这种简单的累加，有可能引起累计误差，而解决这些累计误差的问题应该是数值分析里
: 面的东西，这个应该参考相应数值分析中的方法，自己编写程序。我没做过这方面的东西，只
: ...................

────────────────────────────────────────
  feifeifool (学习、学习、再学习)  于  (Fri May 12 15:53:32 2006)  说道:

quad只能做定上限的吧，一元积分，他的问题做解析解肯定不行。这个不是讨论的关键。
我想累计误差的和，或者说累计误差的级数，如果不收敛的话，这样的数值积分就是错误的。
这个在数值分析里讲到过的，数值计算结果的稳定性就是讲这类问题的。同时，希望有精
通数值积分的网友共同讨论、批评、指正。
【在 zjliu (秋天的萝卜) 的大作中提到: 】
: 这样负责的表达式作符号解不合适, 这样求和或者用
: quad 做数值解还是很实用的方法,只要表达式收敛做
: 累加和还是可以的
: ...................

────────────────────────────────────────
  zhaoxl (小刀)                    于  (Fri May 12 17:32:13 2006)  说道:

谢谢你，秋天的萝卜！从今天开始一定要学习matlab!希望以后还能经常和你请教问题！

【在 zjliu (秋天的萝卜) 的大作中提到: 】
: 这样负责的表达式作符号解不合适, 这样求和或者用
: quad 做数值解还是很实用的方法,只要表达式收敛做
: 累加和还是可以的

────────────────────────────────────────
  baizhf (（小白）)              于  (Sat May 13 09:39:52 2006)  说道:

【在 sammi (非洲老头跳高) 的大作中提到: 】
: 萝卜尊强

萝卜是个好同志~~~

Matlab 版 (精华区)