代码搜索结果

%图形文字标示命令的使用 clear close all clc t=[0:pi/20:5*pi]; plot(t,sin(t),'r:*') axis([0 5*pi -1.5 1.5 ]) %给x轴和y轴命名 xlabel('t(deg)') ylabel('magnitude') %给图形加标题 title('sine wave from zero to 5\pi')

%plot绘图命令的使用 clc %清屏命令 clear %清除工作空间中所有变量 %定义时间范围 t=[0:pi/20:9*pi]; y1=sin(t); y2=cos(t); plot(t,y1,t,y2)

real*8 a, pai, pi_ a=1.d-307; print *, a, 1./a; print *, tiny(a), huge(a) b=1.e-39 ; print *, b, 1./b; print *, tiny(b), huge(b) pi =3.14159265358979_8 ; print *, pi ! pi 无双精度定义，仍为单精度。 a =3.

! ------求圆錐体积和表面积------ PROGRAM Example_1_2 REAL :: pi, r, h, v, s pi = 3.141593 PRINT *,'Input radius r and height h ?' READ *, r, h v = pi*h*r**2/3.0 s = pi*r*(r + sqrt(r**2 + h**2))

function [A0,A,B]=FZZ(func,T, n) syms t; func = subs(sym(func), findsym(sym(func)),sym('t')); A0=int(sym(func),t,-T/2,T/2)/T; for(k=1:n) A(k)=int(func*cos(2*pi*k*t/T), t,-T/2,T/2)*2/T; A

%例程3-4 产生锯齿波形 Fs=1000; t=0:1/Fs:3; x=sawtooth(2*pi*t); y=sawtooth(2*pi*t,0.5); z=sawtooth(2*pi*t,1); subplot(3,1,1); plot(t,x); subplot(3,1,2); plot(t,y); subplot(3,1,3); plo

%例程3-12 求解单位冲击响应的幅频 N=20; a=[0.4 0.3 -1]; b=[1 -0.2 0.6 0.3]; x=[1 zeros(1,N-1)]; k=0:N-1; y=filter(a,b,x); Fs=1000; w=0:pi/Fs:pi; h=freqz(a,b,w); subplot(2,1,1); stem(k,y);

%例程3-1 计算序列抽样和与序列能量 t=(0:0.1:10)'; x=sin(2*pi*t)+2*cos(pi*t); y=x+randn(size(t)); s=sum(x(1:100)); E=sum(abs(y).^2); subplot(2,2,1); plot(t,x); title('X'); subplot(2,2,2); plot(t,y); ti

%例程4－25 基于频率抽样方法2实现Ⅲ型FIR数字带通滤波器 clear all; N=41; T1=0.1095; alpha= (N-1)/2; l=0:N-1; wl= (2*pi/N)*l; Hrs=[zeros(1,6),T1,ones(1,7),T1,zeros(1,11),T1,ones(1,7),T1,zeros(1,6)];

%例程4－24 基于频率抽样方法1实现Ⅰ型FIR数字高通滤波器 clear all; N=61; T1=0.1095; T2=0.598; alpha= (N-1)/2; l=0:N-1; wl= (2*pi/N)*l; Hrs=[zeros(1,22),T1,T2,ones(1,14),T2,T1,zeros(1,21)]; %理想振幅采样响应 Hdr