用matlab与r语言画的两幅图,个人感觉:r语言画图可以设置的项目比matlab多一些,更灵活一些(例如可以分别设置x,y轴的tick的长度),但是matlab画出的图以后可以进行编辑,另外matlab画图的参数比R语言也好理解一些.下面是画图的代码.
- rm(list=ls())
- graphics.off()
- setwd("e:/")
- yi <- read.table("yi.txt")
- a <- yi
- b<-c("北京","天津","石家庄","唐山","秦皇岛","太原","呼和浩特","包头"
- ,"沈阳","大连","丹东","锦州","长春","吉林","哈尔滨"
- ,"牡丹江","上海","南京","无锡","徐州","扬州","杭州"
- ,"宁波","温州","金华","合肥","蚌埠","安庆","福州"
- ,"厦门","泉州","南昌","九江","赣州","济南","青岛"
- ,"烟台","济宁","郑州","洛阳","平顶山","武汉","宜昌"
- ,"襄阳","长沙","岳阳","常德","广州","韶关","深圳"
- ,"湛江","惠州","南宁","桂林","北海","海口","三亚"
- ,"重庆","成都","泸州","南充","贵阳","遵义","昆明"
- ,"大理","西安","兰州","西宁","银川","乌鲁木齐")
- numofmonth<-dim(a)[1]
- for (ii in 1:1){
- windows()
- jpeg(filename = paste(b[ii],".jpg",sep=""),
- width = 3200, height = 2400, units = "px", pointsize = 12,
- quality = 75,res = 300, family = "", restoreConsole = TRUE,
- type = c("windows", "cairo"))
- par(mar=c(4,4,2,0.5)+0.1,xpd=TRUE)
- plot(1:numofmonth,a[,ii],type = "o",col = rgb(runif(1,0,1),runif(1,0,1),runif(1,0,1)), lwd = 1,pch=16,cex.lab=1,
- main = b[ii],xlab="时间",ylab="相对房价",xlim=c(2,numofmonth-1),xaxt="n",yaxt="n",
- cex.main=1,cex=1.2,cex.axis=1,cex.lab=1,bg="yellow",frame=TRUE,fg="black",font.lab=2,font=2)
- axis(1,at=seq(3,numofmonth,by=6),labels=c("11/03","11/09","12/03","12/09","13/03","13/09","14/03","14/09","15/03"),
- tick=TRUE,col.axis="black",col = "black", lwd=1,lty=1,tck=-0.01,font=2,
- padj=-0.9,hadj=0.5,mgp=c(1,1,0) )
- axis(2,col.axis="black",col = "black", lwd=1,las=2,tck=-0.01,font=2,
- xaxs="i",padj=0.5,hadj=0.8,mgp=c(1,1,0) )
- n<-30
- bl<-1/3
- get1<-par()
- qujian<-0.99*(get1$usr[4]-get1$usr[3])
- len0<-qujian/(n+n*bl-bl)
- len1<-len0*bl
- cc<-seq(3,numofmonth,by=3)
- for (jj in cc){
- for (kk in 1:n){
- lines(jj*rep(1,2),c(get1$usr[3]+(len0+len1)*(kk-1)+1/198*qujian,get1$usr[3]+(len0+len1)*(kk-1)+len0)+1/198*qujian,col=rgb(250/255,128/255,114/255))
- }
- }
- dev.off()
- }
matlab:
- clear all
- close all
- load yi.txt
- a=yi;
- b={'北京' '天津' '石家庄' '唐山' '秦皇岛' '太原' '呼和浩特' '包头'...
- '沈阳' '大连' '丹东' '锦州' '长春' '吉林' '哈尔滨'...
- '牡丹江' '上海' '南京' '无锡' '徐州' '扬州' '杭州'...
- '宁波' '温州' '金华' '合肥' '蚌埠' '安庆' '福州'...
- '厦门' '泉州' '南昌' '九江' '赣州' '济南' '青岛'...
- '烟台' '济宁' '郑州' '洛阳' '平顶山' '武汉' '宜昌'...
- '襄阳' '长沙' '岳阳' '常德' '广州' '韶关' '深圳'...
- '湛江' '惠州' '南宁' '桂林' '北海' '海口' '三亚'...
- '重庆' '成都' '泸州' '南充' '贵阳' '遵义' '昆明'...
- '大理' '西安' '兰州' '西宁' '银川' '乌鲁木齐'};
- numofmonth=size(a,1);
- for ii=1:1
- figure
- set(gcf,'visible','on')
- set(gcf,'color',[0.96 0.96 0.96],'position',[300 100 750 500])
- ys=rand(1,3);
- plot(a(:,ii),'-r.','markersize',15,'markerfacecolor',ys,...
- 'markeredgecolor',ys,'color',ys,'linewidth',1)
- hold on;
- set(gca,'color',[0.96,0.96,0.96],'position',[0.08,0.08,0.9,0.86],'fontweight','bold')
- xlabel('时间','fontweight','bold')
- ylabel('相对房价','fontweight','bold')
- title(b{ii},'fontweight','bold')
- set(gca,'xtick',3:6:numofmonth,'xticklabel',{'11/03','11/09','12/03','12/09','13/03','13/09', '14/03','14/09','15/03'},'tickdir','in')
- ylim=get(gca,'ylim');
- n=30;
- bl=1/3;
- qujian=0.99*diff(ylim);
- len0=qujian/(n+n*bl-bl);
- len1=len0*bl;
- for jj=3:3:numofmonth
- for kk=1:n
- plot(jj*ones(2,1),[ylim(1)+(len0+len1)*(kk-1)+1/198*qujian,ylim(1)+(len0+len1)*(kk-1)+len0+1/198*qujian], '-','linewidth',0.8,'color',[250 128 114]/255);
- end
- end
- xlim([0 numofmonth+1])
- set(gcf,'InvertHardcopy','off');
- set(gcf,'PaperPositionMode','auto')
- print(gcf,b{ii},'-djpeg','-r300')
- end
数据文件在附件中.欢迎各位坛友讨论. |
熟悉论坛请点击新手指南
京公网安备 11010802022788号
论坛法律顾问：王进律师
知识产权保护声明
免责及隐私声明
