Question

matlab的问题，大神们帮我看看这个椭圆为什么画不出来啊？

data = inputdlg({'输入跳频周期t(s)：','输入干扰机的反应时间t1(s)：','输入发射机的横坐标x1','输入发射机的纵坐标y1','输入接收机的横坐标x2','输入接收机的纵坐标y2','输入干扰机的横坐标x3','输入干扰机的纵坐标y3'},'inputdlg demo',1)
t = str2num(data{1});
t1 = str2num(data{2});
x1 = str2num(data{3});
y1 = str2num(data{4});
x2 = str2num(data{5});
y2 = str2num(data{6});
x3 = str2num(data{7});
y3 = str2num(data{8});
l=sqrt((x1-x2)^2+(y1-y2)^2);
v=300;
a=0.5*(l+v*(t-t1));
b=0.5*sqrt(v*(t-t1)*(2*l+v*(t-t1)));
c1=(x1-x2)/2;
c2=(y1-y2)/2;

f=@(x,y)(x-c1)^2/a^2+(y-c2)^2/b^2-1;

ezplot(f,200*[-x1-x1/10,x1+x1/10],200*[-x2-x2/10,x2+x2/10]);
hold on;
plot(x1,y1,'r.','Markersize',15);
plot(x2,y2,'r.','Markersize',15);
plot(x3,y3,'r.','Markersize',15);
grid on;
box on;
title('跳频电台的干扰椭圆');
主要是f=@(x,y)(x-c1)^2/a^2+(y-c2)^2/b^2-1;这一行，画椭圆的方程，哪里错了啊？要不要f=@(x,y)(x-c1).^2/a^2+(y-c2).^2/b^2-1;加一个点啊？c1 c2是椭圆的中心

Answer 1

是的，函数句柄的定义方式其实和函数是一样的

要加点的：
f=@(x,y)(x-c1).^2/a^2+(y-c2).^2/b^2-1是对的

追问

我刚才发现这个程序有一个问题
我要做的图的以(x1,y1) (x2,y2)为两个焦点。
a=0.5*(l+v*(t-t1));
b=0.5*sqrt(v*(t-t1)*(2*l+v*(t-t1)));
a和b为长短半轴。
而椭圆的中心不能用下面这个公式来求。
c1=(x1-x2)/2;
c2=(y1-y2)/2;
f=@(x,y)(x-c1)^2/a^2+(y-c2)^2/b^2-1;
怎么画这个椭圆啊

追答

换种画法吧，其实你定义的函数画的是3维的
clear all;clc;
a=5;
b=4;
xx=[-3,3];
yy=[1,2];
c1=(xx(1)+xx(2))/2;
c2=(yy(1)+yy(2))/2;
t=0:pi/40:2*pi;
x=c1+a*cos(t);
y=c2+b*sin(t);
plot(x,y);

追问

那如果我想通过自己输入坐标的值来画这个图呢，并显示出焦点，该怎么改啊

追答

你输入坐标，最后也是要算出a和b来的
a和b算出后，c就知道了：
clear all;clc;
a=5;
b=4;
cc1=-sqrt(a^2-b^2);
cc2=sqrt(a^2-b^2);
xx=[-3,3];
yy=[1,6];
c1=(xx(1)+xx(2))/2;
c2=(yy(1)+yy(2))/2;
t=0:pi/40:2*pi;
x=c1+a*cos(t);
y=c2+b*sin(t);
plot(x,y);hold on;
plot(cc1,c2,'r.','Markersize',20);
plot(cc2,c2,'r.','Markersize',20);