| 查看: 3017 | 回复: 7 | ||
yangyu1002新虫 (初入文坛)
|
[求助]
修改MATLAB语句，实现两个变量的循环 已有1人参与
|
|
clear all;clc; global hx0 hy0 hz0 % 外加磁场参数 hz0=0; hy0=0; i=0; mz3=0; mz4=0; for hx0=0:0.01:5; y0=[0.11;0.01;pi-0.11;pi+0.01]; t0=0:0.1:500; [t,y]=ode45('rhssys',t0,y0); mx1=sin(y(:,1)).*cos(y(:,2)); my1=sin(y(:,1)).*sin(y(:,2)); mz1=cos(y(:,1)); mx2=sin(y(:,3)).*cos(y(:,4)); my2=sin(y(:,3)).*sin(y(:,4)); mz2=cos(y(:,3)); i=i+1; mz3(i)=cos(y(5000,1)); mz4(i)=cos(y(5000,3)); end hx0=0:0.01:5; subplot(2,1,1); plot(hx0,mz3,'LineWidth',2); xlabel('h_x0','fontsize',20);ylabel('m_z1','fontsize',20); title('eta=0.33','fontsize',20); set(gca,'FontSize',20); subplot(2,1,2); plot(hx0,mz4,'LineWidth',2); xlabel('h_x0','fontsize',20);ylabel('m_z2','fontsize',20); set(gca,'FontSize',20); 此循环程序现在实现了hx0和t0的双循环吗?即每一个hx0，此程序算了t0从0到500所对应的mx,my,mz吗？如果没有，如何实现hx0和t0的双循环?最后怎样改变语句，可以读出所有t0对应的y0值，而不是只是最后一组y0. |
回复此楼» 猜你喜欢
计算机、0854电子信息（085401-058412）调剂
已经有4人回复
-
基金申报
已经有3人回复
-
国自然申请面上模板最新2026版出了吗？
已经有9人回复
-
溴的反应液脱色
已经有6人回复
-
纳米粒子粒径的测量
已经有7人回复
-
常年博士招收(双一流，工科)
已经有4人回复
-
推荐一本书
已经有10人回复
-
参与限项
已经有5人回复
-
有没有人能给点建议
已经有5人回复
-
假如你的研究生提出不合理要求
已经有12人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- comsol 5.0 with MATLAB运行的问题
已经有2人回复
- 第一次接触PBS新人求助..
已经有9人回复
- MATLAB程序设计与工程应用
已经有144人回复
- 关于SIFT算法特征点的提取
已经有3人回复
- Matlab心得及学习方法（不断更新）
已经有26人回复
- matlab程序设计题目求解
已经有0人回复
- matlab中if循环不起来，请求大牛帮助！
已经有11人回复
- SCI 我爱你
已经有7人回复
tuburfransth
金虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 2380.1
- 散金: 68
- 红花: 5
- 帖子: 384
- 在线: 124.8小时
- 虫号: 3425987
- 注册: 2014-09-18
- 性别: GG
- 专业: 模式识别
2楼2016-11-14 16:34:43
yangyu1002
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 750.3
- 散金: 83
- 红花: 1
- 帖子: 47
- 在线: 65.8小时
- 虫号: 4223292
- 注册: 2015-11-16
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
3楼2016-11-14 16:37:26
tuburfransth
金虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 2380.1
- 散金: 68
- 红花: 5
- 帖子: 384
- 在线: 124.8小时
- 虫号: 3425987
- 注册: 2014-09-18
- 性别: GG
- 专业: 模式识别
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
yangyu1002: 金币+40, ★★★★★最佳答案 2016-11-14 20:35:00
感谢参与，应助指数 +1
yangyu1002: 金币+40, ★★★★★最佳答案 2016-11-14 20:35:00
|
在循环外面设置一个空数组，在循环里面每计算一次就储存一次。 发自小木虫Android客户端 |
4楼2016-11-14 16:39:04
tuburfransth
金虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 2380.1
- 散金: 68
- 红花: 5
- 帖子: 384
- 在线: 124.8小时
- 虫号: 3425987
- 注册: 2014-09-18
- 性别: GG
- 专业: 模式识别
5楼2016-11-14 16:41:18
yangyu1002
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 750.3
- 散金: 83
- 红花: 1
- 帖子: 47
- 在线: 65.8小时
- 虫号: 4223292
- 注册: 2015-11-16
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
|
clear all;clc; global hx0 hy0 hz0 % 外加磁场参数 hz0=0; hy0=0; i=0; mz3=0; mz4=0; y0=[0.11;0.01;pi-0.11;pi+0.01]; for hx0=0:0.01:5; t0=0:0.1:500; [t,y]=ode45('rhssys',t0,y0); mx1=sin(y(:,1)).*cos(y(:,2)); my1=sin(y(:,1)).*sin(y(:,2)); mz1=cos(y(:,1)); mx2=sin(y(:,3)).*cos(y(:,4)); my2=sin(y(:,3)).*sin(y(:,4)); mz2=cos(y(:,3)); i=i+1; mz3(i)=cos(y(5000,1)); mz4(i)=cos(y(5000,3)); y0=[y(5000,1);y(5000,2);y(5000,3);y(5000,4)]; end |
6楼2016-11-14 16:49:23
yangyu1002
新虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 750.3
- 散金: 83
- 红花: 1
- 帖子: 47
- 在线: 65.8小时
- 虫号: 4223292
- 注册: 2015-11-16
- 性别: MM
- 专业: 凝聚态物性I:结构、力学和
|
clear all;clc; global hx0 hy0 hz0 % 外加磁场参数 hz0=0; hy0=0; i=0; mz3=0; mz4=0; y0=[0.11;0.01;pi-0.11;pi+0.01]; for hx0=0:0.01:5; t0=0:0.1:500; [t,y]=ode45('rhssys',t0,y0); mx1=sin(y(:,1)).*cos(y(:,2)); my1=sin(y(:,1)).*sin(y(:,2)); mz1=cos(y(:,1)); mx2=sin(y(:,3)).*cos(y(:,4)); my2=sin(y(:,3)).*sin(y(:,4)); mz2=cos(y(:,3)); i=i+1; mz3(i)=cos(y(5000,1)); mz4(i)=cos(y(5000,3)); y0=[y(5000,1);y(5000,2);y(5000,3);y(5000,4)]; end 这个程序每次y0为t0=5000时的值，请问如何保存每次循环计算出的y0. |
7楼2016-11-14 16:51:44
tuburfransth
金虫 (正式写手)
- 应助: 9 (幼儿园)
- 金币: 2380.1
- 散金: 68
- 红花: 5
- 帖子: 384
- 在线: 124.8小时
- 虫号: 3425987
- 注册: 2014-09-18
- 性别: GG
- 专业: 模式识别
8楼2016-11-14 17:02:05