Matlab里抛物线如何绕一个点旋转，旋转后新的抛物线方程是什么

Matlab里抛物线如何绕一个点旋转，旋转后新的抛物线方程是什么

同学你好，在平面直角坐标系下，抛物线绕一个点旋转后，表达式是非常复杂的
设原来抛物线表达式为

y = a*x^2+b*x+c

通过这个推导之后可以反解出x和y用x'和y‘表示，然后将x和y代入原来的抛物线方程即可。

syms theta x y x1 y1 a b c
syms xc yc
eq = [x1;y1] - [xc;yc]+...
    [cos(theta) -sin(theta);
    sin(theta), cos(theta)]*[x-xc;y-yc];
[x,y] = solve(eq, x,y);
x = simplify(x)
y = simplify(y)

[x2,y2] = zhuanbian_(pi/2,0,0)这里运行出来的结果是（1/8，1/8），应该是（1/4，1/4）
下面是函数
function [x2,y2] = zhuanbian_(theta1,x3,y3)
syms theta x y x1 y1 a b c
syms xc yc
eq = [x1;y1] - [x;y]+...
[cos(theta) -sin(theta);
sin(theta), cos(theta)]*[x-0;y-1/4];
[x,y] = solve(eq, x,y);
x = simplify(x);
y = simplify(y);%来源

x2=subs(x,{theta,x1,y1},{theta1,x3,y3});%将theta带入
y2=subs(y,{theta,x1,y1},{theta1,x3,y3});%将theta带入