24小时热门版块排行榜    

查看: 2286  |  回复: 10

Darlingzhou

铜虫 (初入文坛)

[求助] 美女求助:matlab中数组计算 画图

问题描述:
f=1:5:50
y=cos(2*pi*f*t)+cos(3*pi*f*t)
f是频率,t是时间,一组f值对应一组y,而y又是时间t的函数,应该怎样得到y的幅值随频率f的变化曲线?

[f,t]=meshgrid(1:5:50,0:0.01:1);
u=cos(2*pi*f*t)+cos(3*pi*f*t);
mesh(f,t,u);
这样为什么画不出来?
哪位高手帮忙指点下,非常感谢!!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

haoshihuang

铁杆木虫 (著名写手)

别梦依稀

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 谢谢应助! 2013-01-16 20:44:56
Darlingzhou: 金币+5, 有帮助, 谢谢,金币有限,多多见谅哈 2013-01-16 21:30:00
[f,t]=meshgrid(1:5:50,0:0.01:1);
u=cos(2*pi*f.*t)+cos(3*pi*f.*t);
mesh(f,t,u);

f后面打个".",就可以了.
爱是因为你,美丽被还原!
2楼2013-01-16 19:34:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

Darlingzhou

铜虫 (初入文坛)

引用回帖:
2楼: Originally posted by haoshihuang at 2013-01-16 19:34:00
=meshgrid(1:5:50,0:0.01:1);
u=cos(2*pi*f.*t)+cos(3*pi*f.*t);
mesh(f,t,u);

f后面打个".",就可以了.

如果t不是0:0.01:1这样变化呢,t是连续变化的时间,应该怎么画呢
3楼2013-01-16 21:31:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tianwaiflyke

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +1
Darlingzhou: 金币+3, 有帮助, 金币有限,请见谅哈 2013-01-18 17:21:23
xiegangmai: 金币+1, 谢谢参与! 2013-01-24 20:40:41
引用回帖:
3楼: Originally posted by Darlingzhou at 2013-01-16 21:31:38
如果t不是0:0.01:1这样变化呢,t是连续变化的时间,应该怎么画呢...

我觉的应该只能是离散的数据组。取步长更小些,应该就可以认为是连续的了。计算机处理貌似就是这样。
立志不存忧世,虽高无益;为文不合大义,虽工弗取!
4楼2013-01-17 04:48:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tider

木虫 (正式写手)

一切都不是连续的。。
5楼2013-01-17 07:05:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Darlingzhou

铜虫 (初入文坛)

引用回帖:
5楼: Originally posted by tider at 2013-01-17 07:05:34
一切都不是连续的。。

没这么绝对吧,我们要用辩证的眼光看实物,嘿嘿~
6楼2013-01-18 17:10:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Darlingzhou

铜虫 (初入文坛)

引用回帖:
4楼: Originally posted by tianwaiflyke at 2013-01-17 04:48:02
我觉的应该只能是离散的数据组。取步长更小些,应该就可以认为是连续的了。计算机处理貌似就是这样。...

问题:f是频率,t是时间,一组f值对应一组y,而y又是时间t的函数,应该怎样得到y的幅值随频率f的变化曲线?

这个程序能处理以上问题吗?t怎么办啊?
7楼2013-01-18 17:23:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tianwaiflyke

金虫 (正式写手)

【答案】应助回帖


xiegangmai: 金币+1, 谢谢参与! 2013-01-24 20:40:50
你的这个问题画出来的应该是一个三维图像吧?可以往三维数组上想想办法,看能不能解决。
立志不存忧世,虽高无益;为文不合大义,虽工弗取!
8楼2013-01-19 01:05:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dreamerEE

金虫 (初入文坛)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
xiegangmai: 金币+2, 谢谢参与! 2013-01-24 20:41:00
楼主这貌似波的衍射图
效果图1
 
效果图2

俯视图

x=0.05

x=0.5

x=1


代码:
clc;
clear;
PI=3.14159265;
Tmin=0;
Tmax=1;
Nt=500;
Fmin=1;
Fmax=50;
Nf=800;
t=linspace(Tmin,Tmax,Nt);
f=linspace(Fmin,Fmax,Nf);
[F,T]=meshgrid(f,t);
u=cos(2*PI.*F.*T)+cos(3*PI.*F.*T);
figure;
mesh(F,T,u);
xlabel('f');
ylabel('t');
zlabel('u');

tx=0.05;
t=linspace(tx,tx,Nt);
fx=linspace(Fmin,Fmax,Nf);
[Fx,Tx]=meshgrid(fx,t);
u=cos(2*PI.*Fx.*Tx)+cos(3*PI.*Fx.*Tx);
ux=u(1,;
figure;
plot(fx,ux);
xlabel('fx');
ylabel('ux');
title(['x=' num2str(tx)]);
9楼2013-01-19 09:46:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dreamerEE

金虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★
Darlingzhou: 金币+5, 有帮助, 我做的不是波的衍射图,不过谢谢你的热心帮助,金币有限,请多见谅哈~~ 2013-01-19 14:51:21
xiegangmai: 金币+1, 帖子中代码部分可以使用“插入代码”功能 2013-01-24 20:41:24
第一次发图片,第一张图片没弄好,但是回复又没有预览功能,重发一次,见谅。
效果图1

效果图2

俯视图

x=0.05

x=0.5

x=1

实现代码:
clc;
clear;
PI=3.14159265;
Tmin=0;
Tmax=1;
Nt=500;
Fmin=1;
Fmax=50;
Nf=800;
t=linspace(Tmin,Tmax,Nt);
f=linspace(Fmin,Fmax,Nf);
[F,T]=meshgrid(f,t);
u=cos(2*PI.*F.*T)+cos(3*PI.*F.*T);
figure;
mesh(F,T,u);
xlabel('f');
ylabel('t');
zlabel('u');

tx=0.05;
t=linspace(tx,tx,Nt);
fx=linspace(Fmin,Fmax,Nf);
[Fx,Tx]=meshgrid(fx,t);
u=cos(2*PI.*Fx.*Tx)+cos(3*PI.*Fx.*Tx);
ux=u(1,;  %这个为:ux=u(1,,系统对特定的符号进行转换自动插入了表情。
figure;
plot(fx,ux);
xlabel('fx');
ylabel('ux');
title(['x=' num2str(tx)]);
10楼2013-01-19 09:55:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Darlingzhou 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ANUsWpSnuWQD 2026-08-06 4/200 2026-08-07 04:46 by nSM3ys5fefM3
[基金申请] 听说今天filecode变了 +19 布布和一二 2026-08-06 34/1700 2026-08-07 03:13 by 虫友是什么虫
[基金申请] 基金中了 +6 laoda193707 2026-08-06 6/300 2026-08-06 23:34 by dragonxp
[基金申请] filecode +10 等待解的谜 2026-08-06 14/700 2026-08-06 21:16 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +3 archvillain 2026-08-06 7/350 2026-08-06 21:05 by archvillain
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 6/300 2026-08-06 20:19 by 苏知砚
[有机交流] 一个有机合成实验室都需要哪些设备? 50+3 kf2781974 2026-07-31 12/600 2026-08-06 15:11 by eddyin
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[教师之家] 咨询面上基金 +4 李长云 2026-07-31 7/350 2026-08-06 11:20 by 李长云
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +9 archvillain 2026-08-04 11/550 2026-08-05 20:06 by wlwhappy
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +22 低垂的野花 2026-07-31 30/1500 2026-08-05 18:03 by 低垂的野花
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
[考博] 【2027博士申请】纳米药物递送方向 20+3 13586093586 2026-08-03 4/200 2026-08-05 09:59 by lfy8008
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 18/900 2026-08-04 14:42 by archvillain
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见