heatmapmatrixplot
if true A = [45.24 16.81 -27.01 -9.614; -19.9 46.77 9.169 -28.33; -27.11 -9.381 55.07 16.98; 8.874 -28.66 -20.19 56.69]% code
end
x=0:0.1:3.14;y=sin(4*x).*sec(3*x);yy = zeros(1,length(x));yy(:) = 0;y1=-20:0.1:20;x1 = ones(1,length(y1));x1(:) = 0.5;area_1 = trapz(x(x<=0.5), y(x<=0.5)); % Calculate ‘Area #1’
xc = find((y .* circshift(y, [0 1])) < 0); % Define ‘y’ Zero-Crossings
idx_rng = xc(1)-1:xc(1); % Index Range Including Zero-Crossing
x0 = interp1(y(xc(1)-1:xc(1)), x(xc(1)-1:xc(1)), 0); % Find ‘x’ At Zero Crossing
area_2 = trapz([0.5 x0], [y(x == 0.5) 0]); % Calculate ‘Area #2’
plot(x,y,x,yy,'--',x1,y1,'--')hold onha1 = area(x(x<=0.5), y(x<=0.5), 'FaceColor','g');ha2 = area([0.5 x0], [y(x == 0.5) 0], 'FaceColor','r');hold offA1str = sprintf('Area 1 = %6.3f', area_1);A2str = sprintf('Area 2 = %6.3f', area_2);legend([ha1 ha2], A1str, A2str)
out = cat(1, matrix{:})
Best Answer