自动化下载解压缩基因组序列脚本-mm9

% 每次下载的基因组序列都很心累, 搞个自动化的
clear all;
clc;
close all;

url_base = 'http://hgdownload.soe.ucsc.edu/goldenPath/mm9/chromosomes/';

chr_list = {};
for i = 1:19
　　chr_list = [chr_list, {['chr', num2str(i)]}];
end
chr_list = [chr_list, {'chrX', 'chrY', 'chrM'}];
% 开始构造需要的url, 具体实例为 'http://hgdownload.soe.ucsc.edu/goldenPath/mm9/chromosomes/chr1.fa.gz'
for i = 1:length(chr_list)
　　if ~exist([chr_list{i}, '.fa'])
　　　　url_full = [url_base, chr_list{i}, '.fa.gz'];
　　　　while true
　　　　　　try
　　　　　　　　urlwrite(url_full, [chr_list{i}, '.fa.gz']);
　　　　　　　　gunzip([chr_list{i}, '.fa.gz']);
　　　　　　　　break;
　　　　　　catch
　　　　　　　　fprintf('try again.....\n');
　　　　　　end
　　　　end
　　end
　　fprintf([chr_list{i}, '.fa.gz\n']);
end