24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1667  |  回复: 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 的主题更新
信息提示
请填处理意见