| 查看: 2113 | 回复: 9 | ||||
wzp0098木虫 (著名写手)
|
[求助]
利用for 语句画图
|
|
K=-10:1:10; for m=1:K; K1=Length(K) w=10; K2=K1+2*w;x=0.2; r1=8;r2=0.001*r1;T=3; A=[-r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B=[-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y=A\B; y=imag(Y(3)); end plot(K,y) 我想利用for循环语句画出图形，可是无法运行，求大虾赐教！ |
回复此楼» 猜你喜欢
时间戳变了，能看出什么问题？
已经有3人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有15人回复
-
filecode=后面第一个是大写字母
已经有9人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有80人回复
-
感觉是下周放榜了
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有6人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
今天系统多次维护，明天很可能放榜！
已经有9人回复
-
售SCI一区T0P文章，我:8O.55.1.O.5.4,科目齐全,可+急
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 这个用matlab怎么画
已经有10人回复
- 急：求用Fortran或者其他语言编写一个程序，处理下面的数据！！！！！！！！！！
已经有37人回复
- 如何用matlab将这个方程代表的图画出来，最好能把画图的代码贴出来，谢谢啊！
已经有11人回复
- 做研究需要学编程吗。学什么语言好
已经有20人回复
- 前处理软件如何生成网格
已经有8人回复
- 如何利用Matlab画图呢？
已经有7人回复
- C语言编程，以便画图。
已经有5人回复
- 拟合方程并作图
已经有6人回复
- 向各位老师求助三维画图坐标轴修改及隐去部分图的问题
已经有7人回复
- matlab画图问题: 怎样用小图标或自定义图形代替+ x . o等做图
已经有10人回复
- fortran程序运行结果输出及绘图
已经有9人回复
- 如何用MATLAB求解脉冲微分方程并画图
已经有5人回复
- 求助Matlab批处理for语句循环
已经有5人回复
- 请教mathematic问题？关于解微分方程组的和画图的一些问题。
已经有19人回复
- matlab中遇到的画图小问题，程序没问题。标题要长@@@@@@@@@@@@@@@@@@@@@@
已经有8人回复
- matlab中如何通过for语句得到已有矩阵的部分数据
已经有3人回复
- 想用C++画一个正方形，该怎么做
已经有20人回复
- matlab画图求助
已经有11人回复
- MCNP5 画图的问题
已经有7人回复
- VC的画图框架中，谁用的比较多？
已经有17人回复
- 【原创】在WORD中使用Mathematica画的图形的体会
已经有13人回复
- 【求助】Mathematica画图问题
已经有7人回复
- 【求助】再次关于Mathematica画图颜色问题（急用）
已经有3人回复
- 【分享】将Mathematica绘图数据转换至Origin软件
已经有134人回复
- 【求助】gnuplot画图软件【已完结】
已经有15人回复
fjtony163
版主 (文坛精英)
米米
- 应助: 361 (硕士)
- 贵宾: 6.876
- 金币: 70612.6
- 散金: 36050
- 红花: 222
- 沙发: 616
- 帖子: 21912
- 在线: 3793.8小时
- 虫号: 1397876
- 注册: 2011-09-11
- 专业: 遗传学研究新技术与方法
- 管辖: 论文翻译
2楼2012-07-08 09:19:50
wzp0098
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 2524.4
- 散金: 39
- 红花: 13
- 帖子: 1218
- 在线: 226.6小时
- 虫号: 853822
- 注册: 2009-09-22
- 性别: GG
- 专业: 光学
|
谢谢，我就是想画出y随K1变化的图像！K1取一个值，就可以得到一个y值，利用for语句循环。做了修改，还是不行！ K1=-10:1:10; for m=1:Length(K); K1=m w=10; K2=K1+2*w;x=0.2; r1=8;r2=0.001*r1;T=3; A=[-r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B=[-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y=A\B; y=imag(Y(3)); end plot(K1,y) |
3楼2012-07-08 11:17:39
fjtony163
版主 (文坛精英)
米米
- 应助: 361 (硕士)
- 贵宾: 6.876
- 金币: 70612.6
- 散金: 36050
- 红花: 222
- 沙发: 616
- 帖子: 21912
- 在线: 3793.8小时
- 虫号: 1397876
- 注册: 2011-09-11
- 专业: 遗传学研究新技术与方法
- 管辖: 论文翻译
4楼2012-07-08 11:30:10
wzp0098
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 2524.4
- 散金: 39
- 红花: 13
- 帖子: 1218
- 在线: 226.6小时
- 虫号: 853822
- 注册: 2009-09-22
- 性别: GG
- 专业: 光学
|
根据你的意见，我做了些修改，具体如下， K=-10:0.1:10; for n=1:length(K); K1=K(n); w=10; K2=K1+2*w;x=0.5; r1=8;r2=0.001*r1;T=3; A=[-r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B=[-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y=A\B; P=imag(Y(3)); end plot(K,P) 但是得不到我想要的图形，望赐教！ 从中间的赋值到end之前的部分是矩阵求解过程。 |
5楼2012-07-09 10:29:59
fjtony163
版主 (文坛精英)
米米
- 应助: 361 (硕士)
- 贵宾: 6.876
- 金币: 70612.6
- 散金: 36050
- 红花: 222
- 沙发: 616
- 帖子: 21912
- 在线: 3793.8小时
- 虫号: 1397876
- 注册: 2011-09-11
- 专业: 遗传学研究新技术与方法
- 管辖: 论文翻译
6楼2012-07-09 10:58:10
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
wzp0098: 金币+30, ★★★很有帮助, 牛人！！ 2012-07-09 11:16:29
臭水沟: 金币+3, 谢谢交流~~ 2012-07-09 21:24:09
感谢参与，应助指数 +1
wzp0098: 金币+30, ★★★很有帮助, 牛人！！ 2012-07-09 11:16:29
臭水沟: 金币+3, 谢谢交流~~ 2012-07-09 21:24:09
|
function parafit K=-10:0.1:10; P=0*K; for n=1:length(K); K1=K(n); w=10; K2=K1+2*w;x=0.5; r1=8;r2=0.001*r1;T=3; A=[-r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B=[-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y=A\B; P(n)=imag(Y(3)); end plot(K,P) |
7楼2012-07-09 11:05:11
duotojh
金虫 (小有名气)
- 应助: 9 (幼儿园)
- 金币: 1066.6
- 红花: 5
- 帖子: 225
- 在线: 137.8小时
- 虫号: 1290636
- 注册: 2011-05-09
- 专业: 生物物理、生物化学与分子
★ ★
臭水沟: 金币+2, 谢谢交流~~ 2012-07-09 21:24:26
臭水沟: 金币+2, 谢谢交流~~ 2012-07-09 21:24:26
|
你的意图应该是每轮循环算出的P连接成曲线，但程序中保留的是最后一轮P值（单值）， 因此画出的图形是单点当然看不到（除非将单点作标记才可以看到，例如三角符标记等等） 只需添加P向量化，并修改循环语句中的P赋值这一句即可，修改程序如下： K = -10:0.1:10; P = []; % <-------------------添加此句分配P为向量 for n=1:length(K); K1 = K(n); w = 10; K2 = K1+2*w; x = 0.5; r1 = 8; r2 = 0.001*r1; T = 3; A=[ -r1 r2-r1 -i*x i*x 0 0 0 0; 0 -r2 0 0 0 0 i*T -i*T; -2*i*x -i*x -r1-i*K1*0.5 0 -i*T 0 0 0; 2*i*x i*x 0 -r1+i*K1*0.5 0 i*T 0 0; 0 0 -i*T 0 -r2-i*(K1+K2)*0.5 0 i*x 0; 0 0 0 i*T 0 -r2+i*(K1+K2)*0.5 0 -i*x; i*T 2*i*T 0 0 i*x 0 -r2-r1+i*(K1-K2)*0.5 0; -i*T -2*i*T 0 0 0 -i*x 0 -r2-r1-i*(K1-K2)*0.5]; B = [-r1;0;-i*x;i*x;0;0;i*T;-i*T]; Y = A\B; P(n) = imag(Y(3)); % <---------------------------P向量中的元素赋值 end plot(K,P) 用freemat画出的图形如下  |
8楼2012-07-09 13:06:49
wzp0098
木虫 (著名写手)
- 应助: 2 (幼儿园)
- 金币: 2524.4
- 散金: 39
- 红花: 13
- 帖子: 1218
- 在线: 226.6小时
- 虫号: 853822
- 注册: 2009-09-22
- 性别: GG
- 专业: 光学
9楼2012-07-09 18:24:16
duotojh
金虫 (小有名气)
- 应助: 9 (幼儿园)
- 金币: 1066.6
- 红花: 5
- 帖子: 225
- 在线: 137.8小时
- 虫号: 1290636
- 注册: 2011-05-09
- 专业: 生物物理、生物化学与分子
10楼2012-07-10 10:51:42
