![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1414839/image.png)
Can anyone plot these function by using Matlab (R2023a version) or (R2023b version) and show me if there is any differenc between my figures and the Matlab (R2023a version)
4 次查看(过去 30 天)
显示 更早的评论
I ploted these two functions in the extended complexplain $\mathbb{C} \cup \{ \infty \}$
\begin{align*}
f_k(z)=(0.10+0.3i) + \sqrt{2}(0.3+0.0012i) z +(0.9-0.001i)z^{2}
\end{align*}
\begin{align*}
f_b(z)=(0.10-0.3i) z^{-1} + \sqrt{\frac{1}{2}}(0.3-0.0012i) z^{-2} +(0.9+0.001i)z^{-3}
\end{align*}
by using Matlab (R2018a version), I have the following figure,
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1414754/image.png)
Where I used the following code
t=0
f0=[0.10 + 0.3i ; 0.3 + 0.0012i ; 0.9 - 0.001i ];
H = [2 -3i 0;3i 2 0;0 0 1];
p=[ ( 0.9000 - 0.0010i) (0.4243 + 0.0017i) (0.1000 + 0.3000i) ]; %t=0
p1=[(0.9000 + 0.0010i) (0.2121 - 0.0008i) (0.1000 - 0.3000i) (0)] ; %t=0
f(:,:) = expm(i*H*t)*f0
re_z = -6.005:.01:6.005;
im_z= -6.005:.01:6.005;
[re_z,im_z] = meshgrid(re_z,im_z);
z = re_z + 1i*im_z;
f_of_z_result = polyval(p,z);
f_of_1_over_z_result = polyval(p1,1./z);
figure();
subplot(2,2,1)
surf(re_z,im_z,abs(f_of_z_result),'EdgeColor','none')
colorbar
title('|f_k(z)|')
xlabel('Z_R')
ylabel('Z_I')
zlim([0 15])
caxis([-6 6])
grid on
subplot(2,2,2)
surf(re_z,im_z,cos(angle(f_of_z_result)),'EdgeColor','none')
colorbar
title('cos(phase of f_k(z))')
xlabel('Z_R')
ylabel('Z_I')
zlim([-5 5])
caxis([-1 1])
subplot(2,2,3)
surf(re_z,im_z,abs(f_of_1_over_z_result),'EdgeColor','none')
colorbar
title('|f_b(1/z)|')
xlabel('Z_R')
ylabel('Z_I')
zlim([0 2])
caxis([-0.15 0.15])
subplot(2,2,4)
surf(re_z,im_z, cos(angle(f_of_1_over_z_result)),'EdgeColor','none'
colorbar
title(('phase of f(1/z)'))
title(('cos(phase of f_b(1/z))'))
xlabel('Z_R')
ylabel('Z_I')
caxis([-1 1])
zlim([-1 1])
grid on
My question is: can anyone plot these function by using Matlab (R2023a version) or (R2023b version) and show me if there is any difference between my figures and the Matlab (R2023a version) or (R2023b version) figures?
Many thanks
I appreciate any help
0 个评论
回答(2 个)
Richard Burnside
2023-6-19
I get this on 2023a:
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1414839/image.png)
2 个评论
Richard Burnside
2023-6-19
编辑:Richard Burnside
2023-6-19
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/1414874/image.png)
x-y view rotated...
The company where I work purchased a 2023a license for me recently. You could maybe use a trial version if you just have this one task to do.
The surf command in your 4th subplot needs a right paren at the end of the line.
Steven Lord
2023-6-19
Looking at the two pictures posted in this thread, I'm guessing the original poster had turned hold on (or had changed the default axes properties to the same effect) before plotting the data. If you do that surf will not override the view (changing it from 2-D to 3-D.)
figure
hold on % axes will stay in 2-D view
surf(peaks)
figure
axes; % not held
surf(peaks)
You could keep the axes held before plotting if after plotting you called view to explicitly change the view to 3-D. hold on just prevents automatic modification of axes properties; if you explicitly change those properties MATLAB will respect the changes.
figure
ax = axes;
hold on
[az, el] = view % 2-D default
surf(peaks)
view(3)
[az2, el2] = view
另请参阅
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!