Question

matlab中如何画出这些数据的图形，要连续的曲线。

x=[58 56 54 48 40 33 20 0 -20 -40 -48 -50 -53 -58 56 53 50 48 20 0 36 -48 -54 -55 -55];
>> y=[1400 1200 1000 500 150 0 -280 -450 -600 -800 -1000 -1100 -1200 -1400 1300 1200 1100 1000 650 500 0 -500 -900 -1100 -1200];

高手们帮帮忙
matlab中如何画出这些数据的图形，要连续的曲线。不需要散点图。听说需要三次样条插值

Answer 1

直接使用plot(x,y)命令，默认的连线是实线，不会出现间断部分的！插值是为了在少量的数据中产生更多能够拟合原来趋势的数据序列，使得它们看起来更光滑，意思就是这样了！
如果你想插值的话：比如步长为0.5，数据将进扩大两倍（其实少一个）
x1=x(1):.5:x(end);
y1=interp1(x,y,x1,'spline');
plot(x1,y1.'-.*'); 为了能让你看清数据点，使用了'-.*'

Answer 2

x(21)有问题，应该是-36，不是36.
clear;clc;close all;
warning off all;
x=[58 56 54 48 40 33 20 0 -20 -40 -48 -50 -53 -58 56 53 50 48 20 0 -36 -48 -54 -55 -55];
y=[1400 1200 1000 500 150 0 -280 -450 -600 -800 -1000 -1100 -1200 -1400 1300 1200 1100 1000 650 500 0 -500 -900 -1100 -1200];
%数据构成循环，并作出反函数：
X=[y(1:14),y(end:-1:15),y(1)];
Y=[x(1:14),x(end:-1:15),x(1)];
%绘制原图形的概貌
plot(X,Y);
hold on;
%拟合成两条三次曲线:
for n=1:3:12,
X1=X(n:n+3);
Y1=Y(n:n+3);
A1=polyfit(X1,Y1,3);
%提取系数：
a1=A1(1);
b1=A1(2);
c1=A1(3);
d1=A1(4);
%画拟合曲线：
X1=min(X1):max(X1);
plot(X1,a1.*X1.^3+b1.*X1.^2+c1.*X1+d1,'r');
end;
for n=15:3:length(X)-3,
X2=X(n:n+3);
Y2=Y(n:n+3);
A2=polyfit(X2,Y2,3);
a2=A2(1);
b2=A2(2);
c2=A2(3);
d2=A2(4);
X2=min(X2):max(X2);
plot(X2,a2*X2.^3+b2*X2.^2+c2.*X2+d2,'r');
end;
plot(X(13:15),Y(13:15),'r');
plot(X(24:end),Y(24:end),'r');
hold off;

Answer 3

你的x和y是一个点的坐标呢，还是分开的两串数据？

追问

x，y是点的坐标，我已经在matlab中输入好了，方便复制

追答

clc;
clear;
x=[58 56 54 48 40 33 20 0 -20 -40 -48 -50 -53 -58 56 53 50 48 20 0 36 -48 -54 -55 -55];
y=[1400 1200 1000 500 150 0 -280 -450 -600 -800 -1000 -1100 -1200 -1400 1300 1200 1100 1000 650 500 0 -500 -900 -1100 -1200];
z=x+i*y;
tx1=1:length(x);
tx2=1:0.01:length(x);
ty1=1:length(y);
ty2=1:0.01:length(y);
xx=interp1(tx1,x,tx2,'cubic');
yy=interp1(ty1,y,ty2,'cubic');
zz=xx+i*yy;
figure(1)
hold on
for idx=1:length(x)
plot(real(z(idx)),imag(z(idx)),'b*')
text(real(z(idx)),imag(z(idx)), int2str(idx));
end
plot(zz,'-b')

**************************************
看这个吧，刚才写的画图的时候有点不好。这个看的清楚点

追问

还是不对，这个是磁滞回线的实验数据，你可以在网上搜下这个图形的样子。
谢谢了。采纳你定了，就算做不出来。
再次谢谢
http://image.baidu.com/i?ct=503316480&z=&tn=baiduimagedetail&word=%B4%C5%D6%CD%BB%D8%CF%DF&in=11394&cl=2&lm=-1&pn=0&rn=1&di=29239117545&ln=1998&fr=&fmq=&ic=&s=&se=&sme=0&tab=&width=&height=&face=&is=&istype=#pn0&-1

追答

我知道了
有两个问题
1我以为你给的数据是一条线，但其实是两条，应该分开画，首尾连接的地方不需要做平滑的。
2.你的数据有个点错了，x的最后第5点，我猜应该是-36
等我2分钟，马上搞定
*********************************
clc;
clear;
x=[58 56 54 48 40 33 20 0 -20 -40 -48 -50 -53 -58 56 53 50 48 20 0 -36 -48 -54 -55 -55];
y=[1400 1200 1000 500 150 0 -280 -450 -600 -800 -1000 -1100 -1200 -1400 1300 1200 1100 1000 650 500 0 -500 -900 -1100 -1200];
z=x+i*y;
x1=x(1:14);
x2=x(15:end);
y1=y(1:14);
y2=y(15:end);

tx1=1:length(x1);
txx1=1:0.01:length(x1);
tx2=1:length(x2);
txx2=1:0.01:length(x2);
ty1=1:length(y1);
tyy1=1:0.01:length(y1);
ty2=1:length(y2);
tyy2=1:0.01:length(y2);
xx1=interp1(tx1,x1,txx1,'cubic');
xx2=interp1(tx2,x2,txx2,'cubic');
yy1=interp1(ty1,y1,tyy1,'cubic');
yy2=interp1(ty2,y2,tyy2,'cubic');
zz1=xx1+i*yy1;
zz2=xx2+i*yy2;
figure(1)
hold on
for idx=1:length(x)
plot(real(z(idx)),imag(z(idx)),'b*')
text(real(z(idx)),imag(z(idx)), int2str(idx));
end
plot(zz1,'-b');
plot(zz2,'-b');