Without really trying to look at your maths, let me suggest that you will have less issues if you treat your 3D objects in a standard way and work backwards from there.
The normal way to define a plane in 3D is to specify 3 vectors: centre, up, and right. The last two define the plane's axis system and are usually unit vectors (to obtain the forward-vector of the axis you take the cross-product of up and right).
These can go into a 3x3 transformation matrix as follows:
T = [right;forward;up];
Normally, you use a 4x4 (or sometimes 4x3) transformation matrix, but let's stick with the 3x3.
You can define your rotation (rotation about Z-axis) as you already have, but in matrix form:
rotZ = @(x) [cosd(x) -sind(x) 0; sind(x) cosd(x) 0; 0 0 1];
And a the tilt is rotation about the X-axis:
rotX = @(x) [1 0 0; 0 cosd(x) -sind(x); 0 sind(x) cosd(x)];
[by the way, your rotations all seem to be clockwise, whereas anticlockwise is more standard -- I've stuck with clockwise here]
Now your initial plane transformation is the unit. Z is up, and the X- and Y-axes are on the plane.
T = diag([1 1 1]); % Unit transform
You're ready to rotate and tilt your plane:
% First tilt by 90 degrees, then rotate by 45 degrees.
S = rotZ(45) * rotX(90) * T;
You could pull your axes out like so:
X = S(1,:);
Y = S(2,:);
Z = S(3,:);
... and then linspace a mesh by projecting onto these new axes (remembering they should be unit-length)...
But remember that an arbitrary transformation will not define a mesh with a 3 simple vectors of x,y, and z co-ordinates. Every single x,y,z vector in a mesh could be unique. So that would be the wrong thing to do.
Instead you need to rotate the whole mesh. First, generate the original mesh:
x = linspace(-xsize,xsize,N) / 2;
y = linspace(-ysize,ysize,N) / 2;
[mx,my,mz] = meshgrid(x,y,0);
M = [mx(:) my(:) mz(:)];
M contains co-ords of untransformed mesh points organised in rows.
Mtrans = M * S;
Finally you add your centre point back to the transformed mesh.
Mfinal = Mtrans + repmat(cp, size(M,1), 1);
Mfinal now contains your whole mesh as a collection of points, rotated and translated however you like. It's up to you to work out what to do with it. =)
