I'm using the following code, but not getting a surface plot (rather, only the wiring underneath). Do you know why this might be? Thanks in advance for your help.
rotation = RotationToRotationStructure_Data(:,1); exactK = ScalingToRotationStructure_Data(:,1); scale_tx = ScalingToRotationStructure_Data(:,2); scale_ty = ScalingToRotationStructure_Data(:,3); rot_tx = RotationToRotationStructure_Data(:,2); rot_ty = RotationToRotationStructure_Data(:,3); scale_ReportX = ScalingToRotationStructure_Data(:,4); scale_ReportY = ScalingToRotationStructure_Data(:,5); rot_ReportX = RotationToRotationStructure_Data(:,4); rot_ReportY = RotationToRotationStructure_Data(:,5); nUniquePert = numel(unique(exactK(1:240,1))); %expected values
thetaHat = zeros(2400,1); rot_scaleHat = thetaHat; for j = 1:length(thetaHat); scale_expected_ReportX(j) = exactK(j)*scale_tx(j) + ... exactK(j)*scale_ty(j); scale_expected_ReportY(j) = -exactK(j)*scale_tx(j) + ... exactK(j)*scale_ty(j); rot_expected_ReportX(j) = 1*cos(rotation(j))*rot_tx(j) + ... 1*sin(rotation(j))*rot_ty(j); rot_expected_ReportY(j) = -1*sin(rotation(j))*rot_tx(j) + ... 1*cos(rotation(j))*rot_ty(j); scale_goodfit(j) = rmse(scale_expected_ReportX(j) - scale_ReportX(j), scale_expected_ReportY(j) - scale_ReportY(j)); rot_goodfit(j) = rmse(rot_expected_ReportX(j) - rot_ReportX(j), rot_expected_ReportY(j) - rot_ReportY(j)); end x = repmat(thetaHat, 1, 2); y = repmat(exactK, 1, 2); z = repmat(scale_goodfit',1,2); colormap hsv; h = surf(x, y, z);
Best Answer