代码搜索结果

wp=0.2*pi;Rp=1;ws=0.3*pi;As=15; [n,wn]=buttord(wp,ws,Rp,As) [b,a]=butter(n,wn); [db,mag,pha,grd,w]=freqz_m(b,a,0.5*pi); figure(3); subplot(2,2,1),plot(w/pi,mag);title('幅度'); subplot(2,2,2),plot(

program ural1192; const g=10; pi=3.1415926535; var v:longint; a:byte; k:real; begin readln(v,a,k); writeln(v*v*sin(a*pi/90)/g/(1-1/k):0:2); end.

% 编写建立自激过程微分方程的M—函数 % 将该M—函数定义为sh_ge_se_ex_ode(shunt_generator_self_excited_ode) function dydt=tongbufadianji(t,y) % 首先指定全局变量 %global a1 a2 a3 a4 Rf % 下面输入电机基本数据： r=2.9069E-03,Rfd=5.9013E-04

function [db,mag,pha,grd,w] = freqz_m(b,a); % Modified version of freqz subroutine % ------------------------------------ % [db,mag,pha,grd,w] = freqz_m(b,a); % db = Relative magnitude in dB comp

% Chapter 7 Example 7.25 % Highpass filter design using PM algorithm % ws = 0.6*pi; wp = 0.75*pi; Rp = 0.5; As = 50; delta1 = (10^(Rp/20)-1)/(10^(Rp/20)+1); delta2 = (1+delta1)*(10^(-As

% Chapter 7: Example 7.14 % Freq. Samp. Tech.: Lowpass,Naive method % M = 20; alpha = (M-1)/2; l = 0:M-1; wl = (2*pi/M)*l; Hrs = [1,1,1,zeros(1,15),1,1]; Hdr = [1,1,0,0]; wdl = [0,0.25

% Chapter 7: Example 7.18 % Freq. Samp. Tech.: Highpass, Optimum method T1 % ws=0.6pi, wp=0.8pi, Rp=1dB, As=50dB % M=33, T1 = 0.1095; T2 = 0.598; M = 33; alpha = (

% Chapter 7: Example 7.19 % Freq. Samp. Tech.: Differentiator % M = 33; alpha = (M-1)/2; Dw = 2*pi/M; l = 0:M-1; wl = Dw*l; k1 = 0:floor((M-1)/2); k2 = floor((M-1)/2)+1:M-1; Hrs = [j*

% Chapter 7: Example 7.15 % Freq. Samp. Tech.: Lowpass,Optimum method % % (a) T1 = 0.39 M = 40; alpha = (M-1)/2; l = 0:M-1; wl = (2*pi/M)*l; Hrs = [ones(1,5),0.3904,zeros(1,29),0.3904,

% Chapter 7: Example 7.22 % Error function plot and counting of extrema figure(1);subplot(1,1,1) h = [1,2,3,4,3,2,1]/15; M = length(h); n = 0:M-1; [Hr,w,a,L] = Hr_Type1(h); subplot(1,