clc,clear,close all d1=load('juli.txt'); x1=0;y1=0; x2=18.9874900741771;y2=0; x=zeros(1,92); y=zeros(1,92); for i=1:92 sqrt((x(i)-x1)^2+(y(i)-y1)^2)==d1(i:1); sqrt((x(i)-x2)^2+(y(i)-y2)^2)==d1(i:2); end x,y 为什么输出全是0啊

在你的代码中，输出全为0的问题出现在循环中对x和y的赋值部分。在MATLAB中，==是用于判断两个值是否相等的运算符，并不是用于赋值的。因此，你需要使用赋值运算符=来给x和y赋值。

另外，在循环中，你需要将计算的结果赋值给x(i)和y(i)，而不是将计算结果与d1(i:1)和d1(i:2)进行比较。修改后的代码如下：

clc, clear, close all
d1 = load('juli.txt');
x1 = 0;
y1 = 0;
x2 = 18.9874900741771;
y2 = 0;
x = zeros(1, 92);
y = zeros(1, 92);

for i = 1:92
    x(i) = sqrt(d1(i:1)^2 - (y(i)-y1)^2) + x1;
    y(i) = sqrt(d1(i:2)^2 - (x(i)-x2)^2) + y2;
end

disp([x; y]);

这样修改后，代码会计算出对应的x和y坐标，并将结果输出。希望这次能够解决你的问题！

相关问题

clc clear close all b = [0.0181,0.0543,0.0543,0.0181]; a = [1.0000,-1.7600,1.1829,-0.2781]; m = 0:length(b)-1; l = 0:length(a)-1; k = 0:500; w = (pi/500)*k; nume = b * exp(-j * m' * w); den = a * exp(-j * l' * w); H = nume ./ den; magH = abs(H); angH = angle(H); subplot(2,1,1); plot(w/pi,magH); title('Magnitude Response'); subplot(2,1,2); plot(w/pi,angH); title('Phase Response');对该代码进行注释说明

这段代码主要是实现了一个数字滤波器的频率响应的计算和绘图，具体注释如下：

clc; % 清空命令行窗口
clear; % 清空工作区变量
close all; % 关闭所有图形窗口

% 滤波器系数
b = [0.0181,0.0543,0.0543,0.0181];
a = [1.0000,-1.7600,1.1829,-0.2781];

% 构造频率向量
k = 0:500; % k为频率采样点数
w = (pi/500)*k; % 计算对应的频率值

% 计算滤波器的频率响应
m = 0:length(b)-1; % 求出分子多项式的次数
l = 0:length(a)-1; % 求出分母多项式的次数
nume = b * exp(-j * m' * w); % 计算分子多项式的频率响应
den = a * exp(-j * l' * w); % 计算分母多项式的频率响应
H = nume ./ den; % 计算滤波器的频率响应

% 绘制滤波器的频率响应曲线
magH = abs(H); % 计算幅频响应
angH = angle(H); % 计算相频响应
subplot(2,1,1); plot(w/pi,magH); title('Magnitude Response'); % 绘制幅频响应图
subplot(2,1,2); plot(w/pi,angH); title('Phase Response'); % 绘制相频响应图

总的来说，这段代码主要是用于计算和绘制数字滤波器的频率响应曲线，其中包括了幅频响应和相频响应。

clc,clear,close all x1=[0,0];x2=[18.9875,0]; f = norm(x2 - x1); v1 = (x2 - x1) / f; u = d*v1; x3 = x1 + u; 函数或变量 'd' 无法识别。 出错 Untitled3 (第 5 行) u = d*v1;

很抱歉，我在之前的回答中遗漏了一个变量d的定义。d代表了x3到x1的距离，你需要根据具体情况给d赋值。在你的代码中，你可以在计算u之前先定义d的值。下面是修改后的代码示例：

clc, clear, close all
x1 = [0, 0];
x2 = [18.9875, 0];
d = 10;  % 假设x3到x1的距离为10，根据实际情况修改该值
f = norm(x2 - x1);
v1 = (x2 - x1) / f;
u = d * v1;
x3 = x1 + u;

disp(x3);  % 输出x3的坐标

在这个示例中，我假设x3到x1的距离为10。你可以根据你的实际情况修改这个值。运行代码后，将会在命令窗口中输出x3的坐标。希望这次能帮到你！