fart

芬兰作者提出的FastICA算法，用以解决盲源分离问题，我编写了一个主程序，调用其已经编写好的功能函数，可以查看仿真结果。(Finland FastICA algorithm proposed by the authors to solve the blind source separation problem, I wrote a main program, calling its already written performance function, you can view the simulation results.)

说明：  最大似然估计频偏，符号定时误差的MATLAB仿真(MATLAB simulation maximum likelihood estimation frequency offset)

一个计算混沌时间序列相关维数的matlab程序(A calculation of chaotic time series correlation dimension of the matlab program)

matlab 代码为图像的分割，执行效果可以看程序的执行(matlab code for image segmentation, the performance of the implementation of the procedures can be seen)

Good application for images compression.with neural networks

说明：  张贤达《现代信号处理（第二版）》习题与解答(Xian-up &quot Modern Signal Processing (Second Edition)&quot problem sets and answers)

一个GPS课程的全部代码，含有数据！很实用(MATLAB GPS SCRIPTS-- Charles Meertens and Doug Miller, U. Utah 3Jan95 These scripts were developed for/during a class in GPS and Active Tectonics Department of Geology and Geophysics, U. Utah, Spring, 1994 EXAMPLE RINEX Files (the only ones guaranteed to work this release) CUR10911.94N CUR10911.94O PAR10910.94N PAR10910.94O Main Program scripts: gpshora.m -- function to calculate Modified Julian Day, Homework#1 HW2.m -- Homework#2 read RINEX observation file and plot C/A code vrs. time HW3_plt.m-- Homework#3 read RINEX navigation file, plot satellite positions HW3_3dplt.m-- Same as HW3_plt.m, but 3d position plots HW4.m-- Homework#4 calculate and plot station position as function of time Functions called: calc_nav.m-- calculate satellite position calc_pos1.m-- calculate station position constants.m-- GPS constants coords.m-- convert x,y,z coordinates to lat, long, height date2j.m-- calculate MDJ eccenr.m-- calculate "E" anomally read_nav.m-- Read RINEX Nav file r)

决策树的ID3算法的源代码，最后是其实验数据，运行很简单！(ID3 decision tree algorithm source code, and finally the experimental data, the run is very simple!)

ber of qpsk modulation

empirical formula with kaiser clc clear all fs=1000 fc=250 df=50 r=0.001 f=fc/fs dw=2*pi*(df/fs) a=-20*log(r) n=floor(((a-8)/(2.285*dw))+1) if a>50 b=0.1102*(a-8.7) elseif a>=21 && a<=50 b=0.5842*((a-21)^0.4)+0.07886*(a-21) elseif a<21 b=0.0 end w=kaiser(n,b) for i=1:n if i~=(n-1)/2 hd(i)= (2*f*sin((i-((n-1)/2))*2*pi*f))/((i-((n-1)/2))*2*pi*f) elseif i==(n-1)/2 hd(i)=2*f end end for j=1:n h(j)=w(j)*hd(j) end subplot(3,1,1), plot(w) subplot(3,1,2), plot(h) subplot(3,1,3), plot(h,n)