| 查看: 900 | 回复: 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客户端 |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Multiwfn支持的弱相互作用的分析方法概览
已经有10人回复
- 程序哪出错了?跪求
已经有2人回复
- 分享matlab的一些代码（入门贴）
已经有19人回复
- 【软件】Mathmatica 软件 7.0版
已经有68人回复
- 【软件】一个科研数据处理软件和一个化学分子量，元素含量，配平等计算软件
已经有83人回复
- 【zhou2009个人文集】关于HF、H2O的电荷计算
已经有15人回复
- 【sobereva个人文集】电子定域性的图形分析
已经有56人回复
- 有意义的硕士研究生阶段如何充实的度过
已经有260人回复
- 催化基础知识普及、探讨帖之一：催化剂表面酸碱性表征
已经有171人回复
6楼2017-06-05 09:31:57
2楼2017-06-04 02:43:46
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31014.8
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 855.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
5楼2017-06-04 22:53:55
7楼2017-06-05 09:45:54
Torreso98楼2017-06-05 10:04:07
9楼2017-06-05 10:08:39
10楼2017-06-05 10:22:43