PI-LNK409 源代码

L0=0;L(1)=0;L(2)=3.72;L(3)=1.62;L(4)=1.05;x0=5;xf=3.5;z0=-2;zf=-2;%挖掘机参数 af=atan((z0-zf)./(x0-xf));%判断所给直线的角度，或者直接给出直线的角度，则更加的方便！…

function phi = dAngle(a,b) d = a-b; phi = mod(d,2*pi); if phi > pi phi = 2*pi-phi; end

%matlab计算程序如下: % Poincare_section[绘制庞加莱截面图] [t,x]=ode45(@Poincare,[0,2800],[0,1.5,0]); x(:,2)=mod(x(:,2),2*pi)-pi; phi0=…

[t,y]=nlbound('c7fun4','c7fun3',[0,pi/2],[-1,1],1e-8); plot(t,y); set(gca,'xlim',[0,pi/2]); y0=tan(t-pi/4); norm(y(:,1)-y0) …

%《数字信号处理教程——MATLAB释义与实现》第七章例7.4.1b程序hc741b % 取不同样本数时的幅频特性 % 电子工业出版社出版 陈怀琛编著 2004年9月 % N=input('N= (N必须为奇数）');wc=input('wc=');…

%《数字信号处理教程——MATLAB释义与实现》第七章例7.4.1a程序hc741a % 理想频率响应的频率样本设定 % 电子工业出版社出版 陈怀琛编著 2004年9月 % N=input('N= (N必须为奇数）');wc=input('wc=')…

%《数字信号处理教程——MATLAB释义与实现》第三章演示程序fgp314 % 数字频率与模拟频率的关系 % 电子工业出版社出版 陈怀琛编著 2004年9月 % T=1; Omega=[-2.5*pi:0.01:2.5*pi]/T; % 设定…

clf clc clear all a=0.1; A=0.2:0.1:1; s=size(A); N1=50; e=8.854e-12; mu=pi*4e-7; for K=1:s(2) AA=A(K); dal_n=2*pi…

function angles = normalize(angles) % Normalizes angles to interval (-pi,pi] %angles = pi - mod(pi-angles, 2*pi); % Same resu…

a=0;b=2*pi;n=1000;h=(b-a)/n; x=a:h:b; f=cos(x-pi/6).*sin(x+pi/6); for i=1:n s(i)=(f(i)+f(i+1)).*h/2; end s=sum(s)