| 查看: 1152 | 回复: 15 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
一个简单的程序求分析 已有2人参与
|
||
|
clc clear b=0.1; sum=0; N=input('请输入N的值:'); for i=1:N x=2.0; y=0.0; r=sqrt(x*x+y*y); while(r>1&&r<3) deg = rand(1); if(deg<0.25) x=x+b; elseif(deg<0.5) x=x-b; elseif(deg<0.75) y=y+b; else y=y-b; end r=sqrt(x*x+y*y); if(r<=1) sum=sum+4;break; end if(r>=3) sum=sum+6;break; end end end 我对程序基本不太懂，论文里有个小程序，前半部分能看懂，后半部分不太明白。 发自小木虫IOS客户端 |
回复此楼» 猜你喜欢
哪位老哥知道今年的国自然具体哪一天放榜？
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有7人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有5人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有50人回复
-
filecode
已经有11人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有5人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有4人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Multiwfn支持的弱相互作用的分析方法概览
已经有10人回复
- 程序哪出错了?跪求
已经有2人回复
- 分享matlab的一些代码（入门贴）
已经有19人回复
- 【软件】Mathmatica 软件 7.0版
已经有68人回复
- 【软件】一个科研数据处理软件和一个化学分子量，元素含量，配平等计算软件
已经有83人回复
- 【zhou2009个人文集】关于HF、H2O的电荷计算
已经有15人回复
- 【sobereva个人文集】电子定域性的图形分析
已经有56人回复
- 有意义的硕士研究生阶段如何充实的度过
已经有260人回复
- 催化基础知识普及、探讨帖之一：催化剂表面酸碱性表征
已经有171人回复
5楼2017-06-04 22:53:55
2楼2017-06-04 02:43:46
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
感谢参与，应助指数 +1
|
clc clear b=0.1; sum=0; N=input('请输入N的值:'); for i=1:N x=2.0; y=0.0; %% 赋初始值 r=sqrt(x*x+y*y); while(r>1&&r<3) deg = rand(1); %%生成 0-1之间的随机数 if(deg<0.25) x=x+b; elseif(deg<0.5) x=x-b; elseif(deg<0.75) y=y+b; else y=y-b; end %%% 一系列判断 赋值 r=sqrt(x*x+y*y); if(r<=1) sum=sum+4;break; %% r<=1 输出 4 end if(r>=3) sum=sum+6;break; %% r>=3 输出 6 end end %% 1<r<3 开始下一轮循环 end 这个给N赋个值 自己运行一下就知道了 |
3楼2017-06-04 20:38:57
|
谢谢，这个程序我用来答辩，我的理解是是用来产生随机数，然后求和的,我不知道如何在答辩时老师问程序我该怎么说给他听，是一句一句像你这样翻译嘛，我一辨没过就是没说好这东西，老师就说我是抄的，没有自己动脑想，换个数据。 发自小木虫IOS客户端 |
4楼2017-06-04 22:51:49
