| 查看: 1109 | 回复: 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客户端 |
回复此楼» 猜你喜欢
面上的水很深
已经有18人回复
-
国社科通讯评审
已经有5人回复
-
5个%2F，非固定段22个字母
已经有9人回复
-
生命口C13面上会评时间
已经有11人回复
-
内心匮乏
已经有16人回复
-
生命口面上会评结束了吗大家有知道结果吗？现在还没有结果正常不？
已经有9人回复
-
投英文期刊如何查重
已经有6人回复
-
不要再数国自然申请书的 filecode 的分隔符个数了
已经有21人回复
-
伯胺上甲基求助
已经有7人回复
-
跨出版社商投稿
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 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
- 金币: 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
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
10