graphing data using webread

Question

I W 2020-11-26

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/663058-graphing-data-using-webread

编辑： I W 2021-12-17

station_id = ["02DD016", "02HA014", "02HC030", "02HC024"]; titles = ["FRENCH RIVER AT PORTAGE DAM", "REDHILL CREEK AT HAMILTON",... "ETOBICOKE CREEK BELOW QUEEN ELIZABETH HIGHWAY","DON RIVER AT TODMORDEN"]; base_url = 'dd.weather.gc.ca/hydrometric/csv/'; province = 'ON'; frequency = 'daily'; file_type = 'csv'; for z = 1:4 my_url = strcat('https://', base_url, province, '/', frequency, '/', ... province, '_', station_id(z), ... '_', frequency, '_hydrometric.', file_type); gauge_data = webread(my_url); depth_data = gauge_data.WaterLevel_NiveauD_eau_m_.'; depth_data(find(isnan(depth_data)))=[]; y = depth_data; x = 1:1:length(y); avg_y_scalar = mean(y); avg_y_vector = avg_y_scalar *... ones(1,length(y)); std_y = std(y); y_plus = y + std_y; y_minus = y - std_y; y_minus_reverse = fliplr(y_minus); figure(1) subplot (2,2,z) % plot the data and average plot(x,y,x,avg_y_vector); % plot std deviation shape patch([1:1:length(y)... length(y):-1:1],... [y_plus y_minus_reverse],... 'b',... 'facealpha',0.05,...% fill colour 'edgecolor','r',... 'edgealpha',0.05) % edge colour % Legend, Title, axis labels. legend('original data',... 'average values',... 'standard deviation'); title(titles(z)); xlabel('Historical Daily Water Level [Day]'); ylabel('Water Level [m]'); end My first graph is supposed to look like this

but instead looks like this

and my other three graphs look fine so I do not know what the problem is

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

VBBV 2020-11-26

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/663058-graphing-data-using-webread#answer_556638

在 MATLAB Online 中打开

	figure(1)
	subplot (2,2,z)
	if z == 1
		% plot the data and average
		
		plot(x(1:100:end),y(1:100:end),x(1:100:end),avg_y_vector(1:100:end));
		%plot(x,y,x,avg_y_vector);
		patch([1:1:length(y)...
		length(y):-1:1],...
		[y_plus y_minus_reverse],...
		'b',...
		'facealpha',0.05,...% fill colour
		'edgecolor','r',...
		'edgealpha',0.004) % edge colour
	else
		plot(x,y,x,avg_y_vector);
		patch([1:1:length(y)...
		length(y):-1:1],...
		[y_plus y_minus_reverse],...
		'b',...
		'facealpha',0.05,...% fill colour
		'edgecolor','r',...
		'edgealpha',0.05) % edge colour
	% Legend, Title, axis labels.
	end