How to plot a smooth graph?
显示 更早的评论
clc,clear all
b0 = 0.015;
b1 = 0.015;
b2 = 0.035;
b3 = 0.050;
v = 40; %velocity
L = 2;
w = (2*pi*v)/L; %angular frequency
F = w*(1/(2*pi)); %in Hz
t = linspace(0,10,18);
b = b0 + b1*sin(w*t) + b2*sin(2*w*t)
+b3*sin(3*w*t);
plot(t,b)
Hi, I want to make this graph as being a smooth graph.
Thank you for your time.
1 个评论
Rik
2020-4-2
Although you could resample your function to draw a smooth curve through your data points, that would not accurately describe your function. You can greatly increase the number of points in your t vector to see that you have an enormous amount of aliasing going on here. See this Wikipedia page for more information.
采纳的回答
You need to increase the step size of your t vector:
t = linspace(0,10,100);
Code:
b0 = 0.015;
b1 = 0.015;
b2 = 0.035;
b3 = 0.050;
v = 40; %velocity
L = 2;
w = (2*pi*v)/L; %angular frequency
F = w*(1/(2*pi)); %in Hz
t = linspace(0,10,100);
b = b0 + b1*sin(w*t)+b2*sin(2*w*t)+b3*sin(3*w*t);
plot(t,b)
5 个评论
donghun lee
2020-4-2
Rik
2020-4-2
if you increase the length of t to 2000 you will see that this will still be aliased. The code below is my attempt (mostly from memory) to apply the Nyquist criterion to determine the minimum number of samples you need for an actual smooth graph.
t_max=10;
n_samples=ceil(2*3*w*t_max);
t = linspace(0,t_max,n_samples);
If you set axis([0 1 -0.125 0.125]) you will be able to see the effect of the number of samples.
donghun lee
2020-4-2
编辑：donghun lee
2020-4-2
Rik
2020-4-2
That is why I suggested to zoom in. The weirdness is not from me, but actually from your function. If you want a true plot of your function you should increase the number of samples until further increasing it no longer changes your graph. For this function that starts at about 2500 samples. With about 10000 samples I no longer notice an improvement in smoothness.
You could avoid having to do this yourself by using fplot:
b0 = 0.015;
b1 = 0.015;
b2 = 0.035;
b3 = 0.050;
v = 40; %velocity
L = 2;
w = (2*pi*v)/L; %angular frequency
F = w*(1/(2*pi)); %in Hz
t_max=10;
b_fun =@(t) b0 + b1*sin(w*t)+b2*sin(2*w*t)+b3*sin(3*w*t);
fplot(b_fun,[0 t_max])
axis([0 1 -0.125 0.125])%zoom in to see detail
donghun lee
2020-4-2
更多回答（1 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Spline Postprocessing 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
