clear; close all; s = tf('s'); xi = 0.1; e = sqrt(10.^(xi/10)-1); tmp = sqrt(1+e^(-2))+e^(-1); theta = 0:0.1:2*pi; x = cos(theta); y = sin(theta); for n=1:9 D = 1; for k=1:n sinhv = 1/2*(tmp^(1/n)-tmp^(-1/n)); coshv = 1/2*(tmp^(1/n)+tmp^(-…

ξ=0.5dB, 1≤n≤9 clear; close all; s = tf('s'); xi = 0.1; e = sqrt(10.^(xi/10)-1); tmp = sqrt(1+e^(-2))+e^(-1); theta = 0:0.1:2*pi; x = cos(theta); y = sin(theta); for n=1:9 D = 1; for k=1:n sinhv = 1/2*(tmp^(1/n)-tmp^(-1/n)); coshv = 1/2*(t…

ξ=0.5d, 1≤n≤9 clear; close all; s = tf('s'); xi = 0.1; e = sqrt(10.^(xi/10)-1); tmp = sqrt(1+e^(-2))+e^(-1); theta = 0:0.1:2*pi; x = cos(theta); y = sin(theta); for n=1:9 D = 1; for k=1:n sinhv = 1/2*(tmp^(1/n)-tmp^(-1/n)); coshv = 1/2*(tm…

修正したlittlebirdy /usr/local/bin/littlebirdy #!/bin/bash set -e # Copyright (c) 2012 Jake Petroules. All rights reserved. # # Redistribution and use in source and binary forms, with or without # modification, are permitted provided that …

clear; close all; s = tf('s'); xis = [0.1 0.5 1]; es = sqrt(10.^(xis/10)-1); e = es(1); tmp = sqrt(1+e^(-2))+e^(-1); for n = 1:10 D = 1; for k=1:n sinhv = 1/2*(tmp^(1/n)-tmp^(-1/n)); coshv = 1/2*(tmp^(1/n)+tmp^(-1/n)); sigma_k = -sin((2*k-…

clear; close all; syms w; chebyshev = @(c1,c0) 2*w*c1 - c0; c0 = 1; c1 = w; cs = [c0 c1]; for n=2:10 cn_prev = cs(n-1); cn = cs(n); cn_next = chebyshev(cn, cn_prev); cs(end+1) = cn_next; end figure(); hold on; for n=1:4 cnw = cs(n+1); x = …