File:Kbd-window.svg
Jump to navigation
Jump to search
Size of this PNG preview of this SVG file: 700 × 525 pixels. Other resolutions: 320 × 240 pixels | 640 × 480 pixels | 1,024 × 768 pixels | 1,280 × 960 pixels | 2,560 × 1,920 pixels.
Original file (SVG file, nominally 700 × 525 pixels, file size: 82 KB)
Summary
| Description |
English: Kaiser-Bessel derived (KBD) window |
| Date | |
| Source | Own work |
| Author | RetoGalli |
| SVG development |
Licensing
RetoGalli, the copyright holder of this work, hereby publishes it under the following license:
 
|
This file is made available under the Creative Commons CC0 1.0 Universal Public Domain Dedication. |
| The person who associated a work with this deed has dedicated the work to the public domain by waiving all of their rights to the work worldwide under copyright law, including all related and neighboring rights, to the extent allowed by law. You can copy, modify, distribute and perform the work, even for commercial purposes, all without asking permission.
|
Source code
This graphic was created by the following Octave script:
pkg load signal
graphics_toolkit gnuplot
clear all; close all; clc
betas = [2 8 24 100];
alphas = betas/pi
colors = [0 0 1; 0 1 0; 0 1 1; 1 0 1]; %blue green cyan magenta
hfig = figure;
grid on; hold on; box on
N = 128; % Relative curve shapes are not sensitive to this number
set(gca, "XTick",[0 : 1/8 : 1]*2*N,...
"XTickLabel",[" 0"; " "; " "; " "; " "; " "; " "; " "; "2N-1 "])
set(gca, "YTick",[0 0.2 0.4 0.6 0.707 0.8 1.0] )
xlim([0 2*N])
ylim([0 1.05])
for k= 1:length(betas)
w = besseli(0,betas(k)*sqrt(1-(2*(0:N)/N -1).^2))/besseli(0,betas(k));
d = zeros(1,2*N); % length 2N
for n=0:N-1
d(1+n) = sum(w(1+(0:n)));
end
for n=N:2*N-1
d(1+n) = sum(w(1+(0:2*N-1-n)));
end
d = sqrt(d/sum(w(1+(0:N))));
plot(0:2*N-1, d, "color", colors(k,:), "linewidth", 2)
end
title("Kaiser-Bessel-derived window functions")
ylabel(" ") % create left margin
text(8, .96, 'd_n', "fontsize", 14)
text(100, -0.05, '\leftarrow n \rightarrow')
%h = legend('\pi\alpha= 2; \alpha=0.64', '\pi\alpha= 8; \alpha=2.55',...
% '\pi\alpha= 24; \alpha=7.64', '\pi\alpha=100; \alpha=31.8');
%But let's do it a less "manual" way:
h = legend(['\pi\alpha= ' num2str(betas(1),'%3i') '; \alpha=' num2str(betas(1)/pi,'%5.2f')],...
['\pi\alpha= ' num2str(betas(2),'%3i') '; \alpha=' num2str(betas(2)/pi,'%5.2f')],...
['\pi\alpha= ' num2str(betas(3),'%3i') '; \alpha=' num2str(betas(3)/pi,'%5.2f')],...
['\pi\alpha=' num2str(betas(4),'%3i') '; \alpha=' num2str(betas(4)/pi,'%5.2f')],...
"location","south");
legend boxoff
set(h, "fontsize",10);
% The following print() converts plain-text Greek characters in text() strings into Symbol font,
% but legend() crashes it. The set() succeeds, but generates warnings, some of which are
% diabled by warning().
% print(hfig,"-dsvg","-color",'C:\Users\BobK\Kbd-window.svg')
warning("off", "Octave:missing-glyph");
set(h, "fontname","Symbol");
File history
Click on a date/time to view the file as it appeared at that time.
| Date/Time | Thumbnail | Dimensions | User | Comment | |
|---|---|---|---|---|---|
| current | 17:36, 24 March 2019 | | 700 × 525 (82 KB) | wikimediacommons>Bob K | replace M with N |
File usage
The following page uses this file:
