Matlab
文章平均质量分 66
zstarstone
研究方向:信息安全,大数据分析
展开
-
function num_lines = specobjid2deg(specobjid_filepath, deg_filepath, len_format
function num_lines = specobjid2deg(specobjid_filepath, deg_filepath, len_format)%Author: shizhixin%Email: szhixin@gmail.com%Date: Dec 11,2011%This function can change specobjid which is the form原创 2011-12-11 23:19:36 · 1388 阅读 · 0 评论 -
观测数据和模型数据根据参数变化的相关性度量
主要是先用三次样条插值对齐两个数据坐标,然后将模型数据进行三阶多项式乘积尽量拟合观测数据,最后求拟合的数据和观测数据之间的内积。内积最大者相关性最大。当然,也可以不通过拟合后内积,直接内积后除以两个向量的模。下面是MATLAB伪代码和源代码。% input: % Mx, My: Model spectrum % Sx, Sy: Observed spectrum% RV_原创 2012-06-06 16:31:47 · 1369 阅读 · 0 评论 -
找出两个将近十万多条记录的文件中的不同记录
% Author:shizhixin% Email:szhixin@gmail.com% Blog:http://blog.csdn.net/shizhixin% Date:2012-05-29% 实验记录:% 找出两个文件a,b中不同的记录,文件中记录的个数分别为m,n% 方法一:最直接的方法是通过两重循环,这样时间复杂度O(m*n)% 方法二:先进行快速排序,然后通过折半查找的原创 2012-05-29 17:19:57 · 2220 阅读 · 1 评论 -
对一维向量中值滤波
function flux_out = meanfilter(lamda, flux, w)% Author:shizhixin% Email:szhixin@gmail.com% Blog:http://blog.csdn.net/shizhixin% Date:2012-03-11% Function:对一维向量中值滤波% Note:输入lamda,flux必须维度一致,并且为一维原创 2012-03-11 14:18:46 · 8212 阅读 · 1 评论 -
通过中值滤波,归一化对比LAMOST与SDSS的光谱
% Author:shizhixin% Email:szhixin@gmail.com% Blog:http://blog.csdn.net/shizhixin% Date:2012-03-12% 通过退红移,中值滤波,然后归一化对光谱进行处理,对比LAMOST与SDSS的光谱。% 光谱文件为:DR8: 4.fits % LAMOST: 1.fits% 参数dered为是否退红移原创 2012-03-12 03:23:45 · 2726 阅读 · 0 评论 -
按照txt中指定的文件名,从src_path中拷贝文件到dest_path(copyfile_from_txt)
function num_processed = copyfile_from_txt(txt_filename,src_path,dest_path,add_num,ext_name)% Author:shizhixin% Email:szhixin@gmail.com% Blog:http://blog.csdn.net/shizhixin% Date:2012-03-02% Func原创 2012-03-02 15:16:35 · 3290 阅读 · 0 评论 -
从并排的两列数组中找到typename1或者为typename2的行号,并去重复
function type_pos = get_type_pos(array1, array2, typename1, typename2)%Author:shizhixin %Email:szhixin@gmail.com %Blog:http://blog.csdn.net/shizhixin %Date:2012-02-29%function:从并排的两列数组中找到原创 2012-02-29 16:58:02 · 939 阅读 · 0 评论 -
读取文件夹中的fit文件名,并且提取第11-12字符的数
function sp_no = get_fit_sp(path)%Author:shizhixin%Email:szhixin@gmail.com%Blog:http://blog.csdn.net/shizhixin%Date:2012-02-29%function:%读取文件夹中的fit文件名,并且提取第11-12字符的数% 如文件名ssy04_1_sp15.fit,原创 2012-02-29 16:59:18 · 1261 阅读 · 0 评论 -
被fclose忽悠了,try...catch...end及lasterr
今天主要把原来的程序整理了一遍,思路清晰多了,程序看起来也美观了。另外一个主要工作,是让程序能够自己一次跑完几千的数据,不会像上次那样,经常跑到几百条,这样那样的原因奔溃了。总结来看,主要有两个原因导致程序的奔溃:1、有些光谱数据,红移太大,退红移后的数据可能还未达到所需要的窗口数据,无法截取数据出来拟合,即所谓的:Index exceeds matrix dimensions.原创 2012-03-01 21:57:46 · 2224 阅读 · 0 评论 -
批量读取文件夹下的DR8光谱数据,并完成显示保存图片以及统计Z和Obj的分布情况
function [Z Obj] = read_dr8_dir(path,option,show_n, save_img, hist_info)%Author:shizhixin%Email:szhixin@gmail.com%Blog:http://blog.csdn.net/shizhixin%函数功能:批量读取文件夹下的DR8光谱数据,%并完成显示保存图片以及统计Z和Ob原创 2012-02-17 18:43:21 · 1677 阅读 · 0 评论 -
matlab重命名文件movefile
function rename_dir(path, option)%Author:shizhixin%Email:szhixin@gmail.com%Blog:http://blog.csdn.net/shizhixin%Date:2012-02-15%option = 'LAMOST','DR8','SPECIAL'%LAMOST:类似spec_release.php@path=%原创 2012-02-15 10:15:52 · 8557 阅读 · 3 评论 -
信号峰拟合的MATLAB程序,包括高斯拟合,多高斯拟合等多种类型
今天准备弄双高斯拟合,看到一个信号峰拟合的MATLAB版本的程序,大体看了一下,很不错,先MARK一下,以后再详细研究。 http://terpconnect.umd.edu/~toh/spectrum/CurveFittingC.htmlhttp://terpconnect.umd.edu/~toh/spectrum/InteractivePeakFitter.htm 下面是其转载 2011-12-28 16:37:44 · 27492 阅读 · 8 评论 -
Get the window data(lamda,flux, true) from a center point and the window width
Function:function [lamda_win, flux_win] = show_window(lamda_cent, w_win, lamda, flux, flag_plot)%Date:Dec 27,2011%Author:shizhixin%Email:szhixin@gmail.com%Blog:http://blog.csdn.net/shizhixin%原创 2011-12-27 15:58:49 · 1033 阅读 · 0 评论 -
MATLAB中小技巧集锦
1、遍历文件夹path = 'photos/';files = dir([path '*.jpg']);n = length(files);for i = 1:n file_path = [path files(i).name]; file_name = files(i).nameendoutput:file_name =DSC_0253.jpgfile_nam原创 2011-12-19 15:48:56 · 1675 阅读 · 0 评论 -
FIT文件格式解析及MATLAB读取程序
FIT文件格式解析(http://blog.csdn.net/shizhixin,Email: zxshi@nlpr.ia.ac.cn,2011年12月) 摘要:本文档主要介绍了用UltraEdit查看FIT文件及简要的说明了FIT文件的数据存储方式,并附上了FIT文件数据读取及保存为TXT的MATLAB程序。此外,通过读取lamda的相关信息,文中还显示了5条FIT文件中读取出来的光谱。原创 2011-12-14 16:58:31 · 16250 阅读 · 2 评论 -
按指定的列col排序数组A,并保持数组行相对不变
function B=sort_array(A, col)%按指定的列col排序数组A,并保持数组行相对不变% Note:此函数仅适合N*2的数组排序,col取 1 or 2,否则B=0% ie:% A =% 1 4% 2 3% 1 2% 4 6% 2 3% B = sort_array(A, 2)% B =% 1 2% 2 3% 2 3% 1 4% 4 6% B原创 2011-06-24 10:31:00 · 772 阅读 · 0 评论 -
获取当前鼠标按键的坐标值,左键输出当前点,右键则删除上一个按键值,中间键退出
clear,close all, clcx = 1:0.1:10;y = 20*sin(pi/2*x);figureplot(x,y);points = get_mouse_point()-----------------------------------------------------function points = get_mouse_point()% Au原创 2013-03-08 13:50:49 · 1986 阅读 · 0 评论