How to plot thisfunction?

Question

Yagnaseni Roy 2012-5-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/37555-how-to-plot-thisfunction

The function to be plotted is:

50*(y*(28.065*cos(60[1/s]*t+0.2094395101999999[1/m]*x)+22.7058*cos(2*(60[1/s]*t+0.2094395101999999[1/m]*x))+15.1368*cos(3*(60[1/s]*t+0.2094395101999999[1/m]*x))+7.014*cos(4*(60[1/s]*t+0.2094395101999999[1/m]*x)))/1[m])[m/s]

As you can see, its a function of three variables, x,y and t.....

However I need it to vary only between 0 to 30 w.r.t x...and the y and t values can be anything...what is the format of the code to be used?

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Walter Roberson 2012-5-6

First you rewrite the expression in terms of MATLAB operations. In particular, in MATLAB, [] is used for building arrays, and is not an arithmetic operator.

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

the cyclist 2012-5-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/37555-how-to-plot-thisfunction#answer_46883

在 MATLAB Online 中打开

First, as Walter mentions, you need to rewrite the definition of your function in MATLAB syntax.

Here, I assumed that you wanted x,v, and t on a simple grid, and plot them for all combinations.

[Also, it looks like the terms you have in square brackets are just units, so I simply took them out.]

xvec = 1:10;
yvec = 1:20;
tvec = 1:5;
[x,y,t] = ndgrid(xvec,yvec,tvec);
f = 50*(y.*(28.065*cos(60*t+0.2094395101999999*x)+22.7058*cos(2*(60*t+0.2094395101999999*x))+15.1368*cos(3*(60*t+0.2094395101999999*x))+7.014*cos(4*(60*t+0.2094395101999999*x)))/1);

But now you have a function three variables, so you want a plot over a four-dimensional space: (x,y,t,f).

What kind of plot are you hoping to see from that?

EDIT: Here's a version of the plotting that selects just one value of y and t, as you suggest:

xvec = 0.1:0.1:5;
yvec = 10;
tvec = 1;
[x,y,t] = ndgrid(xvec,yvec,tvec);
f = 50*(y.*(28.065*cos(60*t+0.2094395101999999*x)+22.7058*cos(2*(60*t+0.2094395101999999*x))+15.1368*cos(3*(60*t+0.2094395101999999*x))+7.014*cos(4*(60*t+0.2094395101999999*x)))/1);
figure
plot(x,f)