Tengo la siguiente transformación de frecuencia utilizada para diseñar un prototipo de paso bajo a partir de una máscara de especificación de paso de banda:
Dado que , me pregunto si al resolver para la ecuación se puede reescribir así:
En caso de duda, trama! codigo matlab:
w_array = -10:1:10;
w0_array = -10:1:10;
wl_array_1 = zeros(length(w_array), length(w0_array));
wl_array_2 = zeros(length(w_array), length(w0_array));
syms wl
for i=1:length(w_array)
w = w_array(i);
for j=1:length(w0_array)
w0 = w0_array(j);
s1 = solve(wl^2 - w0*w*wl/2.45 - w0^2 == 0, wl, 'PrincipalValue', true);
if isempty(s1)
s1 = NaN;
end
wl_array_1(i,j) = double(s1);
s2 = solve(w - 2.45*abs(wl/w0 - w0/wl) == 0, wl, 'PrincipalValue', true);
if isempty(s2)
s2 = NaN;
end
wl_array_2(i,j) = double(s2);
end
end
figure(1)
for j=1:length(w0_array)
plot(w_array,wl_array_1(:,j)), hold on
end
xlabel('\omega'), ylabel('\omega´'), xlim([-10,10]), ylim([-45,5])
title('Not considering the absolute in equation'), hold off,
figure(2)
for j=1:length(w0_array)
plot(w_array,wl_array_2(:,j)), hold on
end
xlabel('\omega'), ylabel('\omega´'), xlim([-10,10]), ylim([-45,5])
title('Considering the absolute in equation'), hold off
Como puede verse, el absoluto asegura que las soluciones solo existen para valores positivos de y . Eliminando el absoluto, existen soluciones para negativos y . Ahora depende de usted decidir qué significa y si esto es importante.
Por mi parte, me resulta más fácil hacerlo en Laplace:
donde BW es el ancho de banda, y es la frecuencia central. Para una raíz real, :
Lo cual es fácil de resolver. Si la raíz es imaginaria pura:
y si es complejo, :
Este es un poco más complicado de resolver, no imposible. Las raíces desnudas y complejas de la ecuación resultan como:
y resolviendo por complejo, obtienes cuatro posibilidades, de las cuales eliges, como antes, solo las frecuencias positivas:
Parece monstruoso, pero funciona.
LvW