这是我的代码:
N=100;
n=50;
tau=0.001;
h=0.01;
lambda=tau/h;
mu=lambda/2;
u=zeros(N,n);
u1=zeros(N,n);
u2=zeros(N,n);
phi=zeros(n,1);
for i=1:n
for j=1:N
u(j,i)=cos(2*pi*i*(j-1)*h);
u1(j,i)=cos(2*pi*i*((j-1)*h-tau));
end
for j=2:N
u2(j,i)=u(j,i)-lambda*(u(j,i)-u(j-1,i));
end
u2(1)=0;
phi(i,1)=2*pi*i/N;
end
uf=zeros(n,1);
uf1=zeros(n,1);
uf2=zeros(n,1);
for i=1:n
for j=1:N
uf(i,1)=uf(i,1)+(u(j,i)*exp(-1i*(j-1)*phi(i,1)))/100;
uf1(i,1)=uf1(i,1)+u1(j,i)*exp(-1i*j*phi(i,1))/100;
uf2(i,1)=uf2(i,1)+(u2(j,i)*exp(-1i*(j-1)*phi(i,1)))/100;
end
end
final=zeros(n,1);
for i=1:n
final(i,1)=-(h/(1i*tau))*(log(uf2(i)/uf(i)));
end
figure;
hold on
z=1:1:n;
b = real(final(z,1));
%plot(phi(z,1),b,'o');
c = imag(final(z,1));
%plot(phi(z,1),c,'-');
%plot(phi(z,1),0,'-');
plot(phi(z,1),b,'ro',phi(z,1),c,'ko',phi(z,1),0,'k-');
legend('Real','Imaginary');
legend ('Location','NorthWest');
xlabel('Reduced Wavenumber')
ylabel('Modified Wavenumber')
我正在y = 0处绘制一条线作为参考。我不希望他们在传说中。但我得到这个数字: 我该如何解决这个问题?
答案 0 :(得分:5)
不是用一个plot
- 命令绘制所有内容,而是这样做:
plot(phi(z,1),b,'ro'); hold on
plot(phi(z,1),c,'ro'); hold on
plot(phi(z,1),0,'k-'); hold off
legend('Real','Imaginary','Location','NorthWest');
xlabel('Reduced Wavenumber')
ylabel('Modified Wavenumber')
那时你不应该有问题。
实际原因是,在这种情况下,图例的调用方式不同。更优雅的解决方案提供answer of Magla.
答案 1 :(得分:3)
matlab中的传说应该像这样调用
plot(phi(z,1),b,'ro',phi(z,1),c,'ko',phi(z,1),0,'k-');
legend( {'Real','Imaginary'} , 'Location', 'NorthWest');
其中图例标签存储在字符串{...}
的单元格数组中。 Location
是一个参数,可以通过图例调用(如上面的代码中)放置,也可以使用set
函数放置在函数外部
h = legend({'Real','Imaginary'});
set(h, 'Location','NorthWest');
这给出了